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ABSTRACT 


Kolmogorov's  principles  provide  a  basis  for  th<^ 
treatment  of  Eulerian  space-time  correlations  for  turbulence 
in  the  universal  range  by  explicit  separation  of  the  kinematic 
effect  of  convection  by  the  large  eddies.  With  reference  to 
unsheared  homogeneous  turbulence,  the  usual  similarity  forms 
are  assumed  for  intrinsic  velocity  correlations  in  a  local  co¬ 
moving  frame.  The  local-convection  approximation,  neglecting 
dispersion  in  this  frame  and  hence  relating  space-time  correla¬ 
tions  to  pure  spatial  correlations,  is  indicated  to  be  a  useful 
one  in  the  inertial  and  viscous  subranges.  For  an  isotropic 
normal  velocity  distribution,  the  structure  functions  of  fluctua¬ 
ting  velocity  are  computed  and  found  to  be  nearly  space-time 
isotropic  in  the  former  subrange  and  exactly  so  in  the  latter. 

The  wavenumber- frequency  spectrum  of  energy  in  the  mean  rest 
frame,  E^(k,ju),  however,  in  the  regime  of  large  uu/vQk,  where  vq 
denotes  rms  velocity,  is  not  given  correctly  by  the  local- 
convection  approximation,  but  essentially  involves  disper¬ 
sion.  Taylor's  hypothesis  relating  wavenumber  (k^)  spectra 
in  the  mean  rest  frame  to  frequency  spectra  in  a  measurement 

1 . 


frame  with  velocity  -uq  is  found,  in  the  local -convection 
approximation  for  the  inertial  subrange,  to  be  exactly  valid 
for  all  vq/uo  provided  a  particular  effective  convection  velocity 
different  from  uq  is  assumed;  for  sufficiently  large  kjL,  where 
L  is  the  scale  of  energy -containing  eddies,  the  dispersive 
correction  to  this  result  is  negligible.  A  plausible  explicit 
form  is  proposed  for  the  intrinsic  energy  spectrum  in  the  co- 
moving  frame,  and  consequent  corrections  to  results  of  the 
local-convection  approximation  for  space-time  correlations 
and  spectra  are  computed.  This  dispersive  correction  to  correla¬ 
tions  is  appreciable  at  Reynolds  numbers  of  typical  grid- 
turbulence  experiments.  An  extension  of  the  basic  separation 
of  convection  and  the  local -convection  approximation  to  shear 
flow  is  suggested. 

From  kinematic  and  similarity  arguments  for  the 
inertial  subrange,  inferences  are  also  made  concerning  pressure 
spectra  in  a  measurement  frame  with  velocity  (-uq,0,0)  relative 
to  the  unsheared  flow.  A  similarity  form  results  for  the 

A  _  _ 

wavenumber- frequency  spectrum  P(K,uj)[K  =  (k^k^)]  in  the  vicinity 

of  the  convective  ridge  ( luu  -  kjUQ|  £  v0K) .  For  vQ/uo«l,  the 

functional  form  is  then  obtained  for  the  point  frequency  spectrum 
A  -7  /3 

P(<d)(«uu  );  the  convective  contribution  to  the  spectrum  Q(<ju)  of 

average  pressure  on  a  moving  circular  area  of  radius  Rq  for 

_3 

uuRo/uo»l  is  smaller  by  a  factor  0.8(u)Rq/uo)  .  On  assumption 

of  quasinormality  of  the  velocity  distribution  and  use  of  the 
non-dispersive  approximation  of  space-time  isotropy  for  the 

A 

longitudinal  correlation,  the  spectrum  P(R,uu)  in  the  vicinity 
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of  Che  convective  ridge  is  determined  explicitly,  as  well  as 
the  cross-spectral  density  and  space- tine  correlation  of 

A 

pressure.  From  quasinormality,  P(R,u>)  is  estimated  also  in 

the  low -wavenumber  domain.  In  the  limit  ujR0/v0»1  the 

ratio  of  the  low -wavenumber  to  the  me  an -convective  contribution 

*  10 /3 

to  the  moving-area  spectrum  Q(tu)  is  ~(vQ/uo)  '  (tuRQ/uo). 
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1. 


INTRODUCTION 


In  previous  descriptive  and  phenomenological  treat¬ 
ments  of  turbulence,  Eulerian  space-time  correlations  of 
fluctuating  velocity  have  received  limited  attention  relative 
to  ihat  accorded  spatial  correlations  at  a  fixed  time.  Yet 
time  correlations  have  been  extensively  measured,  are  amen¬ 
able  to  familiar  kinematic  and  similarity  arguments,  and 
have  essential  importance  in  a  number  of  applications. 

Eulerian  time  correlations,  i.e.,  in  the  general 
sense,  correlations  of  quantities  measured  in  a  coordinate 
frame  whose  uniform  motion  is  defined  independently  of  the 
decorrelating  fluid  motions  in  question,  are  influenced 
by  convection.  Convection,  where  it  has  meaning,  is  a  kin¬ 
ematic  effect  and  c'4n  he  treated  trivially  and  independently 
of  the  turbulent  oynamics.  Moreover,  when  the  large-scale 
convective  effects  are  separated  out,  it  becomes  possible 
to  invoke  similarity  considerations  to  treat  the  effect  of 
the  dispersive  residual  fluid  motion. 

Accordingly,  some  principal  objectives  of  the 
present  paper  are  to  elucidate  the  role  and  consequences  of 
convection,  assess  the  validity  and  explore  the  consequences 
and  limitations  of  the  local-convection  (non-dispersive) 
approximation  relating  space-time  correlations  to  purely 
spatial  correlations,  and  to  propose  and  apply  a  dispersive 
generalization  based  on  similarity  in  the  inertial  subrange. 
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These  questions  are  pursued  in  the  spirit  of  Kolmogorov' s 
principles;  no  effort  is  made  to  deal  with  the  explicit 
dynamics  of  turbulent  flow.  Further  verification  or  mod¬ 
ification  may  presumably  be  expected  from  currently  devel¬ 
oping  treatments  of  turbulent  dynamics.  The  type  of  flow 

treated  is  stationary,  unsheared,  homogeneous,  and  incom- 
★ 

pressible.  A  generalized  form  of  the  non-dispersive  approx¬ 
imation  is  suggested  also  for  shear  flow. 

This  investigation  originated  in  a  concern  for 
the  properties  of  fluctuating  pressure  in  a  turbulent  bound¬ 
ary  layer.  Accordingly,  a  further  objective  is  to  consider, 
in  the  light  of  the  work  on  velocity  spectra,  also  pressure 
spectra  in  unsheared  homogeneous  turbulence,  not  only  rela¬ 
tive  to  a  frame  in  which  the  mean  velocity  vanishes,  but 
also  a  measurement  frame  moving  at  a  velocity  -uQ,  say. 

We  thereby  introduce  a  convection  velocity  as  a  kinematic 
effect  unrelated  to  the  dynamics;  thus  uq  plays  a  role  sim¬ 
ilar  to  that  of  the  local  mean  flow  velocity  in  boundary- 
layer  turbulence  but  without  the  complication  of  a  shear 
flow . 


As  it  was  the  motivating  context,  we  review  briefly 
how  the  problem  of  pressure  fluctuations  in  a  turbulent 
boundary  layer  impels  consideration  of  space-time  correlations 


The  assumption  that  the  turbulence  is  stationary  as  well  as 
homogeneous  and  isotropic  implies,  strictly  speaking,  the 
presence  of  isotropic,  statistically  homogeneous  energy 
sources.  Alternatively,  we  may  regard  the  turbulence  as 
decaying,  but  with  a  characteristic  time  large  compared  to 
any  time  intervals  considered. 
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of  velocity.  Consider  the  wavenumber- frequency  spectrum 
of  pressure  in  the  plane  bounding  a  turbulent  flow,  P(R,a>) 

[R  =  (k^jky].  If  the  turbulent  eddies  were  of  frozen 
shape,  an  eddy  convected  downstream  at  speed  u  would  gen¬ 
erate  pressure  fluctuations  at  a  given  frequency  co  only 
via  its  component  having  streamwise  wave  number  k^  =  ru/u; 
hence  if  the  convection  speed  at  no  depth  in  the  flow  ex¬ 
ceeds  the  as,'~'Dtotic  flow  speed  U  ,  the  spectrum  P(R,cd) 
would  contain  only  wavenumbers  K  2.  ud/U^  (K  -  |  K  |],  In  actu¬ 
ality,  though  the  wavenumber  spectrum  peaks  above  uj/U^,  it 
does  not  vanish  at  lower  wave  numbers  but  has  a  "tail"  there 
on  account  of  non-convective  effects,  i.e.,  distortion  and 
decay  of  eddies. 

In  certain  applications,  however,  some  wavenumber 

range  in  this  tail  with  K<<uuATw  is  heavily  weighted  relative 

to  the  range  K^/U*.  This  circumstance  is  true  of  the  frequency 

spectrum  of  pressure  averaged  over  a  large  area,  e.g.,  over  the 

face  of  a  flush-mounted  transducer  of  radius  R  where  ojR  /U  »tt, 

o  o  * 

by  virtue  of  spatial  averaging  of  the  short-wavelength  compo- 
—  1  irk 

nents  with  K^ttR^  .  It  is  true  also  of  the  frequency  spectrum 
of  pressure  on  an  area  or  at  a  point  shielded  from  the  flow  by 
a  layer  of  material  of  thickness  L>>Uob/uu  and  of  large  lateral 
extent,  at  low  Mach  number,  on  account  of  acoustic  attenuation 
of  disturbances  having  wave  lengch  smaller  than  that  of  sound 
in  the  intervening  layer.  Hence,  to  achieve  a  theoretical 

For  elaboration  of  this  discussion  see  Chase  (1965). 

^,-^1 

The  spectrum  of  average  pressure  due  to  wave  numbers  Kn^/U^ 
in  this  situation  has  area  dependence  as  R"3. 
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account  of  the  effect  of  area  averaging  for  a  flush  area 
and  of  acoustic  averaging  for  a  shielded  area,  we  must  know 
how  the  wavenumber  spectrum  of  pressure  declines  with  de¬ 
creasing  K  in  the  non- convective  tail  where  K< 

Furthermore,  this  tail  is  relevant  to  the  scaling  of 
boundary- layer  pressure  spectra  with  the  flow  parameters,  which 
has  not  been  unambiguously  established  through  the  whole 
domain  of  interest.  Specifically,  by  the  usual  relation 
between  the  pressure  and  its  fluctuating  velocity  sources 
[ Eq . (A- 1 )  ] ,  contributions  to  pressure  are  attenuated  ex¬ 
ponentially  with  source  depth  x ^  and  parallel  wave  number 
K,  i.e.  as  exp(-Kx2).  At  relatively  high  frequencies, 
therefore,  the  convective  contribution  to  pressure,  having 
will  derive  mainly  from  velocity  sources  not  far 
outside  the  viscous  sublayer.  The  characteristic  length 
involved  by  this  contribution  is  thus  likely  to  be  the  sub¬ 
layer-thickness  parameter  v/v*,  where  v*  is  the  friction 
*  ★ 

velocity.  On  the  other  hand,  the  non-convective ,  low- 
wavenumber  contribution  will  derive  from  sources  at  greater 
average  depth,  perhaps  extending  through  a  substantial 
fraction  of  the  boundary  layer.  The  characteristic  length 

if 

In  the  very  low-wavenumber  range  K^co/c ,  wh  re  c  is  the 
sound  velocity  in  the  fluid,  it  becomes  necessary  to 
consider  compressibility. 

This  conclusion  would  be  certain  if  all  the  coordinate 
spatial  scales  of  fluctuating  velocity  components  in  the 
transition  layer  varied  as  wall  distance  X2 ;  if  some 
scales  are  determined  by  the  large  eddies ,  however ,  they 
will  be  of  the  order  of  the  displacement  thickness  in¬ 
stead  of  the  wall  distance.  Recently,  Bradshaw  (1965) 
has  contributed  to  a  grasp  of  this  scaling  dichotomy. 
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involved  by  this  contribution  may  thus  be  rather  the  displace¬ 
ment  thickness  5*.  Therefore,  the  scaling  of  the  resultant 
average-pressure  spectrum  depends  at  each  frequency  on  the 
appropriately  weighted  magnitude  of  the  wavenumber  spectrum 
of  pressure  P(K,ol'),  and  hence  of  fluctuating  velocity,  in  the 

non -convective  range  K<ayTJ  relative  to  that  in  the  convective 

★  00 

range 


The  scaling  of  the  two  contributions  could  be  nearly  equivalent 
if  both  are  independent  of  the  respective  length  parameters; 
the  frequency  spectrum  of  point  pressure  then  would  vary  as  id-1, 
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ROLES  OF  CONVECTION  AND  DISPERSION  AND  THE  LOCAL- 
CONVECTION  APPROXIMATION  FOR  EULERIAN  SPACE-TIME 
CORRELATIONS  AND  SPECTRA  OF  VELOCITY _ 


Let  the  Eulerian  space-time  correlation  function 
between  a  quantity  a(x,t)  measured  at  position  x  at  time  t 
and  a  quantity  0fx+r,t+T)  measured  at  x+r,  t+T  in  a  given  co¬ 
ordinate  frame  be  represented,  on  assumption  of  a  stationary 
homogeneous  process,  by 

(2-i)  <a(x,t)p(x+r,t+t)>  =  Wa^(r,T), 

~  _  i 

where  <>  denotes  an  ensemble  average.  Let  Wap(r  ,t)  represent 
the  correlation  function  for  the  same  two  measurement  events 
where  ?'  refers  to  the  spatial  separation  in  a  frame  moving 
relative  to  the  first  at  the  constant  velocity  -u.  The 
Galilean  transformation  between  fiames  fields  the  familiar 
functional  relation 

(2-2)  (F, T )  =  Wap(r-uT,T). 

This  rudimentary  relation  reflects  the  kinematic  character  of 
convection.  In  the  wavenumber-frequency  domain,  the  corres¬ 
ponding  relation  between  the  respective  four-dimensional  trans¬ 
forms  becomes 

(2*3)  E4a0^,x)  =  E4a0(k»^'k)- 

We  turn  explicitly  to  stationary  homogeneous  turbu¬ 
lence  and  consider  the  Eulerian  correlation  function  (2-1) 
in  the  mean  rest  frame  with  a  and  0  taken  to  be  components  of 


fluid  velocity,  say  and  v^  .  Since  we  shall  consider  the 
universal  equilibrium  range,  we  prefer  ordinarily  to  deal  with 
the  velocity  structure  function  (multiplied  by  If);  we  denote 
this  two-point,  two-component  tensor  by  t^j(r,T): 

(2-4)  ^(r.-r)  =  %  <[vt(x  +  r,  t+r)  -  vi(x,t)  ]  [v^  (x+r ,t+T) 

-Vj(x,t)]>. 

For  i  e  j,  this  quantity  is  the  decrease  of  the  autocorrelation 
from  unity;  we  shall  call  such  quantities  "decorrelations." 

:  *  * 

2.1  SEPARATION  OF  THE  LARGE-SCALE  CONVECTIVE  EFFECT 

We  now  restrict  consideration  to  spatial  separations 
r  and  time  delays  t  that  correspond  to  the  universal  equilibrium  * 

range,  namely 

(2-5)  r  «  L,  vQ  |  t  |  «  L, 

where  vq  represents  the  rms  magnitude  of  fluctuating  velocity 
and  L  the  characteristic  scale  of  the  large,  energy-containing 
eddies.  Together,  the  conditions  (2-5)  insure  that  the  spatial 
separation  of  the  two  correlated  space-time  points  in  a  frame 
convected  with  the  local  velocity  (~vQ)  associated  with  the 
large  eddies  (an  average  of  the  local  velocity  at  the  two  points) 
is  small  compared  to  the  size  of  the  large  eddies.  Later  we 
shall  restrict  consideration  to  Reynolds  numbers  high  enough 
for  existence  of  an  inertial  subrange  where  viscous  effects 
are  negligible. 
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It  is  an  underlying  assumption  of  the  Kolmogorov 
(1941)  theory  of  the  universal  range  that  the  large  eddies 
convect  the  small  eddies  without  directly  distorting  them. 

A  related,  and  probably  consequent,  assumption  is  that,  viewed 
in  the  local  frame  in  which  the  motion  due  to  the  large  eddies 
vanishes,  the  small  eddies  are  statistically  independent  of 
the  large.  These  points  have  been  emphasized  by  Kraichnan 
(1964)  and  were  recognized,  for  example,  by  Heisenberg  (1948), 
Von  Weizsacker  (1948),  and  Silverman  (1957). 

We  may  thus  consider  a  subset  of  realizations  of 
the  velocity  field  all  corresponding  to  approximately  the  same 
local  velocity  v  in  a  given  space-time  region  of  small  size 
(defined  by  (2-5))  and  take  a  statistical  average  of  the  two- 
point  velocity  product  over  this  partial  ensemble.  By  the 
stated  assumptions,  the  resulting  velocity  correlation  tensor, 
if  expressed  as  a  function  of  coordinates  in  the  frame  with 
velocity  v,  will  be  independent  of  v.  We  shall  call  the 
Eulerian  correlation  function  referring  to  this  frame,  which 
is  locally  co-moving  for  any  given  convection  velocity  v  due 
to  the  large  eddies,  the  intrinsic  correlation  function,  and 
designate  the  corresponding  decorrelation  tensor  as 
[  cf .  (2-4)].  For  r=0,  i^j(0,:)  is  the  time  correlation  in  a 

frame  which  for  each  realization  may  be  taken  as  that  in  which 
the  velocity  at  the  initial  time  t  vanishes;  the  Eulerian 
time  correlation  so  defined  is  closely  related  to  the  Lagrar.gian 
time  correlation  and  is  equal  to  it  in  the  limit  of  small  * 


and  perhaps  for  all  t.  We  may  alternatively  view  the  partial 
ensemble  average  at  fixed  v  as  a  spatial  average  over  a  re¬ 
gion  of  dimensions  larger  than  r  (the  separation  in  the  co¬ 
moving  frame)  and  larger  than  vqt,  but  small  compared  to  L. 

In  view  of  relation  (2-2)^  for  fixed  v  the  decor¬ 
relation  tensor  in  the  mean  rest  frame  is  related  to  the 
intrinsic  decorrelation  by 

^ij(r,T)  -f^Cr-vT.T). 

_  T_ 

We  let  P(v)d  v  denote  the  probability  that  the  fluid  velocity 

- 1 

in  the  mean  rest  frame  lies  in  the  three  dimensional  element  ' 

3—  ' 

d~V  of  velocity  space.  Performing  now  the  average  over  the  1 

local  convection  velocity  v  due  to  the  large  eddies,  which 
is  identified  in  the  approximation  in  question  with  the  total 
fluid  velocity,  we  obtain  the  basic  relation  of  the  standard 
Eulerian  space-time  correlation  (2-4)  to  the  intrinsic  decor¬ 
relation  : 

(2-6)  ^(r.i)  =  fd^vTty)  (r-vx ,  t)  . 

Setting  t  =  0,  we  note  the  obvious  relation 

(2-6a)  *1:j(r,0)«  ^(7,0). 

The  pure  spatial  decorrelation  ¥jj(r,0),  we  recall, 
according  to  Kolmogorov's  principles  is  an  isotropic  tensor. 

Likewise,  by  the  related  earlier  assumption,  the  intrinsic 
decorrelation  ¥^j(r,  t)  is  isotropic  even  for  t^O.  It  does 


12 


* 


not  generally  follow,  however,  that  ^j(r,T)  is  isotropic 
for  i=0 ,  in  spite  of  condition  (2-5),  since  in  (2-6)  the 
probability  density  P(v)  referring  to  the  large  eddies  may 
not  be  isotropic,  i.e.,  it  may  depend  on  the  direction  of  v.  * 

Relation  (2-6)  separates  out  the  purely  kinematic, 
convective  effect  of  the  large-scale  eddies  from  the  instrinsic 

kkk 

small-scale  properties  of  the  turbulence.  The  equivalent 
of  relation  (2-6)  was  given  in  the  present  context  by 
Silverman  (1957).  An  instance  of  use  of  a  relation  of  the 
same  sort  where  the  function  having  the  role  of  ^j(r,i)  does 
not  depend  on  the  argument  t  is  provided,  as  Silverman  noted, 
by  the  standard  treatment  of  the  scattering  of  radio  waves  by 
aggregates  of  discrete  scatterers  in  random  motion  (e.g., 

Kerr  (1951)). 

While  the  separation  (2-6)  in  the  universal  range 
is  doubtlessly  valid  to  some  approximation,  it  is  not  entirely 
clear  how  far  this  approximation  may  be  pushed  without  resort 
to  Kolmogorov's  principles  in  a  questionably  strong  form. 


We  could  also  abandon  the  earlier  assumption  of  homogeneity 
in  the  domain  (2-5)  and  write  a  generalization  of_Eq.  (2-6) 
where  P(v)  would  depend  also  on  position  and  'i..(r,T)  would  no 
longer  be  homogeneous  unless  =0.  J 


** 


Since  f^j(r,T)  is  symmetric  in  i  and  j,  according  to  (2-6) 


so  also  is  *jj(r>  )  even  if  P(v)  is  anisotropic, 


★★★ 


This  distinction  between  fluctuations  in  local  convec¬ 
tion  velocity  and  distortion  of  eddies  has  been  empha¬ 
sized  by  Fisher  and  Davies  (1963)  in  connection  with 
intense  shear  turbulence. 


1 
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In  particular,  it  is  not  obvious  that  such  a  relation,  with 
¥ij(r,T)  taken  independent  of  v  and  of  all  properties  of  the 
large  energy-containing  eddies,  is  valid  beyond  the  local- 
convection  approximation  (to  be  considered  in  the  following 
section)  in  which  is  replaced  by  ¥^(r,0).  This 

point  is  discussed  further  in  the  Appendix  [see  also  Kraichnan 
(1964)].  Nevertheless,  the  presumption  that  the  intrinsic 
decorrelation  Hf^(F,  x),  for  a  suitably  defined  local  velocity 
v,  has  the  independence  in  question  has  long  received  tenta¬ 
tive  acceptance  [Von  Weizsacker  (1948),  Heisenberg  (1948), 
Kraichnan  (1959)]  and  is  the  simplest  plausible  explicit 
assumption;  therefore,  we  accept  it  here. 

We  record  for  further  reference  the  relations  that 
follow  from  (2-6)  between  spectra  in  the  proper  and  mean 
rest  frames.  First,  we  introduce  the  trace  of  the  intrinsic 
decorrelation  tensor : 

(2-7)  nr,T)=ii.(Fsr), 

where  a  sum  is  understood  on  repeated  indices  and  ¥,  being 
an  isotropic  scalar,  depends  on  r  only  via  r(^  r  |) .  ?  is 

thus  related  to  the  trace  W(r,x)  of  the  proper  correlation 
tensor  Wj^(r,x)  by 

(2-8)  'F(r.x)  «=  vQ2  -  W(r,x). 

A  proper  wavenumber  spectrum  of  turbulent  energy  generalized 
to  time  delays  i  /0  is  defined  in  the  usual  way  (e.g.Hinze 
(1959))  by 
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(2-9) 


E(k,T)  =  (2ir)"2k2  /  d3?  e‘it,r  W(r,T> 

.1  00 

■  t r  /  dr  kr  sin  kr  W(r,T). 

o 

A  similarly  defined  three-dimensional  spectral  correlation 
Etj(Ic,T),  on  account  of  isotropy  in  the  proper  frame,  is 
related  to  E(k,T)  of  (2-9)  according  to 

(2-10)  =  (2*)'3  /d3?  W^CF.t) 

(2-11)  -  (4t)_1(6  -k1k,k'2)k'2  E (k, t) . 

We  may  define  the  complete  wavenumber -frequency  spectrum  of 
energy  for  the  proper  frame  by 

(2-12)  E4(k,o))  =  (Zir)-1  f  dTeiu)T  E(k,T), 

—  00 

and  likewise  the  spectral  tensor  E^^^  (1c,(jd)  which  is  identically 
related  to  E^(Ic,t).  Eq.  (2-11)  yields 

(2-13)  E4ij(E,a))  =  (47r)_1(51j-ktkjk-2)k-2E4(k,a>). 

From  the  inverses  of  (2-9)  and  (2-12),  the  decorrelation  trace 
is  related  to  the  wavenumber- frequency  spectrum  of  energy  of 

(2-14)  ?(r,T)  =  (2t r)'1/  dco/d3^  [l-ei(^*r”u>T)  ]  k’2  E4(k,oi). 

—  00 

It  is  sometimes  convenient  to  express  this  as 

^  00  00 

(2-15)  i(r,i)  =  4  /  d <jj  /  dk  [  (l-cosioT)sin  kr/kr+(l-sin  kr/kr) ) 

o  o 

x E4 (k,a) . 
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For  the  rest  frame,  analogously  to  (2-10),  the 
three-dimensional  spectral  correlation  is  defined  by* 

(2-16)  E^flt.T)  =  (2ir)'3  /  d3? 

Substitution  into  (2-16)  of  Eq.  (2-6)  for  Wj.  in  terms  of 

#v  ~ 

W ,  transformation  of  the  integration  vector  from  r  to 

_  A* 

r-VT,  and  identification  of  a  factor  E^(7,t)  in  the  result 
yield 

(2-17)  EtJ (k,x)  -  (4ir)"1(61J-k1kJk’2)k"2E(C,T), 

where 

(2-18)  E(E,t)  =  E  (k ,  t  )M  (Ict  ) 
with 

(2-19)  M(Et)  =  /d3^  P(7)  e_i**VT. 

As  in  the  proper  frame,  we  may  define  a  complete  spectral 
tensor  by 

(2-20)  EAij(E,co)  -  (2TT)'1  /  dTe^E^ff.T). 

Inserting  (2-17)  for  E^j  (Ic, we  obtain  the  rest-frame 
analog  of  Eq.  (2-13)’. 

(2-21)  E4ij(E,u))  =  (47r)-1(5ij-k1kjk-2)k-2E4(IT,a)) 

where 

(2-22)  ^(E.uj)  =  (2»)'1/”  dTela,TE(t,T). 


Joint  wavenumber  frequency  spectra  of  the  kind  defined  here 
were  introduced,  for  example,  by  Bass  (1954). 


From  (2-18)  we  then  find,  in  accord  with  (2-3), 


\ 


I 

\ 


t 


(1-23)  E^(lc,co)  =  /d^v  P(v)  E^(k,co-v«t) . 

Eq.  (2-23)  (or  (2-18)]  is  the  spectral  equivalent  of  the 
basic  relation  (2-6).  By  the  inverse  of  (2-16)  and  (2-20), 

(2-24)  ^lj(7,T)=(4Tr)-1/°°cico/d3ir(5ij-k1kjk-2)[l-ei(ir'r-U)T)] 

*  k"2E^(l!,u>)  , 

whence  also  the  trace  f(r,T)  is  related  to  the  energy  spec¬ 
trum  E^(lc,u))  as  in  the  proper  frame  analog  (2-14).  The 
usual  wavenumber  spectrum  of  energy  at  fixed  time,  E(k),  we 
note,  may  be  written  as 

^  00  ^ 

(2-25)  E(k)=E (k,0)  =  /  da)  E4(k,ou) 

-00 

=  E(lc,0)  =  /  du£4(E,u>). 


If  the  distribution  P(v)  of  the  large-scale  turbulent  flow 
is  anisotropic,  E4^  (E, uu),  unlike  E4^.(E,u>),  is  not  an  iso¬ 
tropic  tensor,  but  according  to  (2-21)  it  is  the  product  of  an  iso¬ 


tropic  tensor  and  an  anisotropic  scalar,  E4(lc,u));  by  (2-17)  the 
same  is  true  of  E^^,-:).  It  is  useful  to  define  also  the  fre¬ 
quency  transform  of  the  spatial  correlation  tensor, 


(2-26)  0..(r,O  :  (2-)’1  j  d’eiUJ  Wij(r>  > 

*■  J  -  or 

.  (4-)"1  J  d"*E  (‘ij_kikjk  ^ 


Since  ¥(r,x)  is  isotropic,  we  may  define  longitud- 

inal  and  transverse  decorrelation  functions  ^(r.x)  and 

*■>* 

^g(r ,  x)  ,  referring  to  for  i=j  with  the  direction  i 


respectively  parallel  and  orthogonal  to  r 
(2-27)  ¥(r,x)  .  ¥f(r,x)  +  2  ^(r.x). 

The  continuity  equation  implies 


By  definition, 


(2-28)  ¥g(r,T)  ■  ¥f  (r , t)  +  Ifr  di ^ (r , t )/dr . 

In  terms  of  and  ¥  ,  we  have 

(2-29)  itj(r ,x)  =  rirjr"2(ff(r,i)-fg(r,T)  ]  +  ^gCr.x), 
and  Eq.  (2-6)  becomes 


(2-30)  ^ij (r,T)  "  /d3vP(v)  (PiP^p"2 lif (p, T)-Yg(p , t) ]+5± ig(p, t) } , 
where  ”  =  r-vx . 

Now  we  assume  approximate  isotropy  and  define  P(v)v2dv 
as  the  probability  for  a  velocity  magnitude  between  v  and  v+dv, 
whence 

(2-31)  P(v)  =  (4tt)"LP(v). 


Acceptance  of  (2-31)  in  (2-19)  yields 


(2-32)  M(Ex)^  M(kvQx)  s  /  dv  v2P(v)  sin  kvx/kvt, 

o 

where  the  argument  kvQx  recognizes  that  the  scale  of  v  is  the 
rms  value  v  .  We  may  rewrite  relation  (2-18)  as 

*Thus  the  normalized  longitudinal  correlation  function  f(r,x) 
referi Ing  to  the  proper  frame  may  be  written  f(r,x)  -  1-(3/vq  ):^(r, 
the  function  commonly  denoted  by  £(r)  is  T(r,0). 


18. 


(2-33)  E(k,T)  .  E(k,T)M(kvQT), 

as  given  by  Silverman  (1957)  (also  see  Heisenberg  1948  , 
Wandel  and  Kofoed-Hansen  1962  ).  We  may  similarly  write 
(2-23)  as 

(2-34)  E4(k,o))  =  /d3vP(v)  E4(k,a>-7-E) 

00  _  1 

=  \  j  dv  v  P(v)  /  duE,  (k,u)-kvn)  . 
o  -1  4 

Finally,  in  this  case  of  complete  isotropy,  the  trace  function 
Kr,  i)  =  f^Cr.T)  of  the  decorrelation  tensor  in  the  rest 
frame,  together  with  the  continuity  equation,  suffices  to  de¬ 
termine  this  tensor  (and  depends  on  r  only  via  r),  just  as 
¥(r,i)  determines  ¥^(r,  x).  Hence,  the  relation  (2-6)  be¬ 
comes  equivalent  to  the  relation 

(2-35)  'i'(r.T)  =  /d3vP(v)¥(|  r-vr  |,r). 

The  one-point  velocity  distribution  P(v)  is  ex¬ 
pected  to  be  normal,  on  the  basis,  roughly,  of  the  central 
limit  theorem;  within  experimental  error,  this  result  has 
long  been  experimentally  verified,  as  discussed  by  Batchelor 
(1956).  A  nirraal  distribution  would  be  expected  for  the 
components  v^  of  v  even  in  the  event  of  anisotropy  where  the 
principal  axes  of  the  ellipsoid  of  constant  probability  would 
be  unequal.  In  the  approximation  of  isotropy,  the  normal 
distribution  is  given  by 


(2-36) 


P(v)  =  33' 2  ( 2/r) ^  vo"3exp(-3v2/ 2vq2); 
then  the  convection  factor  M  of  (2-32)  and  (2-33)  becomes 


(2-37)  M(kvo'0:  exp(-k2v2T2/6). 


2.2  LOCAL-CONVECTION  APPROXIMATION  FOR  SPACE-TIME 

CORRELATIONS  BASED  ON  SIMILARITY  AND  THE  CHARACTER 
OF  VELOCITY  DISPERSION 


Having  delineated  the  relations  between  correla¬ 
tions  and  spectra  in  the  proper  and  rest  frames,  we  consider 
the  intrinsic  decorrelation  Y  (  r ,  t  )  on  the  basis  of  the  usual 
dimensional  arguments  for  the  universal  range.  We  also  con¬ 
sider  the  deviation  of  i(p,T)  from  the  zero-delay  value 
'{'(PjO)  for  the  range  of  p  =  |  r-vr  )  that  is  significant  in 
the  integral  (2-6)  with  r  and  t  restricted  by  (2-5).  We 
thereby  tentatively  justify  the  local-convection  approxima¬ 
tion  for  relating  the  rest-frame  space-time  correlation  to 
the  purely  spatial  correlation. 

Because  of  convection  [Eq.  (2-6)],  the  rest-frame 

decorrelation  ^(FjT)  depends  on  the  large  eddies  via  the  rms 

* 

fluctuating  velocity  vq  even  in  the  universal  range.  By  the 

Unless  P(v)  is  isotropic,  YCr,')  depends  on  still  other 
characteristics  of  the  large  eddies. 
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discussion  at  the  beginning  of  Sec.  2.1,  however,  the  intrinsic 

•• 

decorrelation  is  independent  of  the  large  eddies  (Von  Weizacker 
1948,  Heisenberg  1948]  see  also  Kraichnan  1964).  Hence,  by 
Kolmogorov's  similarity  principles,  li(r,T)  in  the  universal 
range  (2-5)  must  have  the  form  (Chandrasekhar  1956) 

(2-38)  ?(r,T)  .  (£v)1/2  F(51,?2> . 


where  F  is  a  universal  function  of  the  arguments 
?  .  £l/4  v-3/4r>  .  g1/2v-1/2t 


with  e  the  energy  dissipation  rate  per  unit  mass  and  v  the 
kinematic  viscosity. 

We  now  restrict  consideration  to  the  inertial  sub¬ 
range  where,  in  addition  to  conditions  (2-5),  either  the  separa¬ 
tion  r  or  the  rms  wander  in  the  proper  frame  during  t  is  large 
compared  to  the  microscale.  This  condition  is  expressed  by 


It  may  suffice  here,  since  i(r,x)  refers  to  the  co-moving  frame, 
to  define  the  universal  range  by  r«L  and  t  in 

place  of  (2-5);  this  condition  on  t  implies  a  pseudo-Lagrangian 
excursion  small  compared  with  L  and  is  slightly  weaker  than 
v  |  t.  k<L,  since  it  may  be  expressed  roughly  as  (vj  t 
In  any  case,  in  one  degree  of  another  it  is  required  that 
vj  t  k<L,  since  if  v  |  t  (-L,  an  initially  co-moving  frame  after 
time  t  no  longer  removes  the  influence  of  even  the  large, 
energy-containing  eddies. 

**  1/3  -2/3 

I.e. ,  if  e  '  x  '  M  ^1,  the  rms  Lagrangian  separation  of 
the  fluid  elements  at  the  correlation  points  (x,t),  (x+r,  t+r) 
is  ~e  |t  if  el/3r-2/3  |T  j  ,  this  separation  is  rather  ~r 

(see  discussion  below).  Condition  (2-39)  thus  insures  that 
the  effective  separation  is  large  relative  to  the  microscale. 
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(2-39) 

where 

(2-40) 


(r2  +  g  |  r  |3)^  »  £q  or  (f^2  +  |  ^  |3)^  »  1, 
£o  «  *3/V* 


Later  we  shall  consider  separately  the  viscous  range.  In 
the  inertial  subrange,  as  well  known,  independence  of  v 
implies  the  more  restrictive  form 

(2-42)  i(r,T)  s  (er)2/3F(z), 

where  F  is  a  universal  function  of 


(2-43) 


,l/3r-2/3 


Introducing  vr  =  (rr)3/3,  which  to  within  a  constant  is  the 
velocity  dispersion  ¥(r,0)  over  a  distance  r  at  fixed  time, 
we  may  write  (2-42)  also  as 

y(r,x)  =  v2  F(vrT/r). 

Referring  to  the  limiting  cases z  -  0  andz-*  ®,we  have  the  stan¬ 
dard  forms  implied  by  similarity  arguments, 


(2-44)  f(r,0)  =  AQ(cr)2/3, 

(2-45)  ¥(0,t)  =  Bg  |  t  |, 

whence 

(2-46)  F (0)  =  Aq,  lim  [z"1F(z)J  :  B, 

oo 

On  assumption  of  an  isotropic  convection-velocity  distribution 
P(v),  we  note,  relation  (2-6)  and  form  (2-42)  imply  a  rest-frame 
decorrelation  of  the  form 

(2-41)  f(r,T)  r  (er)2/3H(v0T/r,z)  =  (er)2/3  H(voi/r,vr i/r) . 


where  Aq  and  B  are  positive  universal  constants,  presumably  of 
the  order  of  unity.  Experimental  results  quoted  by  Corrsin  (1963) 
for  the  constant  N  in  the  ordinary  wavenumber  spectrum 

(2-47)  E(k)  =  Ne2/3k  "5/3 

yield  N  2  1.5,  whence,  Aqz  3.6;  likewise  a  semi-empirical 
analysis  by  him  suggests  4.7.  These  estimates,  especially 
the  latter,  must  be  regarded  as  rough,  but  at  any  rate  there 
is  no  suggestion  of  a  gross  disparity  between  the  values  of 
Aq  and  B  in  (2-44)  and  (2-45). 

F(z)  in  (2-42)  increases  from  the  value  Aq  at 
z  z  0  to  become  Bz  at  large  z.  It  is  reasonable  to  presume  a 
smooth  behavior  for  F^z)  in  between.  On  this  presumption,  since 
B~Aq,  we  may  tentatively  suppose  that  F (z )  differs  from  Aq  by 
only  a  small  fraction  where  z  «  1.  This  supposition  can  be 
confirmed  only  by  suitable  comparisons  with  experiment,  but 
it  is  useful  to  give  a  brief  heuristic  discussion.  ¥(t,i) 
differs  from  ¥(p,0),  i.e.  F(zp),  with  z  =  2//3|  t  |, 

differs  from  Aq,  on  two  conceptual] y  distinguishable  accounts. 

First,  there  is  a  velocity  difference  between  the  points  separ- 

2  ~  2/3 

ated  by  p  whose  mean  squared  magnitude  is  v  *(^,0)  -  AQ(pp)  '  . 


Values  N  =  1.77  and  B  =  17.5  have  recently  been  estimated 
theoretically  by  Kraichnan  (1966). 
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This  velocity  acts  as  a  residual  differential  convection  vel¬ 
ocity  that  cannot  be  removed  by  the  motion  of  the  proper  frame 
and,  during  -^produces  an  rms  displacement  v^i;  hence,  the 
mean  squared  magnitude  of  the  velocity  dispersion  between  the 
two  points  occurring  on  this  account  during  t  is  ~¥(v  |  t  |,0)~ 
AQ(€Vp  I  t  l)2/^  -  AQ(eP)2//^z p^/3 .  Second,  the  velocity  even 
at  a  fixed  point  in  the  proper  frame  does  not  remain  constant, 
so  that,  even  if  p  -*■  0,  some  decorrelation  occurs  during  time 
t.  The  mean  squared  magnitude  of  this  change  is  -  T(0,t)  : 

Be  |  t  |  =  B(ep)^/\p.  Summing  the  two  contributions  in  quadrature, 
we  obtain  an  estimate  of  the  total  velocity  dispersion  due  to 
time  delay  as  a  fraction  of  that  for  zero  delay: 


(2-48)  [¥(p,t)  -  ¥  (p  ,0) )/  UP,0)  =  [F(zp)  -  F(0)]/  F(0) 

-  *p2/3  +  (B/A0)Zp. 

This  fraction,  in  fact,  is  small  if  z^  «  1.  Actually,  by  an 
argument  given  in  the  Appendix  [see  also  Sec.5],  the  contribu¬ 
tion  in  (2-48)  from  spatial  velocity  dispersion  more  likely 
2  2/3 

varies  as  z^  than  as  z  '  for  small  z^,  and  the  contribution 

from  velocity  wander  may  not  enter  in  quadrature  but  in  such 

a  way  as  *>  be  of  still  higher  order  (xzp),  though  at  large 

z  it  must  vary  as  z„  and  predominate.  In  fact,  in  Sec.  3.1 

on  the  basis  of  a  specific  conjectured  intrinsic  spectrum 

E^(k,oj)  and  without  reference  to  the  heuristic  distinction  of 

.spatial  velocity  dispersion  and  velocity  wander,  it  is  found 

that  (F(zd)  -  F (0)  ]/F (0)  tr  0.47(B/A  )2z  2  ,  and  apart  from 
h  op 
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the  value  of  the  numerical  coefficient  the  same  result  holds  for 
a  broad  class  of  possible  E^(k,u>). 

For  the  moment  we  suppose  P(v)  is  isotropic;  it 
then  suffices  to  consider  the  trace  function  i(r,x).  Inser¬ 
ting  (2-42)  in  (2-35)  we  have* 

(2-51)  ¥  (r ,  t)  =  e2/3f  d3^  P(v)  |  r  -  v  x  |2/3F(e1/3  |  x  |  |  r-vx  f2/3). 

In  view  of  the  behavior  F(z)-*  Bz  as  z  -*■  <*,,  there  is  no  singular 
behavior  of  the  integrand  at  r  -  vx  =>  0.  When  the  integrations 
over  angles  of  v  are  performed,  the  spatial  and  convective 
displacement  vectors  r  and  vt  roughly  speaking,  add  in  quad¬ 
rature.  With  the  final  integration  over  v,  the  characteristic 
order  of  magnitude  of  |  r-vx  |  in  the  integrand  is  simply 

(2-52)  R  =  (r2+v  2T2)* 

The  argument  z  involved  in  the  integrand  is  thus  typically  of  tne 


Form  (2-51)  contrasts  with  that  assumed  by  Chandrasekhar 
(1956)  in  his  theory  of  turbulence.  By  his  contention 


(2-49) 


i  (r,  ) 


z  =  ,l/3r-2/3 


where  a  prime  denotes  cydr.  But  integration  over  r  with  neglect 
of  the  contribution  from  the  viscous  range  then  yields 

(2-50)  i(r,~)  =  (tr)2/3  2(z)  +  i(0»")  >  where 

Z(z)  z  (3/2)z  /  dx  x"2o(x) 
z 

[cf  12-42)].  According  to  (2-50),  the  large  eddies  contribute 
to  the  decorrelation  only  additively  and  independently  of  r. 

If  the  large-scale  motion  is  isotropic  and  V ( 0,r)  depends  on 
it  only  via  vq,  for  example,  (2-50)  has  the  form 

i(r,  )  (  r)2/3  Z(z)  +  (evo  |  i  |)2/3  T(  *  i  x 

The  failure  of  form  (2-50)  to  reflect  properly  the  effect  of 

convection  by  large  eddies  was  emphasized  by  Kraichnan  (1959), 
(1959a). 


order  of 
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By  use  of  the  usual  estimate  e~vQJ/L,  we  then  have 


V<vo'  t  |/RC)(RCA)1/3. 

The  first  factor  is  less  than  unity  by  definition  (2-52), 

and  the  second  is  small  relative  to  unity  by  the  restriction 

(2-5)  to  the  universal  range.  Hence  only  arguments  z  «  1 

enter  significantly  in  (2-51),  and  for  such  arguments,  by 

our  previous  discussion,  the  fractional  error  will  be  small 

if  we  replace  F  in  (2-51)  by  F(0)  [ EAq ] .  In  accord  with  * 

(2-6a)  we  identify  (er)^^  F(0)  as  ¥(r,0).  We  must  insert 

Rc  in  place  of  r  as  the  pertinent  argument  in  condition 

(2-39)  for  use  of  the  inertial- subrange  form  (2-42)  in  the 

■a 

integral  (2-51);  since  e  |  i  P  «  vq^t^  by  (2-5),  the  ineitial 
subrange  in  the  present  context  is  then  defined  by 


(2-53)  £Q  «  Rc  «  L. 

We  thus  obtain,  as  a  valid  approximation  in  the  inertial 
subrange  (2-53),  an  explicit  expression  for  the  space-time 
decorrelation  ¥(r,T)  in  terms  of  the  purely  spatial  decor¬ 
relation  T (r,0) : 

(2-54)  *(r,T)  *r  / d3^  P(^)  T(  |  r  -  vt  |,0) . 


In  other  words,  in  the  computation  of  ¥^j(r,~)  by  Eq.  (2-6), 
the  trace  T(r,x),  which  specifies  T^j(r,x)  ir.  that  equation, 
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is  approximated  by 

(2-55)  ¥  (r ,  t  )  -  f  (r  ,0)  =  ¥(r,0). 

We  refer  to  (2-54)  as  the  local-convection,  or 
non-dispersive,  approximation.  This  relation  has  been  com¬ 
pared  with  experiment  by  Favre  (1965),  who  attributes  its 
proposal  to  Kovasznay.  We  shall  consider  these  experiments 
later.  In  Sec.  3.4  we  shall  derive  explicitly  in  certain 
limits  the  dispersive  correction  to  (2-54)  given  by  (2-51) 
for  the  F(z)  derived  from  the  specific  conjectured  form  of 
E40<  ,o>)  mentioned  above.  The  dispersive  correction  to  the 
local-convection  approximation  (2-54)  for  y(r,x)  (or  to  the 
corresponding  approximations  for  i^(r,0)  is  found,  for  the 
class  of  E^(k,(u)  considered,  to  be  of  relative  order  ~(B/Aq) 
x(e  "  (B/Ao)  2(vq  ’/L) 2//3  and,  still  more  generally, 

vanishes  in  the  defining  limit  Rc/L  —  0  of  the  universal 
range.  We  discuss  in  the  next  section  limitations  on  the 
validity  of  the  wavenumber-frequency  spectrum  corresponding 
to  (2-54). 

In  the  more  general  case  of  anisotropic  P(v),  the 
trace  i(r,  i)  does  not  suffice  to  determine  i^(r,i).  Exactly 
as  at  (2-42),  however,  for  the  inertial  subrange  we  may  write 

(2-56)  ^(r.t)  =  Ai(^:r)2/3  ci(z),  (i=  0,1,2) 

where  i Qz  t,  z  and  the  equation  for  i=0  repeats 

(2-42)  with  'q(  z)  =  F(z)/Ao-  We  may  make  gq(0)  :  o^(0)  -  ^(OJrl 
by  taking  (3/11)Aq  ■  (3/4)A2  -  A^  i  A  and  defining  A  accordingly. 
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The  ^^j(r,T)  can  be  formed  directly  by  inserting  y  as 
given  by  (2-56)  into  (2-30).  The  earlier  discussion  of 
Eq.  (2-51)  applies  substantially  unchanged  to  these  equations 
for  '{'^(rjT)  since  o^(z),  02(2)  are  similar  in  character  to 
F(z)  .  Hence  the  local-convection  approximation  (2-54)  may 
be  generalized  to  the  form  required  for  anisotropic  P(v): 

(2-59)  ^(r.i)  -  P(v)'Fij(F-vT,0). 

This  approximation  is  still  equivalent  to  application  of  the 
approximation  (2-55)  in  Eq.  (2-6). 

We  turn  to  consider  whether  (2-54)  land  (2-59)] 
should  be  valid  also  in  the  viscous  range.  In  more  general 
terms,  our  previous  arguments  may  be  summarized  as  follows. 

The  effective  displacement  argument  |  r-v’  j  in  Eq.(2-6)  for 
¥(r,T)  is  characteristically  ~R  .  Hence  approximation  (2-54) 
will  be  valid  if 

(2-60)  if (Rc,  )-*(Rc,0)]/  i (Rc ,0)  «  1. 

Considering  two  types  of  contribution,  we  estimate 

- 

More  explicitly,  the  continuity  equation  (2-28)  with  (2-27)  and 
(2-56)  yields  o2(z)  =  0^(2)-  1/4  zo^(z),  oQ(z)  =  o^z) -(2/11) zo^(z) . 
For  the  coefficients  =  lim  +  a>[o1(z)/z],  these  imply 
(9/ll)b^  *  <12/ll)b2  =  bQ.  From  these  results  we  find 

(2-57)  a1(z)  =  (ll/2)z11/2  /"  dz’  z' '13/2  oQ(z'). 

L 

(2-58)  o,(z)  -  (11/8) [a  (z) -(3/2) z11/2  ™  dz'  z' '13/2o(z' ) ] . 

L  0  z 
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Hence  the  condition  for  (2-54)  becomes 

(2-61)  (U!i(Rc,0)]^  I  t  1,0)  +  i(0,T))/  1(RC,0)  «  1. 

In  the  viscous  range  we  have  the  standard  result 
(Hinze  (1959)) 

(2-62)  l(p,0)  ~  5(t/30i  )p2  for  p  <  £q, 

corresponding  to  F(^1,0)  =  5(1/30) -^2  in  Eq.  (2-38)  for  r,^  ^  1. 

-  o 

Similarly,  for  ^  <!■  we  must  have  ^(O.t)  1  t  ,  i.e. 

2 

F(0,"2)  -  ^c2^2  ’  w^ere  c2  a  constanti  thus 

(2-63)  1(0, t)  -  5c 2(cv)*(t/v)  2  for  ^  (v/c)*. 

This  form  in  the  similar  case  of  the  Lagrangian  correlation 
is  implied  also  by  the  work  of  Uberoi  and  Corrsin  (1953), 
summarized  by  Hinze  (1959),  p.  321.  Identifying  the  two  cases, 
we  would  infer  from  this  work  a  value  for  C2  given  by 

(2-64)  c2  (15)"3/2(12. 75a"1  +  4.46a), 

where  a  denotes  Heisenberg’  s  dimensionless  spectral-transfer 
constant,  presumed  of  the  order  of  unity.  By  insertion  of 
(2-62)- (2-63),  condition  (2-61)  may  be  written 

(2-65)  (1/6) :  2  *  ^Oy2/^  «  1. 

where  ;  -  ^  In  the  v: '*'•  ,s  range  defined  by 


(2-66)  (r2  +  c1/2v1/2t2)1/2  lQ  or  (?2  .  §2)1/2^1, 

2 

we  have  (l/6)?2«^-  (2-66).  Condition  (2-65)  then  reduces  to 


(2-67)  30C2  |(ev)Vv02)ll  +  (r/v^)2)*1  «  1. 

But 

(2-68)  vo2/(ev)^  =  (3//l5)Rex, 

where  Re^  is  the  usual  turbulence  Reynolds  number  based  on 
transverse  dissipation  scale,  and  Re^  >>  1  for  the  assumed 
existence  of  the  universal  equilibrium  range.  According  to 
(2-64),  we  have  C2  ~  Hence  (2-67)  is  satisfied,  and  the 
local-convection  approximation  (2-54)  or  (2-59)  is  expected 
to  hold  in  the  viscous  range  (2-66),  as  well  as  in  the  iner¬ 
tial  subrange  (2-53). 

Conceivably,  (2-54)  is  less  valid  in  the  connecting 
r  gion  between  ^nertial  and  viscous  ranges.  Apart  from  this 
possibility,  the  local-convection  approximation  is  expected 
to  be  valid  throughout  the  universal  range  where  R£  «  L. 


* 
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2.3 


SPACE-TIME  CORRELATIONS  BY  THE  LOCAL-CONVECTION 
APPROXIMATION  FOR  AN  ISOTROPIC  NORMAL  VELOCITY 
DISTRIBUTION!  SPACE-TIME  ISOTROPY;  COMPARISON 
WITH  EXPERIMENT _ ' _ 

We  assume  now  a  flow  having  an  isotropic  and  normal 

distribution  of  total  fluctuating  velocity  P (v)  given  by 

(2-31)  and  (2-36).  In  such  case  we  may  define  longitudinal 

and  transverse  decorrelation  functions  for  the  mean  rest 

frame,  i'^(r,T)  and  ¥g(r,T).  The  relations  (2-27)  -  (2-29) 

must  then  hold  as  well  for  the  rest-frame  functions  (all 

tildes  omitted).  Eq.(2-6;  for  i«j  with  i  respectively 

parallel  and  orthogonal  to  r  then  yields  integrals  for 

¥ 'P  in  terms  of  .  Using  spherical  coordinates  in 

—  3—  — 

v- space  (d  v  s  dvdudf?)  with  r  as  polar  axis,  and  representing 


CMf.¥g) 

by  ('Pq  , 'P  ^ , 'i  2 )  f°r  convenience  (whence  'iQ  +  2^2)  > 

we  have* 

00  _ 

1 

(2-69) 

(r,  t)  -  (1/2)/  dv  v 

P(v)  /  duY.(u), 

0 

-1  1 

with 

(2-70) 

¥Q(u)  =  *f  (P>7)  +  2vg( 

P,  0, 

(2-71) 

Y1(u)  =  (r-v  u)2p‘2['Ff( 

p,  >)  -  *g(p, *)]  +  *g(p, •), 

(2-72) 

*2(^)  s(1/2)<vt)2(1-u2) 

p'2|tf(»  ,T)  -  Yg(p,-  )]  +  Yg(p,l), 

P2  :  r2  +  (vt)2 

-  2rvin. 

* - 

As  a  check,  it  is  verified  directly  that  i^,  ¥  satisfy  the 
continuity  equation  (2-28)  provided  f£,?g  do. 
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We  now  utilize  the  local -convection  approximation 
and  consider  first  the  inertial  subrange.  We  then  have  the 
^(r,!)  given  by  Eq.(2-56)  with  the  o^(z)  replaced  by  unity. 
Inserting  (2-36)  in  (2-69)  and  integrating  over  u,  we  obtain 

(2-73)  ^(r.T)  =  3(6/7r)^A(£r)2//3  /  dx  x2  exp(-3x2/2)Li(xvoVr) , 

o 

where 

Lq(3)  :  (ll/16)3"1(s®/3-  s_8/3), 
s+  :  I  1+3  I,  s _  =  |  1-3  j, 

and  and  L2  are  given  by  similar  but  more  complicated 
expressions . 

We  readily  find  the  limiting  forms  of  (2-73): 

(2-74)  tjOr.O-Aj  (Er)2/3(l+clo@02)  for  0o«  1, 

(2-75)  f1(r,T)  -S1(ev0r)2/3(l+cla>B;2)  for  f>0»  1, 

where 

S  volTi/r’ 

(2-76)  A2  =  (4/3)A  ,  Aq  3  (11/3)A,  =  A, 

S1  z  S2  =  (2/3) 1/3 (22/9)^"^r(ll/6) A,  SQ  =  3S1> 
c10  =  H/27,  C20  r  H/108,  =  5/27, 

cloo  =  1/5.  c2oc  *  2/5,  cocc  =  1/3. 

The  characteristic  dependence  T^(0,t)  -»  S^(evQT)2/3  for  rrO 
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is  well  known  (Landau  aud  Lifshitz  1959,  Eq.  (32.2)). 

We  could  also  derive  and  ^  from  the  simpler 
Yq  by  the  continuity  equation.  Thus,  writing  Eq.  (2-73) 
functionally  as 

(2-77)  ^(r.T)  :  Al(tr)2/3G1(v0Vr), 

(whence  G^vO)«l),  from  (2-27)  and  (2-28)  we  obtain 

(2-78)  Go(t3o)  =  G1Oq)  -  (3/ll)0odG1(3o)/deo 

(cf.(2-58)).  Eq.  (2-78),  together  with  the  asymptotic  form 

Gi(0oMsi/Ai)3o2/3  as  po"*“>  yields 

(2-79)  G1(Pq)  =  (11/3)30U/3  /“  d(30'14/3Go(P). 

Since  the  displacements  r  and  vqt  must  add  nearly  in 
quadrature  for  appropriate  relative  scaling,  we  define 
functions  H^(0q)  by  writing 

(2-80)  Yi(r,x)  =  A.(£Ri)2/3[l-Hi(eo)], 

where 

(2-81)  Rt2:  r2  +  st2  (vqt)2,  Si2:  (S./A^3. 

2 

By  definition  (2-81)  of  the  scaling  factors  sL  and  the  results 
(2-74),  (2-75),  the  functions  )  so  defined  are  such  that 

(2-82)  Ht(0)=  Hi(oo)  =  0, 


•4 
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i.e.  ¥^(r,T)  s  ^(eR^)2/3  if  either  r  .  0  or  t  :  0.  From 
the  numerical  values  (2-76),  we  find 

(2-83)  s2  r  1.455,  s2  r  0.613,  s2  =  0.795. 

The  integrals  (2-73)  have  been  evaluated  numerically,  and 

the  resulting  functions  H^Oq)  of  (2-80)  are  shown  in  Figure 

- 1  - 1  * 

1  as  functions  of  Tan  [«  Tan  (v^r|/r)]. 

We  shall  refer  to  a  function  that  depends  on  r  and 
2  2  2 

t  only  via  r  +  s  (vqt)  for  some  constant  s  as  a  space-time 
isotropic  function.  Figure  1  shows  that  the  ¥^(r,T)  are  all 
closely  approximated  by  the  space-time  isotropic  functions 
defined  by  neglecting  the  in  (280),  i.e., 

(2-84)  ^(r.i)  »  Ai(cRi)2/3. 

The  corresponding  maximum  fractional  errors  are  1.3/!.,  3.5% 
and  2.4 %  for  the  longitudinal,  transverse,  and  trace  functions 
(i  =  1,2,0),  respectively. 

We  turn  momentarily  from  the  inertial  to  the  viscous 

subrange,  applying  once  more  the  local-convection  approximation 
* 

By  expressing  ^(r.T)  by  the  inverse  of  (2-9)  and  E(k,i) 
by  (2-33)  and  (2-37),  with  E(k,x)  replaced  by  E(k)  : 

Ne2/3k”3//3  in  the  present  approximation,  we  obtain 

y  (r,  t)  =  2N€2^3  /“  dk  k‘5//3(l-(sin  krAr)exp(-vQ2T2k2/6)  ] , 
whence,  according  to  Bateman  (1954),  p.74,  (24),  ¥(r,T)  can 
be  expressed  also  in  terms  of  Rummers  confluent  hypergeo¬ 


metric  series. 


(2*54).  In  Chit  InsCtnce  we  have,  at  at  (2-62) 

(2-85)  »1(p)  «  *„p2,  t2(p)  ,  2»„p2,  t0(p)  s  5»vp2, 

where  av  >  e/30v.  Eq.  (2-69),  with  P(v)  given  by  (2-36),  now 
yields 

(2-84)  Y^r.t)  ,  Aq1(c/30v)  R2t 
with 

(2-87)  R2t  =  r2  +.21(v0t)2> 

*01  =  5/i‘  *02  =  V6,  .  2„  s  1,  lef .  (2-83) ) 

A0I  :  A02  =  ^ »  *00  " 

Hence,  in  the  viscous  range  the  decorrelations  fi(r,x)are  ex¬ 
actly  space-time  isotropic. 

A  space-time  isotropic  correlation  function  was 
proposed  previously  by  Lilley  and  Hodgson  (1960),  App.D. 

For  comparison  with  experiment,  we  consider  also 

A  _ 

space-time  decorrelations,  say  Y^j(r,T),  measured  in  a  frame 
having  a  fixed  velocity  -u0  ■  (-uQ,0,0)  relative  to  the  mean 

A 

rest  frame.  By  the  kinematic  relation  (2-2),  ve  have  ^^j(r,T)  - 
f ij  (r  -3q t ,  t ) ,  In  the  present  approuimation  ^(t^t)  is 
isotropic  and  given  in  terms  of  f ^(r, t) »^g(r, t)  as  in  Eq.  (2-29). 
In  particular,  the  decorrelation  for  the  streamvise  component 
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* 


f 


4 


is  expressed  by 

(2-88)  ?n(r,T)  s  ^(r'.r)  +  R2(r*L2  +  R2)"1  [^(r1,  t)-¥£  (r')  t)  ] , 

where  r's  (ri"u0T»r2»r3)  and  R  =  r2  +Ty  In  the  inertial 
subrange  and  H  are  given  by  (2-80),  but  since  a  relative 
error  ~1%  is  negligible,  we  use  the  space- time  isotropic 
approximation  (2-84)  for  ¥£  and  the  consequent  H  [Eq.  (2-28)], 
obtaining 

(2-89)  iu(r,T)  -  A(£R1)2/3(1h^  R2/^2), 

^2  2  2  2  2  2 

where  R^  =  (r^-uQT)  +R  +  vq  ’  v  .  The  correlation  at  a 
point  in  the  mean-rest  frame  after  a  delay  r  according  to 
(2-89)  is  the  same  as  in  the  laboratory  frame  after  a  delay 
-  (s^vo/uo) t~ (v0/u0) t  5  this  result  was  given  by  Lilley  and 
Hodgson  (1960),  App.D,  and  noted  to  be  in  broad  agreement 
with  the  results  of  Favre  (1965)  for  both  grid  and  wall  turbulence. 

The  velocity  component  referred  to  here  is 
parallel  to  uq  by  definition,  and  in  the  measurement  of  present 
concern  r  is  also  nominally  parallel  to  u  .  We  therefore 
suppose  the  angle  0  [  =  sin”^(R/r)]  between  r  and  uQ  is  small, 
as  well  as  VQ/U0«  whence  (2-89)  may  be  written  approximately 
as 

(2-90)  Tn(r,T)  *  A,2/3[a2+  r202  +s2(vQ/uo)2(r+  A)2]1/3 

x  (1+3  rV/[A2  +  r202  +  s2(vo/uo)2(r+  A)2]], 
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where 

A  s  u  (t-t  )  and  t  =  r,/u  -  r/u  . 
ov  nr  m  -  r  o  '  o 

A 

Thus  is  the  time  delay  at  which  d  =  0  for  given 

r,  and  A  is  the  actual  time  delay  in  spatial  units  relative 
to  t  .  Since  flf  «  1  and  s,v  /u  «  1,  the  time  correlation 
at  fixed  separation  is  determined  mainly  by  Ae^^,  but  the 

A 

correlation  envelope  curve  dx  «0  for  varying  separation 

is  determined  mainly  by  the  product  of  and  the  function 

of  fl  and  vo/uQ  given  by  setting  A  :  0  in  (2-90).  The  rate 
of  decrease  of  this  correlation  envelope  is  sensitive  to 
the  misalignment  0,  being  significantly  increased  by  0  unless 

9  «  W 

We  compare  results  computed  from  (2-90)  with  mea¬ 
surements  in  grid-produced  turbulence  by  Favre,  Gaviglio,&  Dumas 
(1953)  , (1954),  summarized  by  Favre  (1965).  Favre  per¬ 
formed  a  similar  comparison  based  directly  on  the  nondisper- 
sive  relation  (2-64 )  (with  r  -*  r  -  u  t)  and  employing  the 
measured  Y-^.O),  Y  (rj  0).  The  present  comparison  entails 

O 

the  further  approximation  of  and  Y  by  the  inertial-sub- 
range  forms  but  eliminates  certain  unassessed  approximations 
made  by  Favre.  The  parameter  A  is  now  regarded  as  adjustable 
on  the  basis  not  only  of  the  measured  correlations  for  i  =  0 
but  the  whole  series  of  measured  space-time  correlations 
Yn(r,t)  and  their  envelope.  Use  of  the  inertial-subrange 
form,  however,  implies  that  the  computed  normalized  space-time 
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correlation  must  be  expected  to  depart  substantially  from  the 
measured  values  at  least  where  these  are  less  than  ~0.6.  In  the 
specific  experiments  referred  to, in  fact,  the  Reynolds  number 
Re^  z  44,  scarcely  sufficient  for  existence  of  an  inertial 
subrange  at  all,  and  small  enough  also  that  the  dispersive 
departure  from  the  local-convection  approximation  may  be 
appreciable,  as  discussed  in  Sec.  3.4  below. 

For  the  geometrical  and  flow  parameters  given  for 
the  Favre  experiments,  we  estimate  from  other  experiments 
on  decay  of  turbulence  downstream  of  a  grid,  notably  Batchelor 
and  Townsend  (1948),  that  the  turbulence  intensity  at  the 

ic 

pertinent  locations  was  v  /u  -  0.027.  Accepting  this 
value,  assuming  no  misalignment  (0rO),  and  adjusting  A  to 
give  a  visually  good  fit  to  the  form  of  the  measured  time 
correlations  at  various  fixed  separations  r,  with  more  empha¬ 
sis  on  agreement  where  the  correlation  is  relatively  high, 
we  arrive  at  A  -  0.55  and  the  series  of  curves  in  Figure  2 
having  the  solid- line  envelope.  The  breadth  of  the  peaks  for 
various  separations  is  fairly  satisfactory  relative  to  the 
experimental  results  shown,  but  the  envelope  of  the  peaks  de¬ 
creases  too  slowly.  This  fact  is  attributed  mainly  to 


The  dissipation  rate  e  is  then  computed  by  Eq.(2-72). 

At  separation  r  =  0,  for  which  the  experimental  curves 
given  in  Figs.  1  and  5  of  Favre  (1965)  fail  to  coincide, 
we  accept  the  latter. 


neglect  of  dispersion  and  to  experimental  misalignment.  If, 

o 

for  example,  the  latter  were  0  =  1.5  for  each  separation, 
with  vo/uq  :  0.027  still  and  A  =  0.55,  the  curves  would  be  those 

represented  -  ;  these  agree  rather  well  with  the  measured 

ones  and  their  envelope.  If  we  suppose  instead  0  -  0  but 
arbitrarily  adjust  the  intensity  as  well  as  A  to 
yield  a  fit,  for  A  >  0.41  and  vq/uq  -  0 . 04 2  we  obtain 
almost  identical  results.  In  the  measurements,  the  nominal 
transverse  separation  was  varied  to  maximize  the  correlation; 
hence  the  error  due  to  misalignment  was  presumably  small.  It 
appears  quite  plausible  that,  if  comparison  were  made  with 
experiments  at  sufficiently  high  Reynolds  number,  the  results  yielded 
by  (2-89)  for  space-time  velocity  correlations  on  the  basis  of  the 
local-convection  approximation  would  prove  in  satisfactory  agreement 
with  measurement,  and  provide  a  useful  determination  of  the  Kolmogorov 
constant  N.  Further  comment  will  be  made  in  Sec.  3.4,  where  a  dis¬ 
persive  correction  is  introduced. 


We  may  consider  also  the  effect  of  possible  anisotropy  of 
the  large  eddies.  Grant  and  Nisbet  (1957)  obtained  experi¬ 
mental  indication  that  the  streamwise  fluctuating  component 
in  grid  turbulence  tends  to  be  the  largest.  Making  the 
extreme  assumption  of  purely  streamwise  fluctuation  in 
place  of  (2-31,we  find,  for  r  parallel  to  uq,  that 
i  f  (0,  :)/i  f(r,0)  =  0. 44(vq  :/r)  2/^  place  of  1.13  (vQ/r)2^3. 
Hence  such  anisotropy  would  be  expected  to  reduce  rather 
than  increase  the  rate  of  decline  of  the  space-time  correla¬ 
tion  envelope.  This  anisotropy  has  no  effect  on  the  ratio 

o  /P  mentioned  below, 
c  s 
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The  measured  streamwise  space- time  correlations 

(R-0)  may  be  characterized  also  by  the  ratio  of  minor  to 

major  axis  for  the  streamwise  isocorrelation  curves,  which 

depends  somewhat  on  the  correlation  value  (Favre,  Gaviglio, 

and  Dumas  1953).  More  generally,  we  define  as  one  half  the 

ratio  of  r^-intercept  (t-0)  of  a  streamwise  isocorrelation 

curve  to  the  value  of  u  t  where  r,  «u  t  on  this  curve:  for  an 

o  i  o  ’ 

ellipse,  p  is  nearly  the  ratio  of  axes.  Transverse  isocor- 
s 

relation  curves  (^-0)  were  also  measured  in  the  experiments 
and  may  be  characterized  mainly  by  the  ratio  p  t  of  the  R- intercept 
(t-0)  to  the  uQT-intercept  (R-O).  Finally,  compensated  iso¬ 
correlation  curves  (r^“u0T)  were  measured  and  may  be  charac¬ 
terized  by  the  ratio  p£  of  the  R-intercept  (r^-0)  to  the  r^- 
intercept  (R-0). 

For  the  viscous  range,  which  applies  only  where  the 

normalized  correlation  is  very  near  unity,  Bass  (1954)  derived 

1  /2 

for  homogeneous  turbulence  the  general  relations  r*(a/b)  /2, 

3 

Pt-1/V*.  pc-V*v  where  a/b  is  a  small  undetermined  quantity. 

Eqs.  (2-86)  and  (2-88)  of  the  local-convection  approximation 

for  the  viscous  range  conform  to  these  results  and  imply  a/b  - 
2 

(5/3)(vQ/uo)  ,  whence  p^O.65  v0/u0-  For  inertial  sub¬ 

range,  on  the  other  hand,  Eq.(2-89)  yields  Psm2s^vo/uo='0.6lvo/uo , 
pt-  |  (1+s12vo2/uo2)1/2-0.65,  -c-(|)3/2s1vo/uo^1.30cs,  which  we  may 
compare  with  experiment.  The  measured  values  of  r  ,  without  dis¬ 
persion,  require  somewhat  too  large  VQ/U0»  as  already  discussed.  The 
measured  p  varies  in  the  vicinity  0.4  to  0.5,  a  value  appreciably 
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smaller  than  the  computed  0.65.  Most  strikingly,  the  measured 
p  / p_  varies  in  the  vicinity  of  4  to  6,  much  larger  than  the 
computed  1.3.  Introduction  of  dispersion,  anisotropy  of  the 
large  eddies,  or  streamwise  inhomogeneity  appears  unsuitable 
to  account  for  this  discrepancy.  An  obvious,  but  unlikely, 
way  to  account  for  it,  if  the  measured  value  is  not  grossly  in 
error,  is  to  suppose  the  spatial  velocity  correlations  at  fixed 
time  are  anisotropic  and  the  streamwise  component  correlated 
over  a  much  larger  distance  in  the  transverse  than  the 
streamwise  direction. 
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2.4  VELOCITY  SPECTRA  IN  THE  LOCAL -CONVECTION 

APPROXIMATION;  IMPORTANCE  OF  VELOCITY  DISPERSION 

We  consider  the  implications  of  the  local-convection 
approximation  (2-55)  for  the  rest-frame  velocity  spectra 
defined  in  Sec.  2.1.  The  corresponding  approximation  for 
the  partial  transform  E(k,T)  is 

(2-91)  E(k,T)  »  E(k,0)  ;  E(k) , 

where  E(k)  ;  E(k,0).  By  definition  (2-12),  then,  we  have 

(2-92)  E4(k,tu)  *  E(k)  5(co)  , 
whence  (2-23)  becomes 

(2-93)  E4(Tc,u>)  -  E(k)  fd^?(y)b((a  -  v)  , 

This  equation  states  that  in  the  approximation  of  local  convec¬ 
tion,  a  fluctuating  velocity  component  with  given  frequency  cjo 
and  wave  number  fc  can  be  generated  only  by  virtue  of  a  local 
convection  velocity  v  of  the  "frozen"  large  eddies  such  that 
the  projection  of  v  in  the  direction  of  k  generates  a  disturbance 
at  frequency  a>  by  motion  of  the  fixed  wave  pattern.  The  restric¬ 
tion  (2-53)  of  (2-54)  to  the  inertial  subrange,  by  the  rest-frame 
analog  of  (2-15),  implies  the  reciprocal  restriction 

(2-93a)  L-1«(k2+  cd2  /v2)  1/2«rL 

on  (2-93)  and  (2-95)  below. 

In  the  case  of  isotropic  P(v) ,  Eq.  (2-93)  is  equivalent 
to  the  following  relation  derived  from  Eq.  (2-34): 

(2-94)  E4(k,ui)-(l/2)k'LE(k)/"  dvvP(v) . 

4  a)/k 
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With  a  normal  distribution  [Eq.  (2-36)]  this  becomes 

(2-95)  E4(k,u>)  ^31/2(2T)'1/2(vok)'lE(k)exp(-(3/2)(a>/vC)k)2] 

Hence,  if  the  normal  distribution  represents  the  true  situation 
sufficiently  veil,  for  large  ^/vQk  the  local-convection 
approximation  yields  an  exceedingly  small  energy  density 
E^(k,oo)  on  account  of  the  minute  probability  of  a  local 
convection  velocity  >U)/k.  Because  of  the  deformation  and 
decay  of  eddies,  however,  we  must  suppose  that  E^(k,U))  in 
reality  is  much  larger;  i.e.,  however  good  an  approximation 
the  local-convection  approximation  may  be  for  ^ (r ,  t)  in  the 
universal  range,  it  is  not  adequate  to  yield  a  good 
approximation  to  such  a  demanding  detail  as  the  frequency-wave¬ 
number  spectrum  at  large  a>/vok. 

An  example  of  the  sensitivity  of  E^(Tc,cd)  in  this 
regime  is  provided  by  comparing  the  spectrum  corresponding 
exactly  to  the  local-convection  approximation  in  the  inertial 
subrange  with  thaf  corresponding  to  the  space-time  isotropic 
approximation  for  the  longitk.  iinal  decorrelation  ?(/ £ ( r ,  t ) 

[Eq.  (2-84)  ];  the  latter  sufficed  to  give  ^(r,-t)  within  1.3% 
of  the  value  given  by  the  former  for  all  values  of  vQT/r 
(Fig.l).  Corresponding  to  (2-56)  with  oQ  =  1,  we  have  the 

E(k)  of  (2-47)  with  N  =  (110/27) [ T( 1/3)  ]' lA  .  Then  (2-33), 
(2-37)  and  (2-95)  yield 


(2-96)  E(k, t)  =  Ne2/3k"5/3exp(-k2v2T2/6) , 

E4rk,0))  =  (3/?-r)1/2Ne2/3v’1k‘8/3exp[-(3/2)(o)/vok)2]. 

The  space-time  isotropic  approximation  (2-84)  for  ^(r^),  on 
the  other  hand,  by  a  standard  relation  (Hinze  1959,  Eq.  (3-70) 
extended  to  r  4  0)  yields 

(2-97)  E(k.-r)  =  2"11/6[r(17/6)  l'1Ne2/3k"5/3(s1voik)17/6 

xK17/6(slvolk)>  .. 

E4(k,<»)  .  ir-1''2(14/ll)[r(7/3)/r(ll/6)]N€2/3(s1v0)'1  '  \ 

xk'8/3[l  +  (a)/Slvok)2)‘10/3 

)  * 

where  s,  was  given  at  (2-83)  and  K  denotes  the  usual  exponen- 

i  V  £ 

tially  decreasing  modified  Bessel  function.  For  large  oo/vok 
the  approximation  (2-97)  for  E^(k,a)  thus  does  not  retain  the 
rapid  rate  of  decrease  given  by  the  exact  local-convection 
form  (2-96) . 

2.5  TAYLOR'S  HYPOTHESIS,  EFFECTIVE  CONVECTION  VELOCITIES,  AND 

CROSS -SPECTRA  IN  THE  LOCAL -CONVECTION  APPROXIMATION 


The  present  methods  afford  a  simple  quantitative 
treatment  of  departures  from  Taylor's  hypothesis  for  the  inertial 
subrange  in  homogeneous  turbulence,  a  topic  previously 

5 - - - — 

If  the  approximation  of  space-time  isotropy  (2-84)  is  applied 
to  the  trace  r,t)  instead  of  r,-r),  the  approximation  to 

(2-96)  is  less  close  than  (2-97),  with  E(k,-r)  and  E^(k,co)  instead 

2  -7/3 

containing  factors  >S^sovo  -k)  an<*  +  (^/sovok)  J 
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treated  otherwise  by  Lin  (1953) ,  Uberoi  and  Corrsin  (1953) . 

We  define  the  hypothesis  here  in  terms  of  the  measured  and 
the  hypothetical  spectra  which  it  is  commonly  invoked  to  relate, 
rhe  hypothesis  approximates  the  one-dimensional  wavenumber 
(k^)  spectrum  of  a  turbulence  quantity  in  the  mean  rest  frame 
by  its  frequency  spectrum,  multiplied  by  uQ  at  a  fixed  point 
in  the  measurement  frame,  the  latter  frame  having  velocity 
-u  relative  to  the  former,  for  a  frequency  o)  =  uQk^.  We 
consider  first,  for  simplicity,  the  spectra  corresponding  to 
the  trace  of  the  velocity  correlation,  W(r,t),  rather  than  to 
the  longitudinal  or  transverse  components  individually.  The 
first  spectrum  of  convern  is  thus  the  one-dimensional  energy 
spectrum  in  the  rest  frame  (multiplied  here  by  -Jr) ,  given  by 

(2-98)  E(1)(k,)  =  (2T‘) ’  1/°°  drie  1  H* (?,,()) 

-00 

where  r^  =  (r^,0,0);  the  second  is  the  measurement -frame  analog 
of  the  trace  of  Eq .  (2-26)  with  F  =  0, 

(2-99)  0(0.)  =  9(0,0))  =(2>r)  Lf  dte1  W(0,i). 

-Oo 

where  0(r,x)  refers  to  the  frame  with  velocity  -u  -  (-u  ,0,0). 

We  wish  to  consider  the  departure  from  unity  of  the  "Taylor  ratio" 

(2-100)  T  =  uo0(uok1)/E(1)(k1)  . 
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Eq.  (2-98)  can  be  straightforwardly  transformed  to 
become,  as  given  by  Hinze  (1959),  Eqs .  (3-72) ,  (3-73); 

(2-101)  E^kj)  =  r  dkk_1E(k). 

kl 

In  Eq.  (3-19)  we  may  use  the  kinematic  relation 
(2-102)  W(7,t)  =  W(?  -  uqt,t). 

By  use  of  the  equivalent  of  (2-14)  written  for  the  mean  rest 
frame  and  of  (2-23),  Eq.  (2-99)  may  be  written 

(2-103)  6(o))  =  (27T)‘1/d3^P(v)/d3irk"2E4(k,o)  -  £.(oo  +  v)). 

Here  we  shall  find  the  Taylor  ratio  T  in  the  local- 
convection  approximation  (2-92)  and  consider  the  dispersive 
effect  in  Sec.  3.3.  From  (2-103)  and  (2-92)  we  find 

(2-104)  0(0))  ,  /d^POOju  +  v!"1/00  _  _  dkk_1E(k)  , 

o;/|uo  +  v| 

whence 

(2-105)  uo0(uQk1)  =  /d3vP(v)I(v/u  A)/”  dkk'lE(k), 

k1I(v/uo,A) 

where 

(2-106)  I(v/u0,A)  -  [1  +  2(v/uo)A  +  (v/uo)2j‘1/2, 
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o 

A  a  V'ir0/vu0#  and  / d'V  may  be  regarded  as  expressed  in  the 

spherical  form  /^dw^/^dfl/^d^ .  In  this  approximation  T  is 
o  o  -1 

given  as  the  ratio  of  (2-105)  to  (2-101)  .  From  (2-105)  and 

^TT 

(2-106),  we  note,  if  the  azimuthal  integral  /  d0p(v)  is 

o 

an  even  function  of  \  then  T  differs  from  unity  only  by  terms 
of  order  (v0/u0)^;  this  will  be  so,  in  particular,  if  P(v) 
is  isotropic. 

We  now  limit  consideration  to  a  power-law  wave- 
number  spectrum: 

(2-107)  E(k)  =  cnk"(n+1) 

In  actual  fact  we  are  concerned  with  the  inertial  subrange  case 
n  ~  2/3.  The  k-integrations  in  (2-101)  and  (2-105)  may  now 
be  performed  to  yield  a  ratio 

(2-108)  t  =  |d3vP(v)(|uo  +  v|/uQ)n  =  /d^vP^  [  I(v/uq  ,*)]  n. 

Form  (2-107)  implies  Tp(r .O^r11,  and  from  this,  proceeding 
from  the  original  forms  (2-98)  and  (2-99)  by  use  of  (2-54)  or 
recasting  (2-108), we  can  express  T  in  the  alternative  form 

(2-109)  T  =  ^(-uot, t)/^(-uqt ,0) ; 

the  right  side, under  the  present  assumptions,  is  properly 
independent  of  t.  For  n  =  2/3  and  an  isotropic  normal  P(v) 
recalling  (2-77)  and  (2-80)  (for  i  -  0),  we  may  write 
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A 


(2-HO)  T  =  G0(Vo/uo),  [1  +  s2(v0/u0)2]l/3[l  -  H0(v0/u0)  ]. 
for  vq/uo  «  1,  Eq.  (2-74)  then  yields  explicitly 
(2-111)  T  *•  1  +  (5/27)(v0/u0)2. 

More  generally,  for  arbitrary  vq/uo  we  can  compute  the  correction 
factor  T  explicitly  from  (2-110)  by  use  of  the  value  of  sq 
[Eq.  (2-83) J  and  the  graph  of  the  function  Hq(0o)  given  in 
Figure  1. 

We  can  proceed  with  the  spectra  corresponding  to 
individual  component  correlations  W^(r,t)  similarly  to  the 
above  for  the  trace  W(F,t).  Thus,  we  let  (k^)  be  defined 
by  (2-98)  with  W(r^,0)  replaced  by  W^^(r^,0),  (co)  by  (2-99) 

A  A 

with  W(0,t)  replaced  by  W^(0,t),  and,  analogously  to  (2-100), 
define 

(2-112)  Tjj  =  -^(u^/E^d^). 

In  particular,  as  at  (2-109),  in  the  local-convection  approxi¬ 
mation  with  a  power-law  spectrum  (2-107),  we  find 

(2-113)  Tij  »=  Vv-  'O/tf'j.j  (-UqT  ,0)  , 

which  once  more  is  properly  independent  of  t.  For  the  inertial 
subrange  and  an  isotropic  normal  P(v) ,  as  at  (2-110),  we  then 
have  as  Taylor's  ratio  for  the  longitudinal  (i  =  1)  and 
transverse  (i  -  2)  velocity  components 
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(2-114)  Tu  =  Gt(v0/u0)  =  [1  +  Si(vo/uo)2)1/3tl  -  H1(v()/uo)  ] 

(no  sum  on  i) ,  where  again  the  are  given  at  (2-83)  and 
the  Hi(v0/U0)  f°r  arbitrary  v/uq  in  Figure  1.  For  v/uq  «  1, 
by  (2-74),  as  at  (2-111),  these  correction  factors  are  given 
by 

(2-115)  Tu  .  1  +  (ll/27)(v0/uo)2,  T22  -  1  +  (11/108)  (vo/u0) 2 

It  is  possible  to  re-interpret  the  computed  correction 

factors.  Suppose  we  define  new  ratios  T1,  Tj^  obtained  by 

replacing  uQ  by  some  other  velocities  u',  u^ ,  respectively,  in 

the  definitions  (2-100)  and  (2-112) .  From  the  power-law  form 

(2-107)  we  readily  find  T’  =  (uo/u')nT,  =  (“o^ij^ij  * 

Since  the  T  and  proved  to  be  independent  of  k^,  we  can 

make  T*  =  T!,  =  1  by  choosing  u'  =  u  T^n,  u.  .  =  u  tH”. 

lj  o  ij  o  1J 

From  (2-114),  for  example,  we  see  that  we  can  provide  an  exact 
Taylor's  hypothesis  by  employing  effective  convection  velocities 

....  ,  2  2  2N  l/2r .  u  r  /  \  ,3/2 

(2-116)  uu  =  (uQ  +  Sivo)  [1  -  Hi(vo/uo)] 

to  compute  (k^)  -  (u^^k^)  .  This  velocity  differs 

for  longitudinal  and  transverse  components.  For  vo/uq  «  l, 
by  (2-115)  we  have 

(2-117)  ua  =  u0fl  +  (11/18) v2/u2  ] ,  u22  =  u0[l  +  (11/72)v2/u2]_ 
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As  for  the  dispersive  correction  to  the  local- 

convection  approximation  applied  here,  in  Sec.  3.3  it  is 

shown,  for  the  E^(k,<D)  considered,  that  dispersion  contributes 

to  Taylor's  ratio  an  added  term  of  the  order  — (eu”^k^)^^ 
2-2/3 

~(v0/u0)  (kiL)  and,  more  generally,  the  added  term  vanishes 
independently  of  vQ/uo  in  the  subject  limit. 

We  note  here  a  result  for  a  measurement  frame 
/\ 

quantity  generalizing  0(o>)  of  Eq .  (2-99),  namely  the  frequency 
transform  ^j(r,ou)  of  the  spatial  velocity  correlation  (cross- 
spectral  density) .  Analogous  to  the  mean-rest-frame  quantity 
0^j(r,ao)  of  Eq.  (2-26),  0^  is  given  by 

(2-118)  0tj(r,a>)  =  (2Tr)_1/“dTeiU,T  W ^  (r  -  uqt,t) 

or  similarly  by  the  second  form  of  (2-26)  with  cu  ->  cu  -  uq  •  Ic 
in  E^.  A  result  will  be  given  only  with  neglect  of  dispersion 
and  with  the  further  approximation  of  space-time  isotropy  for 

the  longitudinal  decorrelation  ^£(r,i).  The  strtamwise  element 

_  _ 

©ll(r,ci-j>  where  uQ  =  (uo,0,0);  is  then  found  from  (2-118)  and 
(2-89)  ,  or  from  (2-97)  and  the  alternate  form  of  (2-118) , 
to  be  given  by 

(2-119) 

®  1 1  (r »UJ)  =  ^  ^[T(l7/6)]  ^Nc^^ue^(ai/ue)  ^^^exp(iar^/uc ) 

xrS'1/6[“r8/ue)K5/6(u,rs/ue)4  ns /ue>  J. 
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where 


d2  2  2  2  2  2  n2 
R  =  r2  +  r3  ,  rfl  =  +  R  , 

(2-120)  ue2  =  u2  +  s2v2,  ve=  slVo/ue> 

2/  /i  2  2 .  2X 

u„  =  u  /u  =  u  (1  +  s,v  /uj . 

C  0  o  o  1  o  o 

Eq.  (2-119)  holds  for  all  v  /u  and  reduces  to  the  rest-frame 

oo 

_  ^ 

quantity  (^(r'.w)  for  vq  =  0.  For  vj uq  «  1,  the  contours 

A  _ 

of  constant  l®n(*>w)  |  are  elongated  along  the  longitudinal 

axis  relative  to  the  transverse  by  a  factor  -  u0/v0*  T^e 

effective  convection  velocity  u„  defined  by  writing  the  phase 
/\ 

of  0^  as  uJr^/uc  thus  proves  to  differ  from  uQ  as  given  by 
(2-120).  We  recall  from  the  discussion  of  Sec.  2.4  the 
limitation  that  when  torg/uc  is  large  (and  hence  l6ul  small) 
the  result  (2-119)  of  the  space-time  isotropic  approximation 

A 

will  conform  neither  to  the  true  0^  including  dispersion  nor 

to  the  exact  local-convection  approximation. 

At  this  point  we  may  usefully  discuss  the  relations 

among  variously  defined  convection  velocities.  Even  in  the 

present  unsheared  homogeneous  flow,  these  differ  from  the 

relative  rest-frame  velocity  uq  and  from  one  another,  though, 

2  2 

when  v  /u  <<  1,  onlv  by  terms  of  the  order  of  v  /u  .  In 
o  O  9  '  J  o  o 

7 - - 

We  note  the  result  of  the  exact  local -convection  approximation 
for  r  -  0:  0u(O,u)  =  (2/3)  1/3r"1/2r  (5/6)NV2/3v‘1  (o>/vo)~5/3. 

This  agrees  with  (2-119)  for  this  case  to  within  a  constant,  as 
it  must  on  dimensional  grounds. 


unsheared  flow,  if  the  measurement  frame  coincides  with  the 
laboratory  frame,  the  intensity  vo/uq  is  small  enough  that 
such  differences  are  unmeasurable.  In  some  types  of  intense 
sheared  turbulence,  or  in  an  envisioned  experiment  where  the 
measurement  frame  has  a  velocity  intermediate  between  that 
of  the  laboratory  frame  and  the  mean  rest  frame*,  on  the  other 
hand,  these  differences  may  be  appreciable.  In  any  event  it 
is  illuminating  to  examine  these  differences  in  the  present 
context  of  a  simple  unsheared  turbulence  with  fixed  uq. 

A  convection  velocity  is  defined  with  reference  to 
some  property  of  a  correlation  function  or  its  transform  in 
the  measurement  frame,  and  this  velocity  in  general  differs 
according  to  the  property  in  question  and  also  according  to 
the  quantity  the  correlation  refers  to,  e.g.  streamwise 
component  of  fluctuating  velocity  (along  uq) .  We  define 
convection  velocities  here  by  four  such  properties  and  refer 
for  definiteness  to  this  streamwise  component;  the  definitions 
ire  summarized  here  in  Table  1: 


Turbulence  measurements  in  a  moving  frame  have  recently  been 
performed  in  a  different  connection  by  Uzkan  (1965) . 


No. 

Defining  Condition 

Convection 
Velocity  Defined 

1 

SW11(r,t)/Sr1  =  0  (Wu  =  max) 

ri/T 

2 

d^i  i  (r ,  -r )  /5  t  =  0  (Wn  =  max) 

rl/T 

3 

9u(r»  »  exp(io)r1/uc)  |  011(r,o3)  | 

uc 

4 

En^kP  =  uiien^unki) 

U11 

Table  1.  Various  Definitions  of  Convection  Velocity 

In  general,  the  convection  velocities  so  defined  might  not  be 

constant  but  dependent  on  the  pertinent  variables.  On  a 

A  _ 

standard  plot  of  the  space-time  correlation  W^(r,t)  vs.  time 
delay  with  the  time-correlation  curve  for  each  streamwise 
separation  r^  displaced  along  the  time  axis  by  r^/u^  condition 
1  corresponds  to  the  envelope  curve,  whereas  condition  2 
corresponds  to  the  curve  joining  the  maxima  of  the  individua! 
curves  .  These  definitions  have  been  discussed  by  Wills  (1964) 
in  the  more  general  context  of  shear  flow,  along  with  other 
definitions  which,  in  the  present  case,  coalesce  with  these. 
Condition  3  defines  a  convection  velocity  based  on  the  phase 
of  a  cross-spectral  density,  as  done  for  example,  by  Corcos 
(1964).  Condition  4  [recall  Eqs.  (2-98)  and  (2-99)]  defines 
a  convection  velocity  based  on  Taylor's  hypothesis  relating 
spectra  in  frequency  and  wavenumber. 

We  give  results  for  the  convection  velocities  of 
Table  1  in  the  local-convection  approximation,  neglecting 
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dispersion,  for  the  inertial  subrange,  and  assume  once  more  an 
isotropic  normal  distribution  of  fluctuating  velocity.  We 
also  compare  the  results  given  by  the  space-time  isotropic 
approximation  to  V'f(r'«T)«  corresponding  to  neglect  of  H  ^  in 
Eq.  (2-80).* 

Since  W^(7,t)  =  W^(7  *  %T»T)  exactly,  the  convection 
velocity  of  condition  1  is  just  up  if  Wu(r',x)  has  a  maximum 
for  all  t  and  all  at  r^  -  0.  Examination  based  on  the 

results  (2-80)  [and  (2-74),  (2-75)]  for  ^(r’.-t)  and 
^g(r'x)  shows  that  in  fact  ^^(7',  t)  (obtained  from  (2-30)), 
as  well  as  and  f  individually,  have  minima  at  r|  =  0.  This 

result  is  given  also  by  the  space-time  isotropic  approximation 
to  ^(r'jT).  Hence  the  convection  velocity  of  condition  1 

is  just  u  . 

J  o 

With  regard  to  the  convection  velocity  of  condition  2, 
we  find  by  use  of  Eq.  (2-75)  and  (2-74),  respectively,  that 
this  is  given  in  two  opposite  limits  by 


(2-121) 


uQ[l  +  (5/3) Vq/u^ ]  for  (v0/uQ)2  +  (R/vot)2«  1, 
UQ[ 1  +  (11/9) v2/u2]  for  (v0/up)2  +  (R/vqt) 2  »  1, 


2  2  2 

apart  from  terms  of  higher  order,  where  R  =  r2  +  r^ .  In 
the  space-time  isotropic  approximation,  on  the  other  hand, 

*We  do  not  neglect  H2  to  find  ^(r'.-r)  but  accept  the  result 
implied  by  the  approximation  to  . 
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from  (2-84)  we  find  this  velocity  to  be  given  independently  of 

r2»r3 >  and  T  by  u0(1  +  8i  vo/uq)  '  uo^  +  1,455vo^uo^- 

The  convection  velocity  defined  by  condition  3 

has  not  been  determined  in  the  exact  local-convection 
approximation.  In  the  space-time  isotropic  approximation, 
however,  it  was  given  by  (2-120).  The  latter  is  the  same 
velocity  that  was  given  by  condition  2  in  this  approximation. 
More  generally,  in  any  approximation  in  which,  for  some  i,  j 
and  7, 

(2-122)  ©^(r.uj)  =  Aij(r,oj)  exp^u^/u^)  (A^  real) 

with  the  phase  convection  velocity  u£ij  independent  of  u),  and 
in  addition  Aij(r^,-ou)=  A^F.u)),  it  follows  from  the  general 
relation 

dWij(r,T)/dT  =-i/do)i'iCL1u5ij  (F.a,) 

that 

[^u(r,T)/3,)ri  =  ucljT  =  0, 

i.e.  in  such  case  the  convection  velocities  defined  by 
conditions  2  and  3  are  the  same. 

The  convection  velocity  defined  by  condition  4 
was  given  at  Eo  (2-116)  or^  for  v0/u0  «  1,  Eq.  (2-117).  In 
the  space- time  isotropic  approximation,  the  result  is  given  by 
omitting  in  (2-116);  for  vQ/uo  «  1,  this  becomes 
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Up  =  uq(1  +  sivq/uo)  "  uQ(l  +  0*728v^/u^).  The  convection 
velocity  defined  by  Taylor's  relation  (condition  4)  depends, 
we  note,  on  the  rest-frame  space-time  decorrelation  ^p(r',T) 
in  a  range  of  the  space-to-time  separation  ratio 
rjyvot[=  (r^-voT)/voT ]  opposite  to  that  which  determines  the 
velocity  defined  by  condition  2.  In  the  former  instance,  by 
(2-113),  the  relevant  ratio  is  | r^| /vqt  =  uQ/vo,  whereas  in 
the  latter,  by  (2-121),  it  is  r{/v0T ~  v0/u0 • 

We  note  that  in  the  exact  local-convection  approxi¬ 
mation  the  convection  velocities  defined  by  conditions  1  and 
4  are  constant  (i.e.,  independent  of  r^.r^jT,  and  of  kp 
respectively,  while  in  the  further  approximation  of  space-time 

isotropy  the  velocities  defined  by  all  four  are  constant.  The 

I  2  2\ 1/2 

results  are  summarized,  for  vq/uo  «  1  and  [r2  +  r^l  '  «  vqt^ 

in  Table  2,  which  gives  the  values  of  aQ  where 

2  2 

the  respective  velocities  are  written  as  uQ(l  +  aovo^uo^  * 


Defining 
Condition 
(Table  1) 

Approximation  Used 

Exact  local  convection 

Space-time  isotropic 
longitudinal  component 

1 

0 

0 

2 

5/3  (v0/u0<U,  R/vot«1) 

s?  =  1.455 

X 

3 

s \  -  1.455 

4 

11/18  (vo/uQ«l) 

s\/2  =  0.728  (vo/uo  <.  1) 

Table  2.  Coefficients  a  for  various  convection  velocities 

2  2° 

u  (1  +  a  v  /u  )  referring  to  the  streamwise  velocity 
o  o  o  o 

component . 
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3. 


RESULTS  WITH  VELOCITY  DISPERSION 


We  wish  to  consider  more  explicitly  the  departure  from 
the  local -convection  approximation  due  to  the  spatial  and  temporal 
dispersion  of  velocity  in  the  local  proper  frame.  This  step 
permits  us  to  consider  the  imporaant  spectral  regime  cu/v0k»l 
where  dispersion  is  the  determining  effect,  and  to  derive 
corrections  to  quantities  where  the  dominant  effect  is  convec¬ 
tion.  We  assume,  as  discussed  in  Sec.  2.1  and  the  Appendix, 
that  Eq.  (2-6)  with  an  intrinsic  decorrelation  having  the 
similarity  form  (2-42)  remains  competent  to  yield  these  correc¬ 
tions. 

3. 1  MODEL  INTRINSIC  VELOCITY  SPECTRUM  INCORPORATING 

DISPERSION - 

The  required  form  (2-42)  for  lf/( r,r)  in  the  xaertial 
subrange  implies  that  in  a  corresponding  domain  the  intrinsic 
wavenumber  frequency  spectrum  E^(k,a>)  has  the  functional  form 

(3-1)  \(k,co)  =  e1/3k'7/3G(n),  n  =  e"1/3k-2/3H 

The  inertial  subrange  for  !j/( r,-r)  is  defined  by 
(3-2)  L  »  (r2  +  c lx  I3) 1/2  »  lQ 

[see  condition  (2-39)  and  footnote  1,  preceding  page];  considera¬ 
tion  of  Eq.  (2-15)  indicates  that  the  conjugate  subrange  of  (3-1) 
for  TJ^(k,o>)  is 

(3-3)  L"1  «  (k2  +  c  _1  loot  3)  1/2  «  l'Ql  . 

The  form  (3-1)  insures  that  *^(r,T)  has  the  requisite  limiting 
forms  (2-46) ,  provided  the  zeroth  and  first  moments  of  G(H)  con¬ 


verge. 
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We  introduce  for  G(fi)  in  (  3-1)  a  conjectural  ex¬ 
plicit  form  that  meets  the  obvious  requirements  on  this  function 
and  has  vertain  other  plausible  properties  (see  Appendix); 
specifically,  we  limit  consideration  first  to  the  following 
class  of  functions  indexed  by  m: 

(3-4)  Gm(n)  -  am<r1/2  (l+ym2n3)‘m 

(mil),  where  a_  and  y  are  constants.  This  form,  with  non- 

vanishing  y  ,  correctly  provides  a  finite  spatial  scale  such 

that  motion  correlated  over  a  distance  is  associated  « 

i 

mainly  with  frequencies  smaller  than  the  similarity  estimate 
-ve^^k2^3.  Some  results  to  be  obtained  are  independent  of  the 
specific  form  of  (3-4),  apart  from  values  of  numerical  coef- 

•A 

ficients,  and  are  really  consequences  of  similarity  together 
with  presumed  asymptotic  properties.  The  type  of  function 
(3-4)  for  integral  m  has  been  chosen  to  lead  to  tractable 
integrations. 

In  illustrating  the  effects  of  dispersion  we  give 
results  on  assumption  that  m  ■  1  in  (3-4),  corresponding  to 

•J 

E^(k,ou)  independent  of  k  in  the  limit  n  ■*  00 ,  and  also  on 

•j  4 

assumption  that  m  -  3,  corresponding  to  E^(k,oo)  varying  as  k 
in  this  limit.  The  latter  assumption  is  favored  by  the 
following  argument.  The  restriction  (3-3)  on  the  domain  of 

Quantities  linearly  related  to  E^(k,u>)  can  be  obtained  for 
m-3  from  those  for  m-1  simply  by  differentiation. 
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the  inertial- subrange  form  (3-1)  suffices  even  when  ?gsL"\ 
at  least  in  the  case  where  even  the  large-scale  turbulence  is 
homogeneous  and  isotropic  so  that  Eq.  (2-15)  provides  a  valid 
basis  for  (3-3)  at  low  k.  (For  k^L~^,  condition  (3-3)  is  equiva¬ 
lent  to  l«(oiL/vo)  3/  )  Now,  in  order  to  have  E(k)  «  k^ 

for  k  as  required  in  this  case,  by  reference  to  (2-25) 

M  A  .1 

we  must  also  have  E^(k,aj)“k  for  k<«L  ;  presuming  that  this 
requirement  on  E(k)  cannot  be  relaxed  with  regard  to  contributions 
from  a>  »vo/Lt  we  must  have  fs^(k,u))«k^  for  sufficiently  small 
k  even  in  the  domain  of  (3-1).  This  result  implies  G(fi)«fi’^^2 

lx 

as  The  requirement  E^(k,a))ock  as  k-K)  is  likewise  implied 

-IQ  /9 

in  the  usual  way,  and  hence  G(fi)“f2  '  as  n-*<=  in  the  domain 
(3-3),  if  it  is  supposed  that  the  intrinsic  velocity  spectral 
tensor  (k,jj)  of  Eq.  (2-13)  should  be  analytic  at  k-0. 

This  limiting  dependence,  in  the  case  of  the  assumed  form  (3-4), 
implies  m«3. 

From  (3-4)  via  (2-15)  we  are  able  to  find  the  limiting 

1c 

forms  of  the  resulting  function  F(z)  in  the  decorrelation  (2-42): 

(3-5)  F(z)  ■*  Ao(l-tamz2)  as  z>0 

(3-6)  F(z)  -  Bz(l+cmz"3/2)  as  z^; 

Aq  and  B  are  regarded  as  given  independently  of  m,  and  in  the 
chosen  cases  m-1  and  m«3  we  obtain  for  the  coefficients  a^j  £m 

See  Appendix, Eqs .  (A-19)  and  (A-22).  In  this  and  other  cal¬ 
culations  below,  we  have  performed  various  quadratures  with 
the  aid  of  Bateman  (1954). 


4 


5Q 


(3-7)  a1  -  (8/3/5n)(B/Ao)2  c  0.881(B/Ao)2 
a3  -  [7(55)/(27)2]a1  -  0.529^ 

£l  -  (2n)'1/2c3/2(B/Ao)'3/2 
e 3  -  (8/3) (27/55)3/2ei  ~  0.9l6£l 


and  for  the  constants  a  and  -y  in  (3-4) 

m  'm  v  ' 


«!  .  n-2c3/2Ao(Ao/B)1/2,  -  c3(Ao/B)3, 

a3  -  7j2  -  (27/55)37l2, 


■  <  * 


I 

l 


where  c  -  5n/6r(l/3)  ~  0.977. 

2 

The  z  dependence  of  the  next- to- leading  term  in 
the  limit  (3-5)  occurs  not  only  for  the  choice  (3-4)  but  for 
any  G (ft)  that  decreases  more  rapidly  than  n  as  Q->®.  The  de¬ 
pendence  of  the  next-to-leading  term  in  (3-5)  is  indicated  in 
the  Appendix  to  be  plausible  also  on  independent  grounds. 

The  dependence  of  (3-4)  as  ft-K)  and  the  related  dependence  of 
(3-6)  as  z-*®  are  open  to  question  (see  Appendix, Eq.  (A- 21), 
for  a  generalization  of  (3-6).  These  dependences,  however, 
are  irrelevant  to  all  results  to  follow,  except  that  numerical 
coefficients  in  Eqs.(3-7)  and  those  following  for  am,  a^,  and 

y  would  be  somewhat  altered  if  the  dependence  in  this  limit  were 
m 

otherwise. 
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3.2 


VELOCITY  SPECTRUM  IN  THE  MEAN  REST  FRAME 


From  the  assumed  intrinsic  energy  spectrum  E4(k,o3)  we 
can  compute  by  (2-23)  the  spectrum  E4(k,oo)  for  the  mean  rest 
frame.  We  give  approximate  results  based  on  two  opposite 
limits.  For  the  moment,  we  assume  an  approximately  isotropic 
distribution  P(v)  .  Eqs.  (2-34)  and  (3-4)  then  yield 

o-8)  f h.-V)  r\,  i^vmi~1/2 

o  -1  (k^+7m^e  M^kvul  3)m 


Two  dimensionless  variables  are  pertinent  here,  namely  t)  k  a)/vQk 
-3  -1  -1 

and  £  *  evQ  k  ..  (kL)  ,  where  the  latter  estimate  follows  from 
3 

€«.vq/l.  We  may  consider  C  «  1,  since  we  are  considering 
the  inertial  subrange.  In  a  domain  where  T]  is  not  too  large, 
the  convective  approximation  represents  a  valid  lowest-order 
result;  assuming  a  normal  distribution  for  P(v) ,  inserting 
E(k)  in  (2-95) ,  and  denoting  this  convective  approximation  by 

E4c’  we  write 

(3-9)  E4c(k,cu)  -  (3/2n)1/2Nvo_1c2/3k’8/3exp[-(3/2)(tc/vok)2] 

with  N  *  [ 10/9r(l/3]Ao ,  in  agreement  with  (2-95).  By  analysis  of 
(3-8)  we  find  the  dispersive  correction  E^-E^  to  lowest  order  in  £. 
The  result  is  given  by 

(3-10)  E4(k,a3)=E4c(k,a)){l-(3/2)boain(cvo_3k-1)  U-3(*/vk)  ]  ) 

where  bQ-( 3/5) 1(1/3) /r( 2/3) *1 . 19  [cf. (3-7)).  As  we  should  expect, 


« 
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dispersion  thus  decreases  relative  to  E^c  for  small  a)/vok  and 
increases  it  for  larger  co/vQk.  In  a  domain  where  o>/vQk  is  large, 
as  discussed  previously,  E^  no  longer  is  related  to  E^c  but  is 
determined  mainly  by  dispersion.  From  (3-8)  in  this  instance,  to 
lowest  order  in  (a)/vQk)"^  we  obtain 

(3-11)  E4(k,o>)  ^ainc3/V7/2(ek2/a)3)m’1, 

the  terms  omitted  being  of  relative  order  (u)/vQk)  (We  con¬ 
sider  only  cosO,  since  E4(k -U))»E4(k,u>) . )  The  approximation  (3-10) 

2/3 

can  be  good  only  if  Q  '  «1  and  in  addition,  if  r&l,  the  value 
(3-10)  greatly  exceeds  (3-11).  Hence,  the  domain  of  validity  of 
(3-10)  may  be  given  roughly  as  that  where 

(3-12)  C2/3«l  if  n  <  1,  or 

(3-13)  (Crf3)  m"1  rr7/2C5/6exp(3r12/2)«l  if  n  >lj 

similarly  the  domain  of  validity  of  (3-11)  is  that  where 

(3-U)  T|»l  and  (Cri'3)m_1  T!"7/2C5/6exp(3Ti2/2)»l . 

The  result  (3-11)  corresponds  to  simply  neglecting  kvp 
22-13  -2-3 

relative  to  u)  and  then  k  /y  e  r  =  y  rj  Q  relative  to  unity  in 
the  integrand  of  (3-8).  Returning  to  the  basic  relation  (2-23), 
we  can  take  these  steps  more  generally  as  follows: 

(3-15)  E4(£,u))  =  /d3vP(v)E4(k,a)  -  v-fc)  ~  t4(k,u)) 

in  the  subject  domain.  The  result  (3-15)  or  (3-11)  is  evidently 
independent  of  the  assumption  of  an  isotropic  P(v) .  More  generally, 
assuming  only  that  E4(k,co)  varies  as  ^(“-1)  as  k-K)  at  fixed  a>/vo, 
and  recognizing  that  E4  in  this  limit  is  determined  by  dispersion, 
not  convection,  and  hence  becomes  independent  of  v  ,  we  can  infer 
the  result  (3-11)  (to  within  a  constant),  along  with  (3-15), 
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directly  on  dimensional  grounds.  Thus  the  result  does  not  depend 
on  the  assumed  explicit  form  (3-4)  for  E4(k,a))  or  on  (2-6). 

It  is  useful  to  extend  the  latter  considerations  also 

A  w 

to  the  spectrum,  say  E^(k,o>)  measured  in  a  frame  having  a  fixed 

_  A 

velocity  -uQ  relative  to  the  mean  rest  frame.  E4  is  given  by 
the  kinematic  relation  [Eq.  (2-3)] 

(3-16)  E4(k,o>)  =  e4(E,od  -  uD  •  IT)  . 

We  consider  the  domain  of  k,o>  where  conditions  (3-14)  are 
satisfied  with  oj  replaced  by  uo  -  uo,Tc*cu'>  so  that 
(3-11)  and  (3-15)  hold  for  E4(k,o>').  Further  restricting 
consideration  to  k  «  ^/uQ,  we  then  have  from  (3-16) 

(3-17)  ^(E.o.)  *  E4(kl<u)  -  E4(k,o,)  am7In'2m€3/V7/2(£k'-fc3)m-1, 


i.e.,  in  the  specified  regime  E4  is  nearly  independent  of  u 


3.3 


DISPERSIVE  CORRECTION  TO  TAYLOR 'S  HYPOTHES IS 


We  proceed  to  consider  the  effect  of  dispersion  on 
Taylor's  ratio  T,  defined  in  Sec.  2.5.  On  assumption  of  our 
standard  form  (3-4)  for  E4(k,u)),  we  may  write  the  basic 
equation  (2-103)  as 


(3-18)  |co-|u  +v  I  k  v 

9(u>)-a  el/\2<m-1>/d37p(V)/”dk/1dv  -7— j4t° 
m  -  -  1  (k  -  * 


-1/2 


+7  c 


10) 


—  TEST 

- 1  uq+v|  kvi  ) 


Likewise,  the  denominator  of  (2-100),  by  (2-47)  and  (2-101), 
is  given  by 
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(3-19)  E('l)(kx)  =  (a/SjNG2/^"573. 


The  ratio  T  is  a  function  of  two  dimensionless  variables, 

-3  -1 

vq/uo  and  eu^  k^  .  Initially,  we  regard  these  as  independent 

and  both  small.  The  departure  of  T  from  unity,  say  aT,  can 

then  be  approximated  as  the  sum  of  a  term  aTc  computed  by  neglect 

-3  -1 

of  dispersion,  with  euQ  k^  =  0  (i.e . ,  T  -  1  as  computed  above 
in  the  local-convection  approximation)  and  a  term  aT^  computed 
by  neglect  of  local-convection,  with  vQ/uo  =  0.  The  actual 
domain  of  validity  may  be  defined  in  the  light  of  the  results. 


To  find  aT^  we  take  vQ  =  0,  whence  (3-18)  yields 


CD-u  kv 

-2, 


-1/2 


(3-20)  3(a>)-a  e1 /2k2(m'1)/°°dk/1dv  T-rr. - -m 

o  -1  (k  +7  e  |U)-U  kv| J) 


-3  -1 

Analysis  of  this  expression  to  lowest  order  in  euG  k^  yields, 
with  some  effort  , 


(3-21)  &Td  =  (5/3)bni(B/Ao)2(€uo-3ki-l)2/3i 

where  b^^  and  b^  were  given  at  (3-10).  In  summary,  we  have 

(3-22)  T  -  1  =  ATc  +  nTd, 

with  aTj  given  in  the  limit  vq/uq  -*  0  and  to  lowest  order  in 
-3  -1 

eu  k^  by  (3-21);  from  (2-108),  we  also  have, 


If  we  consider  instead  T',  defined  by  (2-100)  with  u'  replacing 

- 3  -1,-2 

uQ,  then  uQ  in  (3-21)  is  replaced  by  uQ  u' 


in  the  limit  -*•  0, 

(3-23)  aTc  =  /d37p(7)(|uo  +  vt/uo)2/3  -  1. 

For  isotropic  normal  P(v) ,  moreover ,  by  (2-110)  this  becomes 

ATC  =  [1  +  *o(v0/u0)2]1/2[l  -  H0(v0/u0)]  -  1; 

if  also  vq/uo  «  1,  by  (2-111) 

AT"  (5/27)(v0/u0)2. 

-3  - 1 

Actually,  euo  k^  is  not  independent  of  v0/u0»  but  rather 
e%3k[l  ~(v0/u0)3(k1L)'1,  so  that,  by  (3-21)  , 

(3-24)  ATd-(vo/uo)2(k1L)‘2/3. 

Since  we  are  considering  only  the  regime  k^L  »  1,  we  have 
(3-25)  ATd/ATc^(k1L)"2/3  «  1. 

Hence  the  local-convection  approximation  to  this  extent  suffices, 

and  the  explicit  results  derived  above  on  this  basis  are  valid; 
furthermore,  these  hold  independently  of  the  magnitude  of  vo/uQ. 
In  particular,  by  the  choice  of  convection  velocities  given 
in  Sec.  2.5, Taylor's  hypothesis  becomes  exact  in  the  limit 
1/kjL  -  0. 


« 
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3 


DISPERSIVE  CORRECTION  TO  SPACE-TIME  VELOCITY 
CORRELATIONS 


We  consider  the  dispersive  corrections  to  the  space¬ 
time  decorrelations  ^(r,x)  computed  in  Sec.  2.3  by  the  local- 
convection  approximation.  For  this  purpose  we  assume  the 
standard  dispersive  intrinsic  spectrum  E^(k,u>)  of  (3-4).  In 
the  limits  we  shall  consider,  however,  only  the  limiting  re¬ 
sult  (3-5)  is  actually  pertinent.  Hence  the  results  are  inde¬ 
pendent  of  the  specific  form  of  (3-4)  ,  apart  from  the  relation 

(3-7)  of  the  constant  a  to  b/A  ,  and  are  contingent  sub- 

m  o 

stantially  only  on  the  expected  dependence  F(z)-F(0)*z'  as  z-*0. 

The  trace  function  r,x)  for  the  inertial  subrange 
is  given  by  Eq.  (2-51)  .  We  assume  an  isotropic  P(v)  ,  whence 

(3-26)  Hr,-)  -  j(cr)2/3  J"  dv  v2P(v)  jf1dliT!2/3F(*il'2/3)  . 

where 

i,  =  [l+(vx/r)2-  2(vr/r)M31/.?  z  =  £1/3r'2/;Vl  . 

Eq. (3-26)  is  the  equivalent  of  (2-69),  but  the  dispersive  factor 

oQ(z)  in  (2-56)  is  no  longer  replaced  by  unity.  We  assume  again 

2/3 

a  normal  P(v) .  The  ratio  ^(r,x)/(£r)  '  depends  on  the  two 

variables  0  =v  x/r  and  z,  or  equivalently  on  8  and 
o  o  o 

(3-27)  ’'■Bo1*3'2  -  (ct)1/2/v0  . 

We  have  r(»  i/l)*^«l  by  (2-5);  X  provides  an  appropriate 
r- independent  measure  of  the  effect  of  dispersion. 

The  limit  0  •*  »  is  of  particular  interest,  C<ncf, 
for  v  /u  «1,  it  determines  the  envelope  of  the  space-time 

A 

correlation  curves  W(r,x)  in  a  measurement  frame  having  velocity 
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-u0  (for  r  and  uQ  nearly  parallel).  By  letting  0o+«  and  then 

X-*0,  i.e,  ,  considering  the  limit  where  X-*0,  X0osz3/2— , 

we  are  able  to  evaluate  (3-26),  by  reference  to  (3-5)  and  (3-6), 

to  the  two  lowest  orders  in  with  the  two  coefficients  each 

o 

evaluated  to  the  two  lowest  orders  in  X.  Explicitly,  we  find 

(3-28)  »(r,T)  -  S0(€V0T)2/3|l+(l/3)e^+d0amxi,/3[l-(l/3)p;2lJ, 

where  dQ  *  (5/3)  6"  ^3b2  »  1.30  and  SQ  was  given  at  (2-76)*. 
Likewise,  by  considering  the  limit  0q-K)  at  arbitrary  X(<1) ,  we 
find  from  (3-26)  and  (3-5) 

(3-29)  t(r,T)  -  A0(Sr)2/3|l-«/27)p2+a111),4/3e*/3!l-(l/27)p2ij) 

in  which  the  two  parts  proportional  to  X°  and  X^3  are  evalua¬ 
ted  to  the  lowest  two  orders  in  PQ. 


Since  the  comparison  with  experiment  considered 
here  concerns  the  longitudinal  decorrelation  ^(r,T),  we  wish 
to  obtain  this  also  in  approximations  corresponding  to  (3-28) 
and  (3~29)  by  reference  to  the  continuity  equation.  The 
latter  is  still  expressed  by  Eq.  (2-78)  in  which,  however, 
now  depends  on  X  as  well  as  on  0Q.  We  find 

(3-30)  ^f(r,i)=S1(€v0T)2/3{l+(l/5)s;2+dociI|X4- '[  1-<1  /5)«0‘  2)} 

(x«l,  X0O»1), 

where  was  given  at  (2-76);  in  the  other  limit 

*If,  contrary  to  (3-5)  ,  1>( r,T)  contained  an  additive  component 
B' e|T|for  all  t,  i/( r,x)  would  likewise  contain  this  same 
r- independent  component.  67 

1 


(3-31)  \)f(r,T)  =  A(€r)2/3{l+(ll/27)S^+(ll/7)amX4/3By3 

x[l-(7/27)e2]}  (80  «  1). 

The  transverse  decorrelation  (r,r)  can  ^  obtained  similarly. 

© 

Writing  ^  where  ^c(r,T)  is  the  result  of  the  local 

convection  approximation  and  A^(r,T)  the  dispersive  correction, 

however,  we  may  insert  the  value  of  ^  given  by  Eq.  (2-80)  or 

(2-84)  without  resort  to  the  expansion  of  this  part  in  (3-  30) 

and  (3-31)  in  powers  of  3“^  and  0Q,  and  use  only  the  terms 

4/3 

proportional  to  X  in  these  equations  as  the  expansion  of 
A^  alone. 


We  may  consider  the  effect  of  the  dispersive  correc¬ 
tion  computed  here  on  the  comparison  with  experiment  in  Section 
2.3  and  Figure  2.  For  small  angle  0  between  r  and  uq,  the 

/S  _  A 

dispersive  correction  Av/£^(r,T)  to  as  given  by  (2-89)  or 
(2-90)  may  be  identified  with  A^(|r-uoT |  ,t)  .  We  approximate 


the  latter  for  3Q>land  30<1»  respectively,  by  the  relevant 
terms  in  (3-30)  and  (3-31),  where  3q  r  vDT/|r-voT|;  this 
approximation  is  unjustified  for  r,x  such  that  3q>1  butA3o<l, 
but  this  domain  is  relatively  small,  and  the  envelope  curve, 
in  particular,  will  be  given  correctly.  The  additional  parameter 
B/Aq  now  enters  via  am  [recall  Eqs.  (3-7),  (2-44),  (2-45)]. 

The  magnitude  of  the  quantity  X  of  (3-27)  in  the 
experiment  in  question  is  given  as  a  function  of  the 
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dimensionless  time  n  s  uqt/M,  where  M  is  the  mesh  spacing, 
according  to  (2-68)  by  X  =  (5/3)*^  (yo/uQ)Ke^^Ke^  where 

R«m  *  uqM/v.  For  the  parameters  of  the  particular  experiment, 

assuming  v0/uQ  =  0.027,  we  find  X  =  0.116n^^,  which  is  moderately 

small,  as  required  for  such  n  that  the  correlation  is  high. 

The  space -time  correlation  computed  by  (3-30)  and  (3-31) 

for  n  -2.31  and  A  -  0.55  on  assumption  that  0-0  and  v  /u  -0.027  is 
m  o  o 

compared  in  Figure  3  with  the  experimental  results  shown  previously 
in  Figure  2.  The  value  used  for  am  corresponds  to  B/Ao*1.62  if  m-1 
or  B/Aq-2. 23  if  m-3  in  Eq.  (3-4).  Also  shown  in  Figure  3  (by  solid 
lines)  is  the  result  obtained  by  supposing  that  the  disper¬ 
sive  term  A^  is  not  that  derived  above  but  simply  that 
corresponding  to  the  space-independent  (wavenumber- singular) 
form  ^(r ,t) -^<r,0)  =  Be|T|,  whence  A^f(r,r)  =  (1/3)B e 1 t I ;  the 
result  shown  corresponds  to  B/Aq  *  1.04  with  other  parameters 
as  for  the  other  computation.  There  is  little  to  choose  between 
the  computed  curves.  The  envelope  in  either  case  could  be  made 
to  decrease  more  rapidly  by  increasing  the  assumed  value  of 
B/Aq.  Figure  3  shows  also  the  envelope  curve  (a),  given  previously 
in  Figure  2,  corresponding  to  {jf  =  0,  vQ/uo  =  0.027,  but  without 
dispersion  (B/Aq-0).  For  likely  values  of  B/Aq  the  dispersive  ef¬ 
fect  is  seen  to  influence  the  space-time  correlations  appreciably 
at  the  Reynolds  numbers  of  the  experiments  considered.*  Inclusion 
of  dispersion  necessarily  worsens  agreement  with  experiment  for  the 
envelope  at  the  larger  time  delays,  since  the  inertial- subrange  form 
assumed  already  gives  too  small  a  correlation  at  such  t  or  |r-uoT| 

that  the  correlation  is  far  from  unity. 

w - 

Inclusion  of  dispersion  does  not  affect  the  computed  ratio  p  /p 

-  ■*»  c  s 

pertinent  to  shapes  of  isocorrelation  curves  discussed  in  Sec.  2.3. 


'i 
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4. 


PRESSURE  CORRELATIONS  AND  SPECTRA 


We  turn  to  the  subject  of  pressure  fluctuations. 
Pressure  is  related  to  the  velocity  derivatives  for  an  incom¬ 
pressible  flow  by  the  familiar  Poisson  equation 

(4-1)  V2p(x,t)  =  -pd^v^/aoc^Xj, 

where  v^(x,t)  refers  to  the  total  velocity  or,  in  the  unsheared 
flows  being  considered,  the  fluctuating  velocity.  Prior  to 
the  common  pursuit  of  consequent  relations,  we  progress  by  simple 
considerations. 

4.1  PRESSURE  SPECTRA  IN  A  MOVING  FRAME  BY  KINEMATICS 

AND  SIMILARITY  FOR  THE  INERTIAL  SUBRANGE _ 

We  consider  once  more  a  frame  having  constant  velo¬ 
city  -u  relative  to  the  mean  rest  frame  and  refer  to  the 
o 

pressure  on  a  planar  area  whose  normal  is  orthogonal  to  uq. 

Let  the  spectral  density  of  pressure  in  frequency  oj  and  two- 
component  wave  number  K  =  (k^,k^)  in  the  plane  containing  this 

A 

area  be  denoted  by  P(R,cu)  .  In  the  instance  of  the  rest  frame 

»»  ic 

(uQ=0)  ,  we  denote  this  spectrum  by  P(R,cu).  We  then  have 

(4-2)  P(K,u>)  =  P(R,o>  -  uo-K)  . 

We  assume  that  the  large  eddies  are  roughly  irotropic, 


In  the  terminology  of  Chandiramani  (1965)  and  others,  P(K,u.) 
(nore  precisely  its  wavenumber  integral)  is  the  fixed -transducer 
spectrum  and  P(R,cjo)  the  moving-axis  spectrum. 


so  that  the  local  fluctuating  velocity  Is  adequately  charac¬ 
terized  by  vQ,  and  limit  consideration  to  the  inertial  subrange 
L’1«(K2+ui2/u^)1/2«jt’1  (u0*vQ),  where  K=  iRl  .  P(R,cu) ,  as  a 
function  of  k^,  ,  at  fixed  od,  has  a  peak,  or  rather  a  ridge, 

in  region  of  R  where  the  probability  of  a  total  fluid  velocity, 
uq+  v,  that  satisfies  the  convective  condition  o>-  <v  =0 

is  substantial;  this  region  is  given  by 

(4-3)  la)-uQ*  Rl  <  vQ  K. 

We  recall  in  the  parallel  instance  of  the  mean-rest- 

frame  velocity  spectrum  E^(k,oa)  that  the  result  of  the  local- 

convection  approximation,  E^c(k,a))  [see(3-9)],  was  adequate 

[see ((3-10)]  if  o)/vQk<l  (or  even  somewhat  larger,  according 

to  (3-13),  if  C  is  sufficiently  small).  Similarly  here,  we 

expect  that  a  local -convective  approximation  to  the  mean-rest- 

frame  pressure  spectrum  P(R,cu')  will  be  adequate  roughly  if 

<u*  /vqK<  1.  Now,  the  pressure  is  quadratically  related  to  the 

velocity  fluctuations,  and  the  velocity  spectrum  (or  correla- 

2/3 

tion)  in  the  inertial  subrange  is  proportional  to  e  .  Hence 
P(R,o>')  in  the  local -convection  domain,  being  otherwise  inde¬ 
pendent  of  velocity  dispersion  and  thus  of  •,  must  have  the 
form 

P(R,cd')  =  p2€4/3D(vo,R,u>'), 

where  D  is  a  function  having  dimensions  (length)^^  (time). 
Since  D  must  also  be  a  properly  convariant  scalar  function  of 
the  vector  R,  we  must  have 


0 
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D(v,K,o)')  -  v^1  K'13/3  F(voR/co'); 

on  account  of  isotropy,  however,  F  can  depend  on  R  only  via 
K,  say  F  *  0(co'/voK).  Thus,  in  the  subject  domain,  we  have 

(4-4)  P(R,ou')  =  p2e4^3  v;1  K-13/3  0(cu'/voK),  whence 

(4-5)  P(R,oj)  =  p2€4/3  v”1  K'13/3  0([co-ao-K]/voK) 

We  now  assume  u  »v  ,  choose  coordinates  such  that 
o  o 

uQ  =  (uQ,0,0),  and  consider  the  width  Ak^  in  of  the  convective 

A 

ridge  in  P(R,  oj)  centered  at  k^^co/u^.  According  to  (4-3)  ,  we 

have  Ak,  '"(v  /u  )K.  The  factor  K-13^3  in  (4-5)  with 
loo 

l^^co/u^  2+k2  effectively  cuts  off  contributions  from 
k^>o3/uo.  Hence  A^<  V2 (vQ/uo)  (oa/uo)  for  all  significant  k^, 
so  that  Ak1/k^~vQ/uo«l.  Hence,  when  multiplied  by  a  power 
or  other  sufficiently  smooth  function  of  K,  the  function  0  be 
behaves  roughly  as  a  6 -function: 

(4-6)  0  ( [aj-u0k^  ]  /vQk)  «►  a^(v0/u0)KB(k1-a)/u0), 

where  ad  =  /“  dx0(x) .  This  approximation  equivalent  to 

W  -oo 

Taylor's  hypothesis  (based  on  the  unmodified  convection  velo¬ 
city  uQ)  . 

We  may  now  consider  the  frequency  spectrum  of  the 
point  pressure, 

(4-7)  P(oj)  =  /d2RP(R,co)  , 

where  the  integration  runs  over  the  entire  R-plane.  Since  in 

A  —  ^ 

this  case  P(R,uo)  is  weighted  uniformly  in  R,  P(oo)  is  largely 
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n> 


f 


determined  by  the  domain  of  R  which  contains  the  peak  in  P(R,u}) 
and  where  (4-5)  holds.  From  (4-5)  to  (4-7)  we  readily  obtain 

(4-8)  P(o>)  -  c5/3a^p2€4/3u’1(a)/uo)‘7/3t 

where 

(4-9)  c  =  2  /“  dx(x2+l)"n  -  7r1/2[r(n-l/2)/r(n)  1. 

11  o 

We  thus  infer  uniquely  the  frequency  dependence,  P(uj)«u)"7^3. 

The  conditions  for  (4-8)  are  vo/uo«l,  u)L/uo»l,  and  ujfo/vo«l. 

A 

We  consider  now  the  spectrum,  Q(u>)  ,  of  average 
pressure  on  a  circular  area  of  radius  Rq,  given  by 

(4-10)  Q (  ')  =  /  d2R[2J1(KR0)/KR0]2P(K,u))  . 

This  quantity,  for  the  homogeneous  turbulence  assumed,  is  not 
actually  accessible  to  experiment,  since  the  introduction  of  a 
pressure  sensor  of  the  prescribed  area  would  naturally  disturb 
the  homogeneity  of  the  flow  near  the  surface  of  the  sensor. 

We  form  the  quantity  rather  as  an  edifying  analog  for  the 
spectrum  of  average  pressure  on  a  similar  area  of  a  wall  in 
the  more  difficult  case  of  a  turbulent  boundary  layer. 

Suppose  first  ajRQ/uo<jC  1,  so  that  we  may  approximate 
the  area -averaging  factor  by 

[2J1(KR0)/KR0]2  =  1  -(1/4KKR0)2# 
y - 

By  (2-40),  the  last  condition  may  be  written 
(uuv/v2)  (vqL/v)  1/4«1. 

7  J 


Then  (4-5),  (4-6),  (4-10)  and  (4-8)  yield 


(4-11)  Q(u>)  =  P(aO[l-(a;Ro/uor]. 

Consider  now  the  opposite  limit  a)RQ/uo»  tt  .  In 

A 

this  case  the  range  of  K  near  the  peak  of  P(K,cjd)  in  (4-10), 
namely  K~u>/uo,  is  well  separated  from  the  range  K<£rrR~^  where 
the  area-averaging  factor  is  relatively  large.  We  must  examine 

A  A 

the  contribution  to  Q(u))  from  the  latter  range,  say  Q_(u)),  as 

A 

well  as  that  from  the  former,  say  Q+(co). 

A 

In  the  range  of  K  pertinent  to  Q_,  since  K«cjo/uq, 

A 

we  may  approximate  P(K,cu)  by  the  rest-frame  spectrum  P(R,u)). 

Form  (4-4)  for  P(R,ui)  may  not  hold  for  the  range  of  K  where 

K«o)/v0,  depending  on  the  role  of  dispersion.  At  the  same 

time,  we  have  no  justification  for  presuming  that  P(R,u>)  here 

does  not  depend  explicitly  on  vQ  and  is  independent  of  K,  as 

for  m  =  1  in  the  instance  of  the  form  (3-11)  for  the  velocity 

spectrum  E^(k,u>)(if  these  conditions  held  we  could  uniquely 
9  9  -6 

infer  P(K,u))  ~p  e  u>  ).  In  Sec.  4.2  it  is  indicated  that  in  this 
domain,  even  though  a)/vo K»l,  P(R,x)  in  fact  is  likely  determined  by 
local  convection  and  assumes  the  K-independent  form  of  (4-14): 

(4-12) 


D/v  ^  -a  2  4/3  10/3  -13/3 
P(K,u>)  ApP  c  vq  a)  , 


where  Ap  is  a  constant  of  the  order  of  unity.  The  range  of  K  where 
(4-12)  applies  then  yields,  in  the  limit  where  a>R  /vo>-»l,  the  largest 
contribution  to  the  integral  of  form  (4-10)  for  Q_(.o)  in  the  sura 
(4-13)  Q(u3)  s  §_(u>)+$+(*). 

From  (4-12)  we  thus  obtain 

(4-14)  Q_(o>)  ~4ttApc2C4/3vo10/V13/3ro‘2 

for  ajRo/uo>>l,  o)Ro/vo>>l,  coio/vQ<<l. 
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To  find  §+  we  may  assume  for  P(7,uu)  in  the  pertinent 
region  of  K  the  earlier  form  (4-5)  and,  without  appreciable  error, 
extend  the  integral  of  the  form  (4-10)  for  over  all  K.  Treat¬ 
ing  separately  two  opposite  limiting  cases,  (idR0/uo)(yo/uo)»tt 
and  «7t  ,  using  (4-6),  and  at  the  proper  stages  making  use  of  the 
result  (4/tt)  (KRo)”^  for  the  average  of  [  2J^(KRq) /KRq]  ^  over  an  in¬ 
terval  AK  such  that  n<<'AK«K,  we  find  both  cases 

(4-15;  q»  =  *-‘1e19/«V2«*/: '“'V 
with  cn  given  at  (4-9).  Conditions  for  (4-15)  and  (4-11)  (apart 
from  those  on  wRo/uQ)  are  as  given  for  Eq.  (4-8). 

varies  as  R0  ,  corresponding  to  a  nonvanishing 
effective  area  scale,  whereas  varies  as  R~^  .  Their  ratio,  by 
(4-14)  and  (4-15),  is 

(4-16)  <L/Q+  -  «o(vo/uo)10/3(^o/»o), 

where  aQ-*T  Ap/ci9/6a^‘  As  ^or  t*ie  reduction  factor  due  to  area 
averaging,  by  (4-8)  we  have 
(4-17)  Q(u))/P(o))  =  q_(o))  +  q  +  (u>), 

=  i:aJvo/uo)10/3<»Ro/uo)'2, 

mV'  :  M“Vuo>'3- 

where  a_  -  **A,/c5/3a0  and  a+  =  4c19/6/ttc5/3  -±  0.795. 

For  arbitrary  values  of  u)Ro/U0,  we  may  infer  for  $+(oj)  , 
from  (4-5)  and  (4-10),  a  functional  form 

(4-18)  Q+(o>)  -  o2€4/3uo-1Ro7/3h(vo/ur,a.R0/uo). 

Analogously,  for  boundary-layer  turbulence  we  may  conjecture 
that  q_/Q+  ~(v*/UJ2(u)R0/Ut>>. 
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By  (4-8)  and  (4-15),  limiting  forms  of  the  function  h(x,y) 

_  1  /o 

are  given  by  h(x,y)->  const.  X  y  '  J  as  x-*0,  y-*0  and  h(x,y)  — *  const. 
X  y  ~M>/ 3  ag  x_*o,  y-*®.  Taking  account  of  the  relation 
e«»vo  /L,  we  may  make  explicit  all  velocity  dependence  in  the 
scaling  form  (4-18)  by  writing 

Q+(cd)~P  2R0uq  3  (oi  /uq  )  ’  4  7  3w(  vq  /uq  ,o>Ro  /uq  )  , 
where  w(x,y)  is  a  function  related  to  h(x,y). 
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4.2 


PRESSURE  SPECTRA  AND  CORRELATIONS  BY  THE  MODEL 
OF  QUASINORMAL  VELOCITY  DISTRIBUTIONS _ 


The  results  of  the  preceding  section  can  be  clarified 
and  extended  In  part,  at  the  expense  of  introducing  imprecisely 
assessed  approximations,  by  pursuing  an  explicit  calculable 

A 

model  to  give  P(R,cd)  .  The  model  to  be  considered  will  have 
approximate  validity  only  in  the  convective  range  (4-3)  or 
somewhat  beyond,  and  not  in  the  dispersive  range,  in  particular 
not  where  (4-13)  applies. 

From  the  basic  relation  (4-1)  of  the  pressure  field 
to  its  velocity-derivative  sources,  the  rest-frame  spectrum 
P(K,cu)  is  found  to  be  expressible  in  the  present  case  of 
stationary  infinite  homogeneous  turbulence  as  the  following 
double  integral  over  the  coordinate  normal  to  the  plane  of 

X  (e.g,  from  Kraichnan  1956): 

(4-19)  P(K,w)  =  (l/4)p2K-2  /”  dx2  /“  dx2  exp[-X(|x2  1+ |x2  |)  ] 


where 

(4-20) 


xS(x2  -  x2,K,cd), 

S(C2.R^)  =  (2tt)  *3/  d2(/  dTexp[-i(R-C  -  cur)  ]Q('2, !,t)  , 
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(4-21)  Q(C2,C>t)  -  [  VeW  <<V'@vlv;>-<  W<vivv>>  ] .  t  _  .  , 

X  ®x  ^  f 


va  "  va<*»t>’  vi  "  va(5,»t  +  T)  »  *  =  (x1,x2’x3)»  f  =  ^l^2»f*3> 
£  “  (C i» C3)  ,  dQ  =  b/dxat  da'  =  d/dx^,  and  integrals  extend  over 

the  infinite  domain  unless  otherwise  specified,  P(R,co)  as 
given  by  (4-19)  -  (4-21)  is  so  normalized  that  the  mean  squared 
fluctuating  pressure  is  given  by 

(4-22)  <  >  -  /“ do)P(a))  =  /“doj/d2KP(K,o))  . 

-00  -00 


We  now  assume  the  turbulence  is  isotropic;  regarding  the  iner¬ 
tial  subrange,  this  assumption  corresponds  again  to  isotropic 
P(v)  in  (2-6).  In  this  event  we  may  write  Q  of  (4-21)  as  Q(r,T). 
We  define  L(k,m)  as  the  wavenumber -frequency  spectrum  of  Q: 

(4-23)  L(k,w)  *  /d^f/dxexp[-i(lc- r  -  u)  t)  ]Q(r,x) 


[£  =  ,k2  ,k^)  ] .  Then  S  in  (4-19)  may  be  written  as  a  transform 

of  L  over  k2«  In  the  present  instance  (unlike  that  of  boundary - 
layer  turbulence)  we  can  trivially  integrate  over  X2  and  x£  in 
(4-19)  to  obtain 

(4-24)  P(K.co)  =  (2^'4p2;dk2k'4L(k,x)  , 

in  which, we  recall,  k2  =  k^  +  X2. 

To  proceed  further  we  make  the  common  assumption  of 
a  quasinormal  velocity  distribution  with  reference  to  the 
two-point,  four -component  space-time  correlation  <v^vpv’vy>  in 
Eq.  (4-21)  for  Q(r,v)  ,  i.e.  we  assume  that  the  fc  irth-ordei 
moments  are  related  to  the  second-order  moments  in  the  same 
way  as  if  the  statistical  distribution  of  velocity  in  space-time 
were  a  joint  normal  distribution  (Millionshtchikov  1944, 


4 


<* 


'  • 


1 


\ 
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Batchelor  1951).  This  assumption  has  received  serious  criticism 
(Rraichnan  1957,  1959a)  as  the  basis  for  a  dynamic  treatment 
of  turbulence.  Nevertheless,  its  use  appears  unlikely  to  lead 
to  major  error  in  the  less  demanding  application  to  a  merely 
descriptive  treatment  of  velocity  and  pressure  correlations, 
and  the  assumption  has  often  been  so  used.  It  was  first 
employed  in  the  context  of  pressure  fluctuations  by  Obukhov 
(1949).  From  this  assumption,  Batchelor  (1951)  derived  the 
form  of  the  desired  correlation  function  Q(r,x)  for  vanishing 
time  separation,  i.e.  Q(r,0).  If  the  assumption  of  quasinor- 
mality  is  made  also  for  x  i  0  (Chandrasekhar  1955)  ,  the  deriva¬ 
tion  applies  equally  when  the  two  points  in  question  are 
separated  in  time  as  well  as  space.  In  terms  of  the  (mean -rest - 
frame)  longitudinal  velocity  decorrelation  ^(r.x)  the  relation 
derived  is 

(4-25)  Q(r,T)  -  4[2<^)2+  +10r\V{+  3r'2(*’)2), 

where  a  prime  denotes  d/dr.  Alternatively,  the  spatial  transform 
of  Q(r,x)  defined  by 

q(K  »t)  =  (2tt) '3/d3re~ik  rQ(r,T)  =  j  n  2  k"1  fQ  drrsinkrQ(r.x) 

is  related  to  the  corresponding  transform  E(k,x)  of  f(r?  ) 

(see  Eq.  (2-17))  by 

(4-26)  q(k,x)-(8TT2)"1k^/d3C'E(k,,x)E(|ft  -  £'  |  ,x)sin^0'  !£-£'  1 

where  cose*  =  IMc/kk'  .  L(k,o>)  of  (4-23)  is  then  given  as  a 

time  transform  of  q(k,x).  By  virtue  of  its  definition  (-4-21) 
as  a  fourth-order  derivative,  Q(r,x)  ,  as  shown  for  x  =  0  by 
Batchelor  (1951)  ,  corresponds  to  a  q(k,x)  and  hence  to  a 
spectrum  L(k,o))  that  vary,  for  any  reasonable  function  ^(rfT), 


as  k4  for  k«l,  and  likewise  in  the  case  where  these  quantities 
are  generalized  to  anisotropy  at  low  wave  numbers.  From  (4-24), 

(4-26)  and  (2-22)  we  can  express  P(K,o>)  also  in  terms  of  E^(k,u>)  : 

(4-26a)  P(R,o>)  =  (l/8n2)p 2J  dk2  J  d3fc'  sin4®'  |  I'4 

00 

x  J'  du)'E4(k',uit)E4(|lc-’R'|  ,u)-uu 1 ) 

-  CO 

In  principle,  in  the  inertial  subrange  we  can  assume 
our  standard  form  (3-4)  for  the  intrinsic  energy  spectrum  £4vk,uu) 

(suitably  extended  to  low  k  or  terminated  at  k-L  ^),  write  the  re-  *♦ 

suiting  E4(k,uu)  as  in  (3-8),  or  more  generally  (2-34),  and  compute 
the  desired  P(K,uu)  from  (2-26a)  or  else  proceed  via  (4-25)  or  (4-26). 

We  shall  deal  separately  with  the  domain  uu/voK»l  and  the  opposite 
domain  where  dispersion  may  be  neglected,  starting  with  the  latter. 

In  the  non-dispersive  domain  we  can  simply  approximate 
E(k,T)  by  (2-96)  or  Y^(t,t)  by  (2-80).  On  account  of  the  still 
considerable  analytic  difficulty,  however,  and  the  uncertain  error 
already  introduced  by  the  assumption  of  quasinormality,  we  proceed 
still  more  crudely.  We  approximate  Y^(r,T)  in  (4-25)  as  a  space-time 
isotropic  function.  The  consequent  error  in  the  computed  P(K,uu)  for 
the  inertial  subrance  relative  to  an  exact  local-convective  approx¬ 
imation  will  parallel  the  difference  previously  noted  between  (2-96) 
and  (2-92)  in  the  case  of  the  energy  spectrum  E4(k,uu);  furthermore, 
dispersion  will  limit  the  validity  even  of  an  exact  local-convection 
approximation  for  P(K,<u)  in  a  parallel  manner  to  that  discussed  for 
E4(k,oi>)  in  Sec.  3.2. 

In  the  stated  approximation,  we  have 
t^(r>T)  "  ♦|(8^>0)  s  i  i ( ®^i) 


(4-27) 
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with  Rj  of  the  form  (2-81).  In  the  inertial  subrange,  where, 
as  at  (2-84), 

(4-28)  ^1(R1)  -  A(€R1)2/3  (iQ«  Rx«  L) 

with  R^  given  specifically  by  (2-81)  and  (2-83),  we  find  from 
(4-25) 

(4-29)  Q(r ,t)  -  (16/9) A2€A/3R^8/3[15-(80/3)r2/R2+(112/9)r4/Rj] . 

In  the  viscous  subrange,  where,  as  at  (2-86), 

*1<R1>  -  <6/30v>Rol  <R01  <  V 

With  Rq^  given  by  (2-87)  and  (2-88),  we  find 

(4-30)  Q(r ,t)  =  4(€/15v)2(14  +  r4/R^) . 

We  note  further  the  result  where  a  modified  von  Karman  inter¬ 
polation  form  (Hinze  1959,  Eq.  (3-131)) 

(4-31)  ^1(R1)=  (V3)v2[l-exp(-koR1)] 

is  assumed,  in  which  k  is  a  constant  inverse  correlation 
length: 

(4-32)  Q(r ,t)  -  (4/9)v4exP(-2koRpk2[4k2(rVpJ) -lOk^J1 

x  (r2/R2)  (2-r2/R2>+R“2(15-20r2/R2+8r4/R4)  ]. 

For  the  inertia  1- subrange  form  (4-29),  the  integral 
(4-23)  can  be  performed  to  yield,  in  a  corresponding  range, 
the  function  L(k,u)  in  (4-24) .  We  find 

(4-33)  L(k,cu)  =  A2  v;le4/3k4k;l6/3  (L"1«  k+«  l’1)  , 


T 
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where 

(4-34)  k2  =  k2  +  (o)/ve)2,  ve  =  s^, 

A^  =  24/3(640/81)  7r2[r(2/3)/r(4/3)]A2  . 

Before  proceeding,  we  formally  generalize  to  the 
non-universal  range  where  the  condition  k+  »L*^  is  not 
satisfied.  We  continue  to  use  the  local-convection  and 
space-time  isotropic  approximations  as  at  (4-27).  There  is 
no  justification  for  this  extension,  but  it  may  prove  quali¬ 
tatively  useful  to  point  out  the  resulting  difference.  If 
the  reciprocal  size  of  the  energy -containing  eddies  is  ~ko(~L-^) , 
it  is  natural  to  suppose  that  the  alteration  in  L(k,u))  that 
would  result  from  using  a  ^(Rj)  roughly  appropriate  also 
where  k+  <  kQ  would  be  to  replace  k+  in  (2-54)  by  something 
like  (k^+4k^)^2,  where  the  factor  on  k2  is  uncertain  but  set 
equal  to  four  for  reasons  to  appear* 

To  examine  the  point  explicitly,  we  may  employ 
(4-31)  and  (4-32).  Form  (4-31)  corresponds  to  an  energy 

4 

spectrum  E(k)  that  behaves  appropriately,  i.e.,  as  k  ,  for 

k  «  kQ,  whereas  (4-28)  corresponds  to  one  that  varies  as 
-5 /3 

k  for  all  k.  Form  (4-31)  assumes  that  the  turbulence  is 
isotropic,  however,  even  in  the  energy-containing  range.  It 
does  not  reduce  to  form  (4-28)  for  R^«  L(or  k+  »  kQ) ,  as 
would  the  unmodified  von  Karman  interpolation  form  (Kinze  1959, 
(Eq. 3-136)),  but  suffices  for  the  present  purpose.  Inserting 
(4-32)  in  (4-23),  we  find  after  considerable  labor 

* - 

It  is  convenient  to  evaluate  (4-23)  by  use  of  four-dimensional 
spherical  coordinates. 
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k4 [40a4+45a2k2+8k4+5a (k2+a2) 1/2 (8a2  +5k2) ] 
(k£+a2) 5/2 [ (k2+a2) l/2+a ]5 

with  a=2kQ.  In  the  high-k  limit  where  u)/vck  +  0  and  a/k -*■(), 
(4-35)  becomes 

(4-36)  L(k,o>)  -  (32  »/9)v;1v*a2k'2  ; 

similarly,  in  the  limit  where  a/k+->C^  it  becomes 

(4-37)  L(k,o>)  -  (32  7r/9)v;1v4a2k4k;6, 

which  may  be  compared  with  the  unmodified  inertial -subrange 
result  (4-33)  that  pertains  to  the  same  limit.  A  simplified 
form  that  approaches  the  same  limits  (4-36)  and  (4-37)  and 
broadly  preserves  the  character  of  (4-35)  in  the  general  domain 
of  k  and  cu  is  given  by 

(4-38)  L  (k  ,oj)  =  (32=r/9)v;1v^a2k4(k2+a2)-3  . 

In  view  of  the  approximations  already  made  and  the  non¬ 
universality  in  the  range  in  question,  we  may  accept  (4-38)  in 
place  of  (4-35).  Since,  however,  (4-37)  and  (4-38)  justify 
the  initially  suggested  replacement  k“-*k  +a  for  general  a/k+, 
we  can  preserve  the  universal  character  for  kf»a  by  accepting 
not  (4-38)  but  the  generalized  form  of  (4-33): 

(4-39)  l.(k  ,<i>)  =A2v;1e4/3k4(k2+a2)  -8/3, 

where  2a~^  roughly  measures  the  size  of  the  energy-containing 
eddies. 

Inserting  (4-39)  in  (4-24)  ,  we  obtain  the  wavenumber- 


frequency  spectrum  of  pressure  in  the  mean  rest  frame. 


(4-40)  P(K,w)  =  CpP2€4/3v;1K^13/3, 

where 

-  K2+(oD/ve)2+a2, 

C2  -  (8/9)21/37r'3/^r(13/6)/r(4/3)]A2. 

A 

To  obtain  from  (4-40)  the  corresponding  spectrum  P(K,oo)  in 
the  measurement  frame  having  relative  velocity  -uq,  we  need 
only  use  (4-2)  to  find 

(4-41)  P(K,uO  -  C2p264/3v;1K^13/3,  where 

(4-42)  =  K2  +  (a:-uok1)2/v2+a2 

and  we  choose  uQ  =  (uqI0,0).  On  account  of  our  use  of  the 

A 

local -convection  approximation,  the  result  for  P(K,u>)  can  hold 
only  where  |u)-uQK^|/ve  <  K,  corresponding  to  the  region  of  the 

A 

convective  peak  (4-3).  Comparing  (4-40)  for  P(K,cu)  with  (4-4) 

(and  setting  a=0  for  the  inertial  subrange)  we  see  that  accord¬ 
ing  to  the  present  model,  the  function  0  there  introduced  is 
given  by 

(4-43)  0<x)  -  CpS^(l  +  s'2xV13/6  . 

In  view  of  the  previous  comparison  between  (2-96)  and  (2-97) 
for  E^(k,a))  ,  however,  we  should  expect  that  an  exact  calcula¬ 
tion  of  the  local  convection  approximation  for  P(K,cu)  by  means 
of  (4-26)  (still  on  the  basis  of  a  normal  isotropic  velocity 
distribution)  would  yield  a  function  0(x)  different  from  (4-43); 
in  any  case,  the  result  is  not  to  be  used  for  xv>l.  Eq  (4-43)  yields 
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for  the  value  of  the  coefficient  a^  defined  at  (4-6) 

('>-«!  a*  -  CpC13/6. 


where  c^y^  given  by  (4-9). 

/s 

In  the  present  model,  the  frequency  spectrum  P(cu) 
can  be  found  exactly  from  (4-7)  for  arbitrary  v0/u0;  from  (4-40) 
with  -►  I^,  we  obtain 

(4-45)  P(w)  -  (6.r/7)CpP2€4/3u;l(m2/u2  +  a2)'7/6, 

where  ufi  is  the  sum  of  the  mean  convection  speed  uQ  and  the 
effective  fluctuating  convection  speed  vg  in  quadrature,  given 
at  Eq.  (2-120)  (vg  -  s^vq)  .  *n  v  /u  -*0  (with  a=0)  , 

this  reduces  to  the  form  (4-8)  obtained  in  this  limit  more 
generally.  Eq.(4-45)  requires  the  nonviscous  condition  w£o/vo«l 
We  may  consider  Taylor's  hypothesis  for  pressure, 
as  in  Sec.  3.3  for  velocity,  by  regarding  the  ratio 

(4-46)  VUpP^pkp/P^), 

where  P  ^ ( k^)  is  the  ire  an- rest -frame  spectrum  in  streamwise 
wave  number  at  a  fixed  time  and  cross-stream  position  (x^) : 

P  (k.)  =  r  dia  r  dk,P(R,a), 

XX  —  OO  —  00  J 

and  u  is  an  effective  convection  velocity  open  to  choice. 

P 

From  (4-41)  and  (4-45)  we  find  Tp=l  exactly  provided  up=ue, 
given  at  (2-120). 

By  (4-22) ,  Eq.  (4-45)  yields  a  mean  squared  pressure 
<  P2>-  (8/9)21/3:r(2/3)/r(4/3)  ]A2P2e4/3a'4/3’ 


properly  independent  of  uo  and  in  accord  with  the  usual 


2  1/2  2 

order -of -magnitude  relation  <  p  >  ~p<  vq>. 

Concerning  the  properties  of  the  convective  ridge 

a  _  a^ 

in  P(K,u))  as  given  by  (4-40)  ,  at  fixed  k^  the  quantity  has 
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its  minimum  and  hence  P(R,uj)  its  maximum  at 

(4-47)  kL  =  klm=  (o)/uo)(uo/ue)2  =  tu/uc  ; 

this  minimum  value  of  is  given  by 

9  9  9 

*tm  "  Wue)  +k3+a  ; 

A 

the  half  width  Ak^  of  P(R,uj)  at  one-fourth  maximum  is  given 
for  Akx<  klm  by 

(4-48)  Akt  -  (ve/ue)(a)2/u^k^+a2)1/2 

(cf.  paragraph  following  (4-4)).  The  contours  of  constant 

P(R,cu)  in  K-space  are  defined  by  -  constant;  by  (4-42)  they 

are  ellipses  centered  at  k^  *  ±.kim»  ■  0  and  have  semiaxes 

along  k^  and  k^  in  the  ratio  vg/ue  , 

From  (4-41),  (4-47),  and  (4-48)  we  see  how,  as  uQ 

increases  from  zero,  the  ridge  in  P(R,u;)  sharpens  and  moves 

from  k.  ~  0  out  to  k,  =  uj/2v  at  u  =  v  and  then  again  toward 
1  1  e  o  e 

lower  k^;  similarly,  from  (4-45)  we  see  how,  as  uq  increases, 

A 

the  peak  in  P(oj)  at  co  =  0  becomes  broader  and  lower,  with 

A 

P(uj)  decreasing  for  oj<(  vr3/2)aug  and  increasing  for  co  v  (s'3/2)au 
We  consider  now  the  frequency  transform  of  the 
spatial  correlation  (cross-spectral  density)  in  the  measurement 
frame, 

(4-49)  ep(^)  =  /d2KP(K,^)eiR'  '  . 

Defining  the  normalized  magnitude  M  and  effective  convective 
velocity  u£  by  setting 

9p(’,w)  =  M(",uj)P(uJ)exp(iuri/uc)  , 


(4-50) 


from  (4-40)  and  (4-45)  we  find  that  u£  is  again  as  given  by 
Eq.  (2-120)  and 

(4-51)  M(c.co)  =  t21/6r(7/6)]“l(a2/U2+a2)7/12(:2+Nte2''2)7/12 

x  K7/6((o)2/u2+a2)1/2(C3+^2'2)1/2), 

where  ve«  ve/ue  as  at  (2-120^  The  elliptical  contours  of 
constant  M(",uj)  have  axes  along  and  '3  in  the  ratio  u£/ve 
inverse  to  that  pertaining  to  P(R,u))  ;  this  ratio  gives  also  the 
relative  scales  of  the  principal  coefficients  M(s|,0>o)), 

B(C3,u))2  M(0,C3,a>). 


Taking  the  liberty  of  comparing  with  experimental 
results  for  boundary-la  yer ,  as  opposed  to  homogeneous,  flow, 
with  the  rms  fluctuating  velocity  vq  in  vy  =  s ^vq  attributed 
a  value  characteristic  of  the  constant-stress  layer,  say 
vQ  =  0.1uq  In  the  regime  where  the  friction  velocity  v*  is 
given  by  v*/uq-  0.033,  we  obtain  for  lle/ve  a  value  in  good 
agreement  with  the  exper imenta  1  ly  observed  ratio  of  scales 
(-7.5)  (Corcos  1963,  Willmarth  and  Wooldridge  1962).  Other 
rough  correspondences  are  evident  between  present  results 
(e.g.  ,  (4-45) , (4-54) )  and  measurements  for  boundary-layer 
turbulence;  in  view  of  the  great  differences  between  these 
types  of  flow  these  correspondences  will  not  be  pursued, 
pending  a  careful  extension  of  the  present  approach. 

The  space-time  correlation  of  pxessure  in  the 
measurement  frame,  W  (”,~)  can  be  found  from 

(4-52)  Wp(-?,T)  =  fdLe~lJZ  ‘ 


or  be  replacing  '  ^  by  C^-uqt  in  the  rest-frame  correlation 


« 
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Wp(',t)  similarly  computed  from  9p(^,u>).  From  (4-50)  and 
(4-51)  we  find 

(A -5  3) 


Wp(C.T)  -  C^p2£4/3a-2/3t('1-u0T)2+^+v^2]1/3 


x  K2/3(«[(C1-u0t)2  +  C2«2t2]1/2), 

where  C2  =  2  3^3 (16/9) [ r (4 /3) ]~3A2 ,  in  which  the  argument  in 

2  2  2  2  2 

square  brackets  may  also  be  written  as  ue(T-"^/uc)  +^+\fe  . 

In  the  inertial -subrange  limit,  where  the  argument  of 

A 

becomes  small,  we  obtain  for  the  decorrelation  ^p(^, t) 

[-  Wp(0,0)-Wp(C,T)] 

//  c/  \  ?  x  o a2  2  4 / 3 r  /  r  ,2  .2^  2  2,2/3 

(4-54)  ^p( ^,-r)  =  2A  p  e  [(,fuoT)  +^3+veT  3 


As  apparent  also  from  (4-33)  and  (4-24)  ,  the  space-time 
isotropic  approximation  for  the  longitudinal  decorrelation 
^(r,T)  in  the  inertial  subrange  and  the  quas inormality  relation 
(4-25)  thus  lead  to  a  space-time  isotropic  form  (in  two  space 
variables)  also  for  the  (mean-rest-frame)  pressure  decorrela¬ 
tion.  Since  the  result  of  the  exact  local -convection  approxi¬ 
mation  to  ^(r,x)  was  well,  approximated  by  space-time  isotropy, 
as  seen  by  (2-80)  and  Fig.  1,  the  exact  approximation  to 

A 

tp(~, f)  within  the  validity  of  the  quasinormality  assumption, 
is  probably  well  approximated  by  (4-54).  In  the  approximation 
of  (4-27)  and  (4-28)  ,  by  (4-54)  we  have 

K'p(t.t)  =  2p2^(£  ,t)  ; 


this  result,  for  r=  0,  reduces  to  that  obtained  previously 
from  the  quasinormality  assumption  by  Batchelor  (1951)  and 
Obukhov  (1949) . 
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We  note  the  variously  defined  convection  velocities 
for  pressure.  That  defined  by  '^/t  where  dWp/cK^  0  is  exactly 
uQ.  According  to  approximation  (4-53),  that  defined  by  C ^/t 

A 

where  d'W^/b'i  =  0  is  uc,  given  at  (2-120),  the  same  as  the  one 
defined  at  (4-50)  by  the  cross-spectral  density.  The  convection 
velocity  defined  by  imposing  Taylor's  hypothesis  (4-46),  computed 
by  use  of  (4  41),  was  ug ,  given  by  (2-120).  In  short,  these 
convection  velocities  referring  to  pressure  are  just  the  same 
as  the  corresponding  ones  referring  to  the  streamwise  velocity 
component  in  the  same  approximation  of  space-time  isotropy  of 
the  longitudinal  velocity  correlation  (see  Table  2,  last 
colunn)  . 

We  turn  now  to  P(lT,u;)  in  the  low-wavenumber  region, 

K«au/v  ,  where  we  must  abandon  approximation  (4-27)  and  consider 
o 

dispersion.  For  this  purpose  we  employ  (4-26a).  The  integral 
(4-26a)  derives  contributions  from  domains  where  both,  one,  or 
neither  of  the  factors  E^  has  arguments  corresponding  to  the  iocal- 
convective  peak  in  this  function,  as  opposed  to  the  dispersive  tail; 
the  peaks  embrace  roughly  the  ranges  |it'|/vok'<  1  and 
Iu-oj'I/v  Ifc-lt'k  1,  as  seen  in  Sec.  3.2. 

We  consider  first  the  contribution  to  P(K,x),  say 
Pc(K,uO,  from  the  non-dispersive  domain  where  both  E^  factors  are 
determined  largely  by  local  convection  and  are  given  approximately 
by  Eq.  (3-9).  Examination  of  Eq.  (4-26a)  for  wL/vo»l  yields  in 
order  of  magnitude 

(4-55)  Pc(*,»)  ~c2»(€2/3vo'1K'8/3)2(o,/voK)'16/3 

=  pVA^/V13'3. 

P£  thus  has  a  form  roughly  independent  of  K  and  consistent  with  the 
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scaling  (4-4).  The  domain  of  integration  in  (4-26a)  that  contributes 

to  Pc  has  k2  w /vo~(»K) ,  so  that  local  convection  in  the  normal 

direction  can  generate  velocity  fluctuations  at  frequency  uu , 

though  the  local  convection  in  the  plane  of  K  cannot.  To  be 

justified  in  regarding  Pc  as  determined  by  the  inertial  subrange 

we  must  restrict  consideration  to  the  nonviscous  regime  where 

uu i  /v  «1,  in  accord  with  condition  (2-93a). 
o  o 

On  the  basis  of  (4-55)  we  obtain  for  the  contribution 

of  Pc  to  the  spectrum  (4-10)  of  average  pressure  on  a  moving  area 

of  radius  R  where  uuR  /ttv  »1  and  uuR  /nu  »1,  the  order  of  magnitude 
o  o  o  o  o 

estimate 


(4-56) 


q.»-o2£4/3v„10/V13/3r  -2. 


The  ratio  of  Qc  to  the  high -wavenumber  mean-convective  contribution 

A 

Q+  was  given  at  Eq.  (4-16).  It  is  noteworth  that,  on  account  of 
the  sufficiency  of  condition  (2-93a)  for  the  inertial  subrange  in 
the  local-convection  approximation,  it  was  not  necessary  to  assume 
KL  and  L/Rq  large  to  obtain  the  results  (4-55)  and  (4-56). 

We  now  attempt  to  estimate  in  order  of  magnitude  the 
contributions  to  P(K,uu)  from  the  domains  where  one  or  both  factors 
are  determined  largely  by  dispersion.  To  do  so  we  consider  crudely 
that  in  the  domain  of  the  local -convective  peak  E^(k,uu)  is  again 
given  by  (3-9)  and  outside  that  domain  (but  where  condition  (3-3) 
is  still  satisfied)  by  the  dispersive  form  (3-11).  Assuming 
uuL/vo»l,  we  join  the  domains  at  a  value  of  uu/vQk  somewhat  larger 
than  unity.  Examination  of  (4-26a)  for  uu/vqK»1  on  this  basis 
shows  that  the  results  depend  essentially  on  the  value  of  m  in 


(3-4)  and  (3-11). 
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In  particular,  for  m  •  3  both  the  convective-dispersive 
and  dispersive-dispersive  contributions,  defined  with  reference 
to  the  arguments  of  the  factors  in  (4-26a),  diverge  on  account 
of  the  dependence  of  E^(k *  ,uu*) ,  as  given  by  (3-11),  when  k'-K)  and 

It  is  inferred  that  the  actual  values  of  these  contributions, 
if  m  ■  3,  depend  on  E^k',^)  in  the  non-universal  range  where 
k'L^l,  uu'L/vo£l  and  the  inertial  subrange  form  (3-11)  no  longer 
applies  even  for  uu'/vok'»l.  In  the  range  where  conditions  (2-93a) 
(left  inequality)  does  not  apply,  we  estimate  the  order  of  magnitude 
of  E^(k,uj),  taking  L  as  the  wavenumber  scale  and  vq/L  as  the 
frequency  scale,  as 
(4-57)  E4(k,«»>-v0L2, 

i  2 

or,  if  kL«l,  as  (kL)  v  l  .  We  may  then  estimate  the  previously 
divergent  convective-dispersive  and  dispersive-dispersive  contri¬ 
butions  to  P(iC,uu)  with  the  integration  domain  properly  restricted 
to  the  universal  range  (2-93a),  and  estimate  the  contribution 
from  the  non-universal  domain  separately  by  use  of  (4-57).  We 
have  contented  ourselves  with  placing  upper  limits  on  the  orders 
of  these  contributions.  We  find  the  non-universal  contribution 
is  at  most  of  order  (u)L/vq)”^^  relative  to  the  universal  convec¬ 
tive  contribution  Pc  of  (4-55),  and  the  convective-dispersive  and 
dispersive-dispersive  contributions  are  likewise  of  higher  order 
relative  to  Pc  with  regard  both  to  (tt)L/vQ)”^  and  to  the  other 
small  expansion  variable  (uu/vqK)  Hence,  if  m  =  3  in  (3-4) 
and  the  quasinormality  form  (3-26a)  is  roughly  valid,  at  least, 
in  the  subject  domains  P(K,uu)  is  given  approximately  by  the 

*v 

convective  part  P  of  (4-55)  and  Q_  by  Q  of  (4-56),  as  stated  in 

c  ■*  c 


(4-12)  and  (4-14). 


4 


91 


Assumption  that  m  -  1  in  (3-4)  and  (3-11)  leads  to 
different  and  rather  curious  results  for  the  convective-dispersive 
and  dispersive-dispersive  contributions.  Denoting  these  contri¬ 
butions  to  P(K,id)  respectively  by  Pcd  and  Pj,  from  (4- 26a)  once 
more  for  u>L/vo»l  we  find  these  increase  with  decreasing  K  and 
estimate  their  orders  of  magnitude  relative  to  P  of  (4-55)  as 
given  by 

(4-58)  Pcd/Pc  -  („,/v0K)5/6(KL)‘5/6,  Pd/Pcd  -  (KL)'5/6. 

Thus  Pd  is  small  relative  to  P  ^  for  KL»1,  and  in  the  non- 
dispersive  limit  where  KL-*»  at  fixed  uu/vqK,  Pcd  in  turn  becomes 
small  relative  to  Pc-  In  a  domain  where  u>/voK^KL,  however,  Pcd 

A 

exceeds  P  .  As  for  the  corresponding  contributions  to  Q_ (uu )  of 
(4-13)  ,  where  a)Ro/vo»l»Ro/L  and  u>Ro/uo»l,  the  contribution  Qd 
from  Pd(K,uu)  is  controlled  by  non-universal  behavior  at  K£L~*  but 
is  estimated  as  ~-(Rq/L)  '  ,  and  hence  small,  relative  to  the  con¬ 
tribution  Qcd  from  Pccj(R,w).  Relative  to  of  (4-56)  we  then 

estimate  (^d/(^(u)Ro/vo)5/6(Ro/L)5^* ;  Qcd  varies  as  rather 

.  o 

than  Rq  ,  and  depends  on  vq  only  via  e.  Again,  where  Rq/L-*0  the 

contribution  Q  from  the  K-independent  P  thus  predominates,  but 
c  c 

where  uuRo/vo^L/Ro  the  largest  contribution  is  Qcd.  The  condition 
L/Ro»idRo/vo(»1)  for  applicability  of  (4-12)  and  (4-14)  thus 
obtained  for  m  ■  1  is  much  more  restrictive  than  the  condition 
uuRo/vo»l  obtained  for  m  -  3.  We  have  already  given  reason  to 
suppose  that  m  -  3  is  the  appropriate  value. 
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Assumption  that  m  -  1  in  (3-4)  and  (3-11)  leads  to 
different  and  rather  curious  results  for  the  convective-dispersive 
and  dispersive-dispersive  contributions.  Denoting  these  contri¬ 
butions  to  P(K,uu)  respectively  by  P£d  and  P^,  from  (4- 26a)  once 
more  for  uuL/vq»1  we  find  these  increase  with  decreasing  K  and 
estimate  their  orders  of  magnitude  relative  to  P£  of  (4-55)  as 
given  by 

(4-58)  Pcd/Pc  -  (o,/voK)5/6(KL)'5/6,  ?d/Pc<)  '  (KL)‘5/6. 

Thus  P^  is  small  relative  to  Pcd  for  KL»1,  and  in  the  non- 
dispersive  limit  where  KL-*»  at  fixed  uj/voK,  Pcd  in  turn  becomes 
small  relative  to  Pc<  In  a  domain  where  a/v^^KL,  however,  Pcd 

A 

exceeds  Pc<  As  for  the  corresponding  contributions  to  Q_(uu)  of 
(4-13)  ,  where  a)Ro/v0»l»Ro/L  and  u>Ro/uo»l,  the  contribution  Qd 
from  P(j(R,uu)  is  controlled  by  non-universal  behavior  at  K£L  *  but 
is  estimated  as  ~(Rq/L)  '  ,  and  hence  small,  relative  to  the  con¬ 
tribution  Qcd  from  PC(j(R,tu).  Relative  to  of  (4-56)  we  then 
estimate  Qcd/%r(u*0/''o)5 /S  (RofL)5 /(> ;  varies  as  R"1/3,  rather 

than  Rq  ,  and  depends  on  vq  only  via  e.  Again,  where  Rq/L-H)  the 
contribution  Qc  from  the  K-independent  Pc  thus  predominates,  but 
where  wR0/vo^L/R0  the  largest  contribution  is  <w  The  condition 
L/Ro»uuRq/vo(»1)  for  applicability  of  (4-12)  and  (4-14)  thus 
obtained  for  m  -  1  is  much  more  restrictive  than  the  condition 

ujR  /v  »1  obtained  for  m  ■  3.  We  have  already  given  reason  to 
o  o 

suppose  that  m  *  3  is  the  appropriate  value. 
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LOCAL -CONVECTION  APPROXIMATION  FOR  SHEAR  FLOW 


It  is  suggested  in  further  work  to  treat  space-time 
correlations  of  velocity  and  pressure  also  in  shear  flow  and, 
in  particular,  a  turbulent  boundary  layer  by  extension  of  the 
present  approach;  the  treatment,  however,  will  become  more  a 
guide  to  construction  of  a  heuristic  model  and  less  a  systematic 
application  of  kinematic  and  similarity  principles.  We  take 
one  step  here  by  suggesting  the  generalization  to  shear  flow 
of  the  separation  (2-6)  of  large-scale  convective  effects,  •  • 

and  of  the  local-convection  approximation  (2-59). 

We  assume  for  the  present  a  uniform  turbulent  shear 
flow,  taking  the  flow  in  the  x.  direction  and  velocity  gradient 

1  A 

u'  in  the  X2  direction;  relative  to  the  mean  rest  frame  for 
the  fluid  layer  at  X2  =0,  the  mean  velocity  is  thus 

(5-1)  u(x2)  =  u,x2I1, 

where  denotes  a  unit  vector.  We  may  denote  the  decorrelation 
tensor  of  the  fluctuating  velocity  in  this  frame  between  space- 
time  points  (x,t),  (x+r,t+x)  by  ^  (X2  ,x£  ,C>T)  »  with 
?  =  (';1,x^-x2,;3)  and  ?  *  (^,£3),  by  virtue  of  homogeneity 
in  the  x^-x3  plane.  Likewise,  the  decorrelation  tensor  measured 
in  a  frame  in  which  the  mean  velocities  at  x  and  x'  are  equal 
but  opposite,  i.e.,  the  frame  with  velocity  (1/2) '^(x^) +u  (x2)] 
relative  to  the  frame  of  definition  of  u^)  ,  by  the  homogeneity 
in  the  X2  direction  is  independent  also  of  X2  and  x^  except 
via  x^-X2  and  may  be  denoted  by  where  ?’  denotes  the 
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separation  vector  in  the  specified  frame.  By  the  usual  kine¬ 
matic  transformation  these  correlations  are  related  exactly  by 

(5-2)  (x2,x£,;,t)  =  <^j(r-(L/2Xu(x£)+u(x2)  ]t,t)  . 

As  in  the  instance  of  unsheared  flow,  we  may  expect  that 
provided  r  and  t  satisfy  certain  conditions,  namely  the  previous 
conditions  (2-5)  with  possibly  another  to  be  adjoined,  we 
may  regard  the  eddies  that  mainly  determine  ^°j(r,x)  as  being 
statistically  independent  of  the  large  eddies  that  contain 
most  of  the  turbulence  energy.  Then,  as  at  (2-6) ,  we  may 
approximately  separate  out  the  effect  of  the  large  eddies  by 
defining  an  intrinsic  decorrelation  tensor  ^j(r,x)  relative 
to  the  locally  co-moving  frame  which  is  independent  of  this 
motion : 

(5-3)  ^°j(?,t)  =  /d3vP(v)  (?-vt,x)  . 

is  still  affected  by  the  existence  of  shear. 

Now  we  distinguish  a  pseudo-convective  effect  of  the  shearing, 
just  as  we  did  of  spatial  fluctuating-velocity  dispersion  in 
the  unsheared  case  in  the  discussion  of  Eq.  (2-48)  and  in  the 
appendix,  Eqs.  (A -5)  ,  (A-7) .  Specifically,  we  imagine  the 
fluctuating  velocities  at  the  correlated  points  with  their 
differing  mean  velocities  as  tending  to  be  preserved  as  the 
medium  between  the  points  is  distorted  by  shear.  Thus,  apart 
from  residual  dispersive  effects,  we  expect  that  in  the  local 
proper  frame  in  question  the  velocity  correlation  for  the  two 
space-time  points  will  not  change  if,  for  any  change  6?'  in 
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the  separation  vector  r '  ,  a  corresponding  change  5-r  is  made 
in  the  time  separation  t  such  that  points  moving  with  the 
respective  mean  velocities  and  coinciding  with  x+r' ,  x  respec¬ 
tively  at  times  t-K ,  t  would  be  separated  by  ?'+6r'  at  the 
respective  times  t-K+6T  ,  t+Bt.  Hence,  if  we  define  a  decor- 
relation  tensor  ip^  with  this  pseudo-convective  effect  removed 
by  writing 

(5-4)  i'ij  (f  ’  ,t)  =  ^j(r'-[u(xJ)-u(x2)]T,T>  , 

then  ^j(r",T)  will  depend  on  t  (at  fixed  r")  only  on  account 
of  the  residual  dispersive  effect.  Eqs.  (5-3)  and  (5-4)  yield 

(5-5)  r,T)  =  /d3vP(v)^j(r-[v+u(x^)-u(x2)]x,T) 

and  (5-2)  then  yields  in  terms  of 

We  define  a  non-dispersive  approximation,  generalizing 
the  local-convection  approximation  for  unsheared  flow,  by 

(5-6)  iij(r",T)  =  ^(*",0)  . 

Since,  by  (5-5)  without  approximation,  ^°j(r,0)  =  ^j(r,0), 

Eqs.  (5-5)  and  (5-6)  yield 

(5-7)  J®j(r,T)  -  /d3  vP(v)  ^  (r-[v+u(xp -u(x2)  ]r  ,0) 

as  the  expression  of  the  space -time  correlation  in  terms  of 
the  space  correlation  in  the  non-dispersive  approximation. 

In  the  frame  referred  to  by  (5-2)  this  similarly  becomes 

(5-8)  ^/i j  (x2 , x2  , C  »t )  - 

/d3^P(v)\ti:j  (x2  ,xJ,C-  [  (3/2)u(x2)  - (1/2) u(x2) +v]t ,0) . 
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If  the  shear  is  truly  uniform,  as  assumed  at  (5-1) ,  the 
square -bracketed  effective  local  convection  velocity  in  (5-8) 
may  be  written  also  as  (1/Z)[u(x2)+u(x2)  1+u'  ^i(x2“x2)  The 

consequences  of  this  non-dispersive  approximation  can  be 
extended  to  the  treatment  of  pressure  also  in  the  present 
anisotropic  situation  by  assumption  again  of  a  quasinormal 
velocity  distribution.  (The  mean  shear-turbulence  interaction 
contribution  naturally  must  now  be  included.) 

As  for  the  conditions  on  ?  and  t  defining  the 
universal  range  where  (5-3)  may  be  applied,  we  apparently 
require,  apart  from  the  conditions  (2-5),  that  [u^)  -u^)  ]t«L; 
since,  however,  in  shear  flows  typically  vo~u'L,  the  added 
condition  is  normally  already  implied  by  conditions  (2-5) . 

As  written,  Eq.  (5-7)  or  (5-8)  is  meaningful  also 
when  applied  to  a  non-uniform  shear  flow,  e.g.  ,  a  boundary 
layer.  If  it  is  so  applied,  the  crudity  of  the  approximation, 
at  best,  must  be  greater  or  the  domain  of  applicability  smaller. 
In  any  case,  it  remains  true,  as  we  saw  in  the  corresponding 
instance  in  unsheared  flow,  that  the  non-dispersive  approxi¬ 
mation  does  not  suffice  to  give  a  good  approximation  to  the 
tail  of  the  wavenumber -frequency  spectrum  where  oj/v^k  is 
large. 

Other  well  known  complications  are  involved  in  treat¬ 
ing  pressure  fluctuations  in  a  boundary  layer,  as  opposed  to 
the  unsheared  homogeneous  flow  considered  in  the  preceding 
sections,  even  with  regard  to  the  domain  of  the  convective  peak 
in  the  wavenumber  spectrum  to  which  we  might  tentatively  apply 


(5-8).  These  effects  are  attended  by  the  entrance  of  the  length 
parameter  v/v*  characterizing  the  viscous -sublayer  thickness 
and,  with  reference  to  spatially  dependent  quantities,  also 
the  distance  ^  from  the  wall  . 


Among  other  consequences  of  the  wall,  the  spectrum  P(R,u>)  of 
pressure  on  the  wall  in  the  approximation  of  incompressibility, 
streamwise  homogeneity,  and  vanishing  viscosity,  varies  as  K* 
when  K«6"l,  where  6  is  the  boundary- layer  thickness,  and 
hence  vanishes  a**  K+0,  in  contrast  to  the  nonvanishing  limit 
P(0,cu)  in  the  unsheared  homogeneous  case  (Kraichnan  1956, 
Phillips  1956). 
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6. 


SUMMARY  AND  CONCLUSION 


Kolmogorov's  principles,  as  expressed  by  Eq.  (2-6), 
provide  a  basis  for  the  treatment  of  Eulerian  space-time  corre¬ 
lations  in  the  universal  range  (2-5)  in  homogeneous  turbulence 
by  permitting  the  explicit  separation  of  the  kinematic  effect 
of  convection  by  the  large  eddies.  From  the  related  presumption 
that  the  space-time  correlation  in  the  local  co-moving  frame 
has  the  commonly  accepted  similarity  character  independent  of 
the  velocity  distribution  P(v) ,  and  from  the  order  of  magnitude 
of  the  coefficients  involved,  it  was  indicated  that  the  local- 
convection  (non-dispersive)  approximation  (2-59)  has  validity 
in  the  inertial  and  viscous  subranges. 

In  the  local -convect ion  approximation,  with  an  iso¬ 
tropic,  normal  velocity  distribution,  the  space-time  structure 
functions  of  velocity  were  computed  explicitly  for  the  inertial 
subrange  [Eq.  (2-80)  and  Fig.  1]  and  for  the  viscous  subrange 
[Eq.  (2-86)].  These  functions  were  found  to  be  nearly  space- 
time  isotropic  in  the  inertial  subrange  and  exactly  so  in  the 
viscous,  with  differing  velocity  scales  for  longitudinal  and 
transverse  components.  The  result  in  the  former  case  was 
coopared  with  the  space-time  correlation  in  grid  turbulence 
measured  by  Favre  and  associates.  The  computed  envelope  of 
the  time  correlations  at  fixed  spatial  separations,  for  the 
presumed  value  of  turbulence  intensity  in  the  experiment, 
decreases  somewhat  too  slowly  to  agree  well  with  the  measured 


one.  This  discrepancy  is  presumably  due  mainly  to  neglect  of 
dispersion,  which  is  significant  at  the  Reynolds  number  of  the 
experiment.  At  higher  Reynolds  number,  such  fitting  of  the 
space-time  correlation  may  usefully  supplement  measurements  of 
pure  spatial  correlations  and  wavenumber  spectra  to  yield  the 
Kolmogorov  coefficient  [Eq.  (2-56)]. 

For  the  unsheared  homogeneous  turbulence  assumed  in 

this  work,  Taylor's  hypothesis  relating  wavenumber  spectra 

in  the  mean  rest  frame  to  frequency  spectra  in  the  measurement 

freme  (velocity  -uq)  was  examined  for  the  inertial  subrange. 

In  the  local-convection  approximation  the  hypothesis  was  found 

to  be  exact  for  arbitrary  turbulence  intensity  VQ/U0*  provided 

an  appropriate  effective  convection  velocity  is  used  to  relate 

wavenumber  to  frequency  for  the  respective  spectra  [Eq  (2-116)], 

Ft'  small  v  /u  ,  this  convection  velocity  differs  from  u  by  a 
o  o  J  o  J 

2  2 

term  of  relative  order  vq/u^  depending  also  on  the  quantity 
whose  spectra  are  in  question,  and  other  commonly  defined 
convection  velocities  also  difrer  by  such  terms  [Table  2]. 

The  related  cross-spectral  density  of  streamwise  fluctuating 
velocity  in  the  measurement  frame  was  also  computed  in  the 
space-time  isotropic  approximat  ion  .  Eq .  (2-119)]. 

The  local -convect ion  approximation  for  fixed  finite 
scale  L  of  the  energy-containing  eddies  does  not  yield  the  wave¬ 
number-frequency  spectrum  of  turbulence  correctly  up  to  arbitrarily 
large  values  of  iL'/v^k;  rather,  velocity  dispersion  in  space 
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and  time  due  to  local  eddy  distortion  enters,  and  it  is  necessary 
to  consider  the  frequency  spread  of  the  intrinsic  (proper -frame) 
spectrum  E^(k ,w) .  The  standard  similarity  dependence  of  E^  for 
the  inertial  subrange  was  accepted  (Eq.  (3-1)],  and  an  argument 
given  for  the  dependence  E^(k,«)ak4  as  k  e/w  -K).  At  small  vok/uj, 
the  intrinsic  and  mean-rest- frame  spectra,  E^  and  E^,  could  then 
be  inferred  uniquely  within  a  constant  factor  [Eqs.  (3-11),  (3-15) 
with  m  -  3] .  An  explicit  form  having  plausible  properties  was 
suggested  to  be  used  for  E^(k,uu)  [Eq.  (3-4)].  Its  most  conse¬ 
quential  property  is  to  have  a  finite  second  moment  with  respect 
to  frequency  [Eq.  (A-18),  n  =  2] .  The  resulting  dispersive 
correction  to  the  rest-frame  spectrum  as  given  by  the  local- 
convection  approximation  was  computed  for  moderate  uL)/vQk  [Eq.  (3-10)]. 
At  given  k,  dispersion  increases  the  energy  E^(k,(u)  in  the  higher 
range  of  u)/vQk  at  the  expense  of  that  Jn  the  lower. 

From  the  assumed  intrinsic  spectrum,  the  dispersive 

departure  from  validity  of  Taylor's  hypothesis  was  found.  The 

departure  from  unity  of  the  pertinent  ratio  of  spectra  at 

2  -2/3 

wavenumber  k^  is  of  the  order  of  (v0/u0)  (k^L)  and  hence, 

for  k^L»  1,  small  relative  to  the  non-dispersive  departure 
compensable  by  the  above  mentioned  redefinition  of  convection 
velocity  [Eqs.  (3-21),  (3-24)].  The  dispersive  correction  to 
the  local -convection  approximation  to  the  space-time  velocity 
correlation  was  also  found  in  two  limits;  it  is  of  the  relative 
order  of  (vox/L)2/3if  (vor/L)1/2«l  and  (vqt) 3/2/rL1 /2»1 ,  and  the 
order  of  (vQT)2/r4/3L2/3  if  vQT/r«l  [Eqs.  (3-28)  -  (3-31)].  The 
effect  on  the  comparison  with  measured  results  for  grid  turbulence 
was  computed.  The  magnitude  of  this  effect  depends  on  the  ratio 
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B/Aq  (Eqs.  (2-44),  (2-45)]  measuring  the  relative  decorrelating 
influence  of  time  delay  and  of  spatial  separation  in  the  proper 
frame;  for  a  likely  magnitude  of  B/Aq  and  the  Reynolds  number  of 
the  experiment  the  computed  dispersive  effect  is  substantial. 

Inferences  from  kinematic  and  similarity  arguments  for 
an  inertial  subrange  can  be  drawn  also  concerning  pressure  spectra 

A 

in  the  measurement  frame.  The  wavenumber- frequency  spectrum  P(K,a>) 
referring  to  a  plane  parallel  to  uq  was  indicated  to  have  the  simi¬ 
larity  form  (4-5)  in  the  neighborhood  of  its  convective  ridge 
(|  uu-uo*K|  <vqK).  For  vo/uq«1  ,wL/uo»1  ,  and  ^^0/v0«l,  the  point 

A 

frequency  spectrum  P(ju)  could  be  inferred  from  (4-5)  to  have  the 

-7/3 

dependence  (4-8),  varying  as  i  .  The  spectrum  Q(uu)  of  average 
pressure  for  a  circular  area  of  radius  Ro,  though  not  directly 
measurable,  was  considered  with  a  view  to  future  treatment  of 
a  turbulent  boundary  layer.  For  a  small  area  (jlRo/uo<x<1)  ,  the 
area  correction  is  given  by  (4-11).  For  a  large  area  (ojR0/uo^2tt)  , 
contributions  associated  with  the  convective  peak  of  P(K,oj)  and 
with  the  low-wavenumber  region  (Kv2~Rq^)  are  distinguished,  the 
former  varying  as  Rq  .  From  the  form  of  P(K,jl)  in  the  convective 
region  the  former  contribution  is  given  by  (4-15). 

With  assumption  of  quasinormality  of  the  velocity 
distribution  and  the  non-dispersive  approximation  of  space-time 
isotropy  for  the  longitudinal  fluctuating  velocity  correlation, 

A 

the  spectrum  P(K,cj)  applicable  in  the  neighborhood  of  the  con¬ 
vective  ridge  was  computed  explicitly  for  the  inertial  subrange 
[Eq.  (4-40),  a  =  0]  (and  formally  extended  into  the  non-universal 
range,  assumed  isotropic) .  In  this  approximation  the  point 

/s 

spectrum  P(oj)  was  computed  explicitly  for  arbitrary  vq/u^ 
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[Eq.  (4-45)],  and  the  cross-spectral  density  and  space-time  corre¬ 
lation  of  pressure  were  also  obtained  [Eqs.  (4-50),  (4-53),  (4-54)]. 
The  variously  defined  convection  velocities  are  the  same  as  those 
for  the  streamwise  velocity  component  in  the  same  space-time 
isotropic  approximation.  From  the  assumption  of  quasinormality , 

A 

P(K,id)  was  estimated  also  in  the  low-wavenumber  part  of  the 
inertial  domain  where  u)/vqK»1,  L"1«tu/vo«£^1 ,  and  ai/uoK»l. 
Throughout  this  domain  (assuming  m  -  3  in  (3-11)),  P(K,uj)  is  given 
approximately  by  the  local-convective  wavenumber- independent  form 
(4-55).  Correspondingly,  the  low-wavenumber  contribution  to  the 

/s 

moving-area  spectrum  Q(id)  for  ujRq/vo»1  and  uuRo/uo»l  assumes 

_  o 

the  form  (4-56),  varying  as  Rq  .  The  ratio  of  the  low-wavenumber 

/\ 

to  the  mean-convective  contribution  to  Q(u)  is  given  by  (4-16). 

The  kinematic  separation  of  large-scale  convective 
effects  basic  to  ths  work,  and  the  related  local-convection 
approximation  for  space-time  correlations,  were  plausibly 
extended  to  the  more  general  case  of  shear  turbulent  flow 
[Eqs.  (5-5),  (5-8)]. 

The  present  work  provides  a  basis  for  further 
comparisons  with  measurements  on  grid-produced  and  other 
relatively  unsheared  turbulent  flows  and  a  springboard  for  an 
attack  on  boundary- layer  turbulence  with  reference  to  proper¬ 
ties  not  yet  adequately  explored. 
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APPENDIX 


PROPERTIES  OF  THE  INTRINSIC  SPACE-TIME 
CORRELATION  OF  VELOCITY  IN  THE  INERTIAL  SUBRANGE 

If  we  propose  to  calculate  explicit  corrections  to 
the  local -convection  approximation  (2-54) [or  (2-59)]  in  the 
inertial  subrange,  i.e. ,  to  consider  terms  of  next  lowest  order 
in  R  /L  (with  R  given  by  (2-52)  and  L  the  scale  of  the  energy- 
containing  eddies) ,  we  should  examine  further  the  degree  of 
validity  of  the  basic  separation  (2-6)  of  large-scale  convec¬ 
tion.  We  therefore  consider  the  possibility  that  the  intrinsic 
decorrelation  ^(r,i)  for  given  v  depends  significantly  on  v  and 
likewise  on  the  scale  L,  the  rms  velocity  fluctuation  vq ,  and 
any  parameters  of  large-scale  anisotropy.  We  denote  such 
possible  dependence  in  a  partially  suppressed  fashion  by  writing 

ip  =  L~^,r,T).  Later  we  must  specify  more  precisely  the 

v 

definition  of  v.  Eq.  (2-6)  is  thus  generalized  to 

(A-l)  <//(  r,T)  =  /  d3  vP(v)  ^_(L_1,r-VT  ,t)  . 

v 

Now  let  ip( r,T)  again  denote  a  suitable  function  of 

the  form  ^/(?,t)  =  (er)2/f,3F(ei^3r"2^3tTl)  [Eq.  (2-42)]  independent 

of  v  and  the  large  eddies;  ^(r,T)  is  to  be  so  chosen,  if  possible, 

that  when  inserted  in  place  of  ^_( r,-r)  in  (A-l)  ,  the  computed 

v 

correction  A^(r ,t)  =  ^(r ,t) -^c(r,T)  to  the  local -convection 
approximation 

(A -2)  ^c(r, t)  =  /  d3vP(v)^(r-vT  ,0) 

is  correct  to  lowest  order  in  Rc/L.  For  any  tf/(r,x)  we  can 


rewrite  (A-l)  as  an  (exact)  equation  for  A  ip: 

(A-3)  A^(r,x)  =  /  d^  vP(v)  [^/_(L”\r-vx  ,0)-^(r-vx,0)  ] 

v 

+  /  d3vP(v)  [^_(L"I,r-VT,x)-^_(L"1,r-VT,0)  ]  . 
v  v 

The  first  integral  represents  an  increment  due  to  non-vanishing 
L  ^  at  zero  dispersive  time  delay,  and  the  second  represents 
an  increment  due  to  the  dispersive  time  delay  at  fixed  v 
and  L.  The  objective  is  then  to  define  v  and  choose  ip( ?,t) 
such  that  A  ip  in  (A-3),  to  lowest  order  in  Rc/L,  is  equal  to 

(A-4)  /  d  J  vP(v)  [^(r-v-r  ,t)  -tf/(r-vx  ,0)  ] . 

With  regard  to  the  first  integral  in  (A-3)  ,  by  (A-l) 

O  _  _  ~  _1  _  _  _ 

we  have  /  dJvP(v)^-(L  ,r-vx,0)  =  ip( r,0);  furthermore,  to 
the  lowest  two  orders  in  r/L,  tf-(r.O)  has  the  form 

'/ (r  ,0)->F(0)  (er)  ^3-ALv^  (r/L)  ^  ,  (A^  =  constant) 

as  exemplified  by  the  result  for  the  von  Karman  interpolation 
form  (Hinze  1959,  Eq.  (3-136)).  Hence,  the  contribution  to 
A <j,  from  the  first  integral  in  (A-3),  relative  to  the  magnitude 
of  ^c,  is  of  the  order  -v^(Rc/L)^-:  (eRc)^^~-(Rc/L)^^.  On  the 
other  hand,  the  value  of  (A-4)  on  the  basis  of  a  ;/(r,x)  charac¬ 
terized  by  (3-5)  is  estimated  by  reference  to  Eqs.  (3-28) , 

_i  2/3 

(3-29),  to  lowest  order  in  L  ,  as  of  relative  order  (Rc/L) 
and  hence  of  lower  order  than  the  first  term  in  (A-3) .  It 
remains  then  to  consider  whether  the  second  term  in  (A-3) 
alone  agrees  with  (A-4)  in  this  lowest  order. 
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We  are  free  to  choose  v  in  (A-l)(and  hence  (A-3) 
and  (A-4))  to  include  contributions  to  fluctuating  velocity 
from  all  eddies  down  to  some  minimum  size  ~L*  ,  say  (cf.  Inoue 
1951);  for  example,  we  might  set  L'  2  Rc(L/Rc)5  for  some 
0  <  6  <  1  before  letting  Rc/L  ■*  0.  The  smaller  L'  is,  the  less 

these  eddies  are  removed  from  the  size  ~RC  of  eddies  mainly 

determining  (L_1  ,r-vx ,x)  (with  |r-vx|~Rc),  and  hence  the 
v 

more  likely  is  to  depend  significantly  on  v  on  account  of 
v 

distortion  of  these  eddies  by  eddies  of  size  >  L* .  On  the 

other  hand,  unless  L'~R  (i.e.  6  =  0),  then  ^_(L~^ ,r '  ,x)  will 

c  v 

have  to  include,  by  means  of  its  dependence  on  the  dispersive 

time  delay  t,  the  decorrelating  effect  during  x  due  to  the 

residual  convection  velocity  arising  from  eddies  in  the  size 

interval  from  ~L'  down  to  ~R  .  Consideration  on  this  basis 

c 

indicates,  indeed,  that  the  second  term  in  (A-3)  is  of  lower 

order  in  R£/L  than  (A-4) ,  and  hence  cannot  reduce  to  the  latter 

unless  L'~RC.  If  L'~Rc,  the  desired  reduction  can 

occur  if  the  intrinsic  decorrelation  ^_(L  ^,r-vx,x)  (with 

v 

ir-vx|~Rc)  is  statistically  independent  of  the  motion  due  to 
eddies  of  size  ~R  (L/R  )^  even  for  5-K).  We  leave  further 
clarification  to  explicit  dynamic  treatments  (e.g.  ,  Kraichnan 
1959,  1966)  and,  as  stated  in  the  text,  accept  the  similarity 
form  (2-42)  for  use  in  (A-l)  where  a  definite  form  is  required. 

As  a  guide  to  the  possible  properties  of  the  intrinsic 
decorrelation  $  (r ,  x)  [  =  ^(r  ,x)  ]  ,  assumed  of  the  form  (2-42),  we 
append  discussion  of  certain  conjectures.  Previously  we  dis¬ 
tinguished  two  sources  of  decorrelation  in  ^(r,x)  relative  to 
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^(r ,0)  [=^/(r ,0)  ] :  (1)  differential  convection  due  to  velocity 
dispersion  over  the  distance  r  acting  during  time  t; 

(2)  pseudo -La grangian  velocity  wander  during  x ,  existing  even 
when  r  =  0  With  regard  to  the  latter,  since  ^(0,t)~€|tI,  this 
effect  may  be  considered  a  random  walk  of  the  velocity  at  a 
given  point  in  the  proper  frame,  corresponding  in  the  related 

•fa 

Lagrangian  case  to  white  noise  in  acceleration  . 

For  identifiable  effects  contributing  to  velocity 
decorrelation  with  time  in  a  given  Eulerian  frame,  one  approxi¬ 
mate  mode  of  inclusion  to  be  considered  for  possible  validity 
is  this:  the  effect  is  regarded  as  producing  during  x  a 
displacement  vector,  with  some  independent  probability  density, 
and  for  each  realization  of  this  displacement  the  velocity 
correlation  is  supposed  to  be  formed  as  the  Eulerian  correlation 
function  that  applies  in  the  absence  of  the  effect  with  a  dis¬ 
placement  argument  that  is  the  vector  sum  of  the  previously 
applicable  one  with  the  new  one.  This  mode  is  exemplified  by 
the  basic  Eq.  (2-6)  that  treats  the  effect  of  large-eddy  motion 
on  time  correlation  in  the  mean  rest  frame.  A  second  simple 
mode  of  inclusion  of  a  given  decorrelating  effect  may  be 
considered:  the  effect  is  regarded  as  producing  a  statistically 
independent  contribution,  not  to  the  displacement  vector  in  the 
correlation  function,  but  to  the  velocity  decorrelation  itself. 
The  former  approach  is  surely  indicated  (through  statistical 
independence  may  be  doubtful)  where  the  effect  in  question  is 

though  it  is  useful  to  distinguish  these  effects  we  cannot 
expect  that  the  actual  dynamics  permit  a  clean-cut  separation 
of  the  two. 
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due  to  an  identifiable  convection  of  a  presumed  relatively 
static  eddy  configuration,  either  a  convection  common  to 
both  spatial  points  in  the  frame  of  the  Eulerian  correlation 
in  question  or  (in  cruder  approximation)  one  that  arises  from 
a  relative  velocity  between  the  two  points  . 

Of  the  two  effects  distinguished  in  the  first  para¬ 
graph,  the  first,  due  to  velocity  dispersion  over  r,  should 
thus  be  treated  by  the  convection  view  (recall  also  the  dis¬ 
cussion  of  Eq.  (5-4)).  Which,  if  either,  treatment  is 
appropriate  for  the  second  effect,  velocity  wander,  is  less 
clear.  The  view  that  velocity  wander  contributes  to  an  effec¬ 
tive  displacement  vector  is  akin  to  the  independence  hypothesis 
of  Corrsinl96l,  Saffman  1963,  referring  to  the  relation  of 
Lagrangian  and  Eulerian  correlation  functions.  If  we  support 
this  view,  we  incline  to  a  rough  conjectured  form 

(A-5)  ^(r  ,x)  -Aqc2/3  /d3xx  /  d3  x2P(xL  ,x2)  jr-c^  |2/3, 

where  v^  =  (e  |t  |)  '  and  vr  was  defined  after  (2-43).  P(x^,x2) 

is  the  joint  probability  of  contributions  to  the  effective 
displacement  vector,  for  points  separated  by  ?,  of  c^v^x^  due 
to  differential  convection  and  C2VXT*2  ^ue  t0  vel°city  wander  ; 
the  probability  variables  x^,x2  are  assumed  normalized,  so 
that  the  rms  values  of  the  two  displacement  contributions  are 
clvr I  t  I  and  c2vt|t|,  and  c^,c2  are  of  the  order  of  unity 
since  vf  is  of  the  order  of  the  differential  velocity  over  r 

x - ' - 

We  note  that  the  assumption  of  Chandrasekhar  (1956)  corres¬ 
ponded  to  the  opposite  view  that  even  the  decorrelation  with  time 
due  to  local  convection  by  large  eddies  appears  in  the  form 
of  an  independent  additive  component  in  the  velocity  decorrela¬ 
tion  [see  Eq.  (2-50)]. 


1 


* 


and  of  the  order  of  the  velocity  wander  during  t.  In  the 
spirit  of  the  conjectured  approximation,  we  regard  P(x^,x2)  in 
(A-5)  as  independent  of  r,  t  and  of  the  directions  of  x^,  x2, 
except  that  it  may  depend  significantly  on  where 
Pi=  x^*r/lx^lr,  i.e.  ,  the  probability  of  a  relative  velocity 
between  points  (x,t)  and  (x+r.t-Hr)  may  well  depend  on  the 
direction  with  respect  to  the  separation  vector  r.  If  in 
addition,  the  processes  associated  with  x^  and  x2  may  be 
considered  independent,  we  may  set 

(A-6)  PCxpX^  *  P1(x1)P2(x2)  =  (4^“2p1(x1,p1)p2(x2). 

If  we  support  rather  the  view  that  velocity  wander 
contributes  directly  and  independently  to  the  velocity  decor¬ 
relation,  we  write  in  place  of  (A-5) 

(A-7)  J(r,T)  *  Aq€2/3  /d3i1P1(i1)|f-c1vrTi1|2/3+  (£r)2/3G(z), 

where  z  was  defined  by  (2-43).  The  second  term  represents 
velocity  wander;  by  (2-45)  G(z)«z  as  z-+»,  but  G(z)  may  vanish 
rapidly  enough  as  z-*0  so  that  the  first  term  alone  gives 
f (r ,t) -^(r,0)  to  lowest  order  in  z.  The  simplest  assumption 
for  G(z) ,  however,  would  be  simply  G(z)  *  Bz  whence 

(A -8)  (cr)2/3G(z)  =  7 (0 ,t)  =  Be  I t l  . 

This  form,  being  independent  of  r,  yields  a  contribution  to 
E(k.r)  or  E^(k,u))  having  singular  dependence  as  6  (k)  and  hence 
could  be  useful,  if  at  all,  only  for  quantities  entailing 
integration  over  wave  number. 
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If  we  could  suppose  that  P^(x^)  in  (A-7)  is  isotropic 
(i.e. ,  independent  of  ^j)  and  also  normal,  we  obtain,  as  at 
Eq.  (2-80)  (i-0) ,  a  contribution  to  (Kr,x)  from  the  first  term 
in  (A-7)  given  by 

A0(€R)2/3tl-H0(clZ)] 

where 

2  2  2  2  2  2 
R  =  rSs;(vrx)z  =  r  (l+s^zZ)  . 

An  analysis  of  the  corresponding  contribution  to  the  spectral 
density  E^(k,oo)  appears  to  yield  a  negative  value  (varying  as 
c5/2k2af 13/2)  in  the  limit  0  *  cu/t1^2/3-**,  however,  an 
unacceptable  result . 

Similarly,  examination  of  the  E^(k,cu)  given  by  (A -5) 
in  the  simplified  case  of  c^=0with  isotropic  P2OC2)  yields 
E^(k,(u)«  -  E(k)€k2oT^(«-  as  Q+®>,  i.e.,  a  negative 

spectrum  again  in  this  limit.  It  is  natural  to  consider  also 
the  variant  of  (A-5)  ,  in  the  simplified  case  where  c^=0,  obtained 
by  adding  displacements  r  and  c^vrx  in  quadrature: 

(A -9)  i(r,T)=  A062/3(r2+c^€lTl3)l/3. 

This  form  also  yields 

(A- 10)  E4(k,oj)« -€  5/3k1/3uf4  as 

in  fact,  by  virtue  of  the  discontinuous  third  derivative  with 
respect  to  t,  any  form  leading  to 

E(k,T) -►  A(k) -B(k)  |t  I  3  as  t^O  , 

where  the  function  B(k) >  0  for  some  k,  yields  this  same  negative 
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result  for  such  k  in  the  limit  £}-*■». 


A  similar  variant  on  the  first  term  of  (A-7)  (or  on 
(A-5)  with  C2  ■  0) ,  namely 

(*•11)  i(«.0-A0«2'V ^62/3r2/3,.y/3  , 

on  the  other  hand,  yields  a  positive  result 
(A-12)  E4(k,oj)  «  €5/2k2oT13/2  as  ft  ■*  °o  . 

The  quadrature  variant  of  the  complete  form  (A-5)  is 

c*-i3)  ;(r,o-v2/V+c2l£2/v/v*c2£MV/3 . 

3 

It  appears  that,  despite  the  c^  term, the  term  «  [t|  will 
cause  this  form  to  yield  once  more  a  negative  E,.  (k ,gj)  in  the 
limit  of  large  ft,  so  that  this  form  also  fails  in  this  regime. 
Though  (A-13)  is  thus  unacceptable  even  as  a  tentative  approxi¬ 
mation,  it  may  correctly  indicate  the  dependence  of  ip(r,i)  in 
the  limits  of  large  and  small  z. 

In  the  limit  z  ■*  0,  if  /  d  xiP(x^  >*2) =  0,  farm 
(A-5)  yields  to  lowest  order  in  z 

(A-14)  1  (r  ,t)  (rjOh  1  i(r,0)z2  «  e4/3r-2/3x2  . 

In  particular,  if  (A -6)  holds  and 


(A- 15) 


00  3  1. 

dxx  /  dp.  M-P ^  (x  ,p)  =  0  , 

O  “1 


More  generally,  an  assumed  form 

;(r.T)  -  A0«2'3(rl'V2«2'V'— */3tV'3 

2/3 

(which  has  the  required  form  (er)  F(z))  also  yields  the 
result  (A-12K  If  m  =  1/2,  this  form  reduces  to  (A-ll);  if 
m  *  3/4,  it  has  instead  the  required  behavior,  (*  |t|)  as  z  •*  <x , 


no 


then  (A -14)  is  the  limiting  form.  But  Eq.  (A-15)  must  be 
valid,  in  fact,  since  otherwise  the  separation  of  all  pairs 
of  fluid  elements,  on  the  average,  would  be  increasing, 
contrary  to  the  condition  of  compressibility.  Likewise  form 
(A-7),  if  /  d  =  0  and  G(z)  -»  0  faster  than  z  as 

z  -»  0 ,  also  yields  this  result;  so  also  does  form  (A-13)(cf.  Eq. 
(2-48)).  On  the  other  hand,  if  G(z)«  z  as  z  ->  0,  which  would  be 
so  if  form  (A-8)  is  valid  in  this  limit^  form  (A-7)  yields  instead 

(A-16)  \p(r ,t)  -^(r,0)-^i/(0,T)«  c|t|; 

in  such  case,  velocity  wander  rather  than  differential 
convection  represents  the  leading  contribution. 

In  the  opposite  limit  where  z  •*  <»,  if  (A-6)  holds, 
form  (A-5)  yields  to  lowest  order 

tf/(r,T)-(H0,x)«  ((0,t)z-1  «  (er)2^3; 


form  (A-13)  also  yields  this  result.  On  the  other  hand,  form 
(A-7)  in  this  limit  yields 


*'(r,T)-*(0.-.)-(cr)2/3z2/3 


£8/9r2/9|T,2/3_ 


We  can  avoid  constructing  functions  r,x)  that  lead 
to  negative  E^(k,uj)  by  starting  from  the  latter.  In  the  inertial 
subrange  (3-3)  the  requisite  form  for  t^(k,tu)  is  given  by  (3-1) 
in  the  text.  Form  (3-1)  yields 

(A-17)  i(r,0)  -  (9/10)r(l/3)No(er)2/3,  7(0, t)  =  (3/2)-N,e  nl  , 

where 

(A-18)  Nn  -  2 /“  donnG(n) , 

i.e.,  the  similarity  form  (3-1)  implies  the  requisite  limiting 


.  i 
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forms  (2-44)  and  (2-45)  independently  of  G(Q),  provided  only 

that  the  moments  N  and  N,  exist  and  do  not  vanish. 

o  I 

Expressions  (A- 17)  represent  the  leading  terms  of  the 
general  similarity  form  lf(r,T)  -  (er)^^F(z)  in  the  respective  limits 
z-K)  and  z-*®,  and  we  may  consider  the  z -dependence  of  the  next- to- 
leading  terms  in  these  limits.  The  dependence  as  z-K)  is  deter¬ 
mined  by  the  behavior  of  G(fl)  as  n-*®.  First,  suppose  G(n)  in 

.  O 

this  limit  decreases  faster  than  n  .  From  Eqs.  (2-15)  and  (3-1), 

by  expansion  of  l~cosu)T,  we  readily  find 

(A-19)  F(z)-(9/10)r(l/3)No-K3/2)r(2/3)N2z2  as  z-0, 

the  first  term  having  been  given  at  (A-17).  Now  suppose  instead 
G(n)-*gaBn~^2+B^  as  n-*®,  0<6<1 ,  g^  constant. 

In  this  instance  we  find 

(A-20)  F(z)^(9/10)r(l/3)No-fbogooz1+6  as  zK), 

b  =  2tt  r(26/3)cos(;,T6/3) 

°  1  -  26/3  r(2+6)cos(n6/2) 

Convergence  or  non-convergence  of  the  second  moment  N2  of  0(0)^ 
which  led  respectively  to  (A-19)  or  (A-20),  is  equivalent  to 
decrease  of  §(r,<ju) ,  the  frequency  t  ansform  of  the  proper-frame 

_3 

spatial  correlation,  more  or  less  rapidly  than  w  as 

-1/32/3  .  * 

c  r  u.'^30  • 

*We  may  consider  also  a  hypothetical  space-independent  velocity- 

wander  term  (A-8)  in  ?(r,x),  which  implies  a  next- to-leading  ten.. 

<*z  in  F(z)  as  z-K),  in  place  of  those  of  (A-19)  or  (A-20).  Th~. 

corresponding  contribution  to  E/(k,uu)  would  be  TT-iBc6  (k)x"  ,  thi- 

singular  function  is  the  limit  of  cne  of  the  form  (3-1),  formally 

corresponding  to  G(o)  «  -r_ ^BC"  '  ' ^5 (C^/  *) .  Such  an  H/(k,j  )  may  be 

regarded  as  the  limit  of  that  given  by  form  (3-4)  with  m  -  1  as 

v  -*0,  in  the  sense  that  TT-^x/(k2+x2)-»^ (k)as  x-*0. 
m 
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1 


Similarly,  the  dependence  of  F(z)  as  z is  related 

to  the  behavior  of  G(n)  as  Q-H).  In  particular,  if 

-1+A 

G(fi)-*gQn  as  Q-*0,  0<a<1,  gQ  constant, 

we  find 

(A-21)  Fcwi;i).v^il^  ggtfjffiSffi  V‘a 

as  z -►<*>,  the  first  term  having  been  given  at  (A-17).  The  condition 
A>0  is  required  for  finiteness  of  NQ  in  (A-17).  If  A  -  1/2, 
this  result  may  be  written 

(A-22)  F(z)-(3/2)ttN1z  +  1^/2) 1/2  gQz"1/2  as  z+®. 

The  type  of  function  assumed  for  G(Q)  at  Eq.  (3-4) 
meets  the  conditions  yielding  the  limiting  forms  for  F(z)  in  (A-19) 
and  (A-22)  above,  and  Eqs.  (3-5)  and  (3-6)  accordingly  represent 
particular  cases  of  these  results.  Apart  from  the  coefficient 
am  ,  form  (3-5)  is  thus  correct  provided  only  that  the  second 
moment  Nj  of  G(Q)  is  actually  finite.  Further,  the  form  (A-19) 

(and  (3-5))  agree  with  the  limiting  form  suggested  by  (A-5) 
and  (A-13),  or,  more  generally,  with  the  limiting  form  obtained 
if  differential  convection  determines  i|i  (r ,t)  (r ,0)  as  z  +  0. 

At  the  same  time,  in  the  opposite  limit  where  z  •*  ®,  if  differen¬ 
tial  convection  contributed  to  effective  displacement  in  quad¬ 
rature  with  the  contribution  from  velocity  wander,  as  in  (A-5) 
and  (A-13) ,  the  next-to-leading  term  in  F(z)  would  be  constant 
(A  -*  0  in  (A- 21)),  not  «z~^2  as  in  (A-22). 


113 


(A- 18)  E4(k,cu)  -  «1/3k_7/3G(ft)  ,  a  =  e_1/3k”2/3|aj | , 

in  the  domain  (3-3).  Form  (A-17)  is  the  limit  of  such  a  farm, 
formally  corresponding  to  G(ft)  *  Tr'^Bfi ~^26(ft"3/2)  ,  Form 
(A-18)  implies 

(A- 19)  J(r,0)  -  (9/10) T (l/3)No(sr) 2/3 ,  ^(0,t)  -  (3/2)ttN1«  |t|  , 

where 

(A -20)  Nn  -  2  ^°°cktt2nG(ft)  , 

i.e.,  the  similarity  form  implies  the  requisite  limiting  forms 

(2-44)  and  (2-45)  independently  of  G(fi)  ,  provided  only  that 

the  moments  N  and  N,  exist  and  do  not  vanish.  It  is  natural 
o  1 

to  examine  the  consequences  of  replacing  6(k)  in  (A-17)  by  a 

-1  2  2 

common  function,  tt  3/ (k  +6  )  ,  which  reduces  to  it  in  the  limit 

3  -►  0,  but  introduces  a  finite  scale  of  spatial  correlation  of 

1/2  -3/2 

appropriate  order  (~e  '  oj  ).  In  view  of  (A-18)  we  therefore 
assume 

(A-21)  (k ,oj)  =  7r'2Beco"2.7e“1/2ico|3/2/(k2-Ky2e"1|co|3)  , 

(A-22)  i.e.,  G(fi)  -  7r‘2B7ft"1/2  (l+y2^3)  _1. 

The  limit  (A-17)  would  correspond  to  ”>  ->  0.  Eq.  (A-21)  further 
yields,  by  use  of  (2-15),  the  limiting  forms  of  ^(r,x)  given 
at  (3-5)  and  (3-6).  Hence,  the  dependence  of  the  two  leading 
terms  of  f(r,T)  as  given  by  (A-21)  in  both  limits  z  -*•  0  and 
z  •*  oo  agrees  with  the  limiting  forms  suggested  by  (A -5)  and 
(A-13)  ,  or  more  generally  with  the  limiting  forms  obtained  if 
differential  convection  determines  tf  (r  ,t) -tf  (r ,0)  as  z  -*  0  and. 
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as  z  -*  oo,  contributes  to  effective  displacement  in  quadrature 

with  the  contribution  from  velocity  wander. 

We  may  inquire  how  general  is  the  dependence 
2 

F(z)-F(0)«cz  in  the  limit  z  -►  0 ,  implied  by  the  specific  form 

(A-21)  as  stated  at  (3-5),  ie.,  the  dependence  ^(r ,t) -^(r,0)*T^ 

1/3  -2/3 

as  e  7  r  7  t-*0.  The  singular  form  (A-17)  ,  for  example, 
corresponds  rather  to  *|t|.  In  general,  we  have  by  the 
inverse  of  the  proper-frame  analog  of  (2-26) 

A  =  2/°0dci)(l-cosarr)6(r,ao)  . 

T  o 

Hence,  provided  only  that  the  frequency  transform  of  the 

.3 

spatial  correlation  trace  falls  off  more  rapidly  than  m  for 

all  r  ±  0,  i.e.,  that  G(fi)  falls  off  more  rapidly  than  0“^, 

2  1/3  -2/3 

we  have  ill  fact  A^<*t  as  6  r  7  t  0,  since  then 

approaches  the  result  obtained  from  the  leading  term 

2 

1-cosrn  ->(1/2)(u}t)  ;  explicitly, 

(3/2)r(2/3)N2(£r)2/3(e1/3r"2/3x)2 


in  terms  of  the  assumed  convergent  moment  N2  defined  by  (A-20). 
On  the  other  hand,  if  0(r,oo)  falls  off  as  ,  0<  6<  1, 

i.e.  , 

G(fi)  gQ  as  fi  ^  «  (go  constant)  , 


we  obtain 


b0g0(er)2/3(£l/3r-2/3lT,)1+6 


*  I 


1-26/3  r (2+6) cos (tt6/2) 
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is  finite 


It  seems  likely  that  the  second  moment  ^  of  G(fl) 
through  we  cannot  exclude  the  contrary.  If  it  is,  then  assump 
tion  of  (A-21)  can  do  no  worse  than  yield  numerical 
constants,  e.g.  ,  CQ  in  (3-5),  whose  values  may  differ  somewhat 
from  the  correct  ones. 


4r 
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FIGURE  1 


FIGURE  2 


FIGURE  3 


Functions  ^(vjTj/r)  of  Eq.  (2-80).  These  functions 
express  the  fractional  departure  from  space-time 
isotropy  of  the  velocity  decorrelation  functions 
computed  in  the  local-convection  approximation  for 
the  inertial  subrange  [(1)  longitudinal,  (2)  trans¬ 
verse,  (0)  trace]. 

.  Space-time  correlations  for  streamwise  fluctuating 
velocity  in  laboratory  frame.  Experimental  curves 

- -  as  given  for  grid  turbulence  by  Favre 

(1965),  Figs.  1,5;  computed  curves  for  inertial  sub-  « 

range  by  local  convection  approximation,  Eq.  (2-90):  '  . 

- >  turbulence  intensity  vQ/uo  =  0.027, 

spatial  separations  aligned  with  mean  flow  (cp  =  0) ; 

- ,  vQ/uo  =  0.027  but  misalignment 

cp  =  1.5°,  or  ^  =  0  but  intensity  arbitrarily  adjusted 
to  VQ/U0  =  0.042.  Computed  curves  for  vanishing 
separation  coincide.  Results  are  omitted  at  correl¬ 
ations  <0.2. 


Space- time  correlations  for  streamwise  fluctuating 
velocity  in  laboratory  frame.  Experimental  curves 

- ,  as  in  Fig.  2,  computed  curves  for 

inertial  subrange  with  dispersion: - , 

envelope  from  non-dispersive  computation  ( - ) 

of  Fig.  2; - ,  vQ/uo  =  0.027,  cp  -  0,  disper¬ 
sion  included  by  Eqs .  (3-30),  (3-31)  with 

am  -  2.31;  — - - ,  vQ/u0  =  0.027,  cp  =  0,  dispersion 

included  by  space -independent  addition  to  decorrelations 
with  B/A0  =  1.04. 
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local-convection  approximation  for  the  inertial  subrange,  to  be  exactly 
valid  for  all  vQ/u0  provided  a  particular  effective  convection  velocity 

different  from  uq  is  assumed;  for  sufficiently  large  k^L,  where  L 

is  the  scale  of  energy -containing  eddies,  the  dispersive  correction 
to  this  result  is  negligible.  A  plausible  explicit  form  is  proposed 
for  the  intrinsic  energy  spectrum  in  the  co-moving  frame,  and  con¬ 
sequent  corrections  to  results  of  the  local-convection  approximation 
for  space-time  correlations  and  spectra  are  computed.  Tnis  dispersive 
correction  to  correlations  is  appreciable  at  Reynolds  numbers  of 
typical  grid-turbulence  experiments.  An  extension  of  the  basic  sep¬ 
aration  of  convection  and  the  local-convection  approximation  to  shear 
flow  is  suggested. 


From  kinematic  and  similarity  arguments  for  the  inertial  sub¬ 
range,  inferences  are  also  made  concerning  pressure  spectra  in  a 
measurement  frame  with  velocity  (-u  ,0,0)  relative  to  the  unsheared 
flow.  A  similarity  form  results  foP  the  wavenumber -frequency  spectrum 


F(K,u>)[K  =  (k 

0“  -  ki%l  < 
obtained  for 


p  ky ]  in  the  vicinity  of  the  convective  ridge 
vqK) .  For  vq/uo  «  1,  the  functional  form  is  then 
the  point  frequency  spectrum  P(ai)(o«^  ;  the  convec¬ 


tive  contribution  to  the  spectrum  Q  (u>)  of  average  pressure  on  a 

moving  circular  area  of  radius  R  for  u)R  /u  »1  is  smaller  by  a 

o  o  o  J 

factor  0.8(uoR  /u  )  "  .  On  assumption  of  quasinormality  of  the  velocity 
distribution  8nd°use  of  the  non-dispersive  approximation  of  space-time 

isotropy  for  the  longitudinal  correlation,  the  spectrum  P(K,^')  in  the 
vicinity  of  the  convective  ridge  is  determined  explicitly,  as  well 
as  the  cross-spectral  density  and  space-time  correlation  of  pressure. 

From  quasinormality,  P(K,U))  is  estimated  also  in  the  low -wavenumber 

domain.  Ia  the  limit  <a>R  /v  »1  the  ratio  of  the  low-wavenumber  to 

o  o  a 

the  mean-convective  contribution  to  the  moving-arta  spectrtim  Q(u-)  is 
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l>  ABSTRACT 

Kolmogorov's  principles  provide  a  basis  for  the  treatment  of 
iulexian  space-time  correlations  for  turbulence  in  the  universal 
range  by  explicit  separation  of  the  kinematic  effect  of  convection 
y  large  eddies.  With  reference  to  unsheared  homogeneous  turbulence 
he  usual  similarity  forms  are  assumed  for  intrinsic  velocity  cor¬ 
elations  in  a  local  co-moving  frame.  The  local -convection  approxi- 
ation,  neglecting  dispersion  in  this  frame  and  hence  relating 
pace-time  correlations  to  pure  spatial  correlations,  is  indicated 
to  be  a  useful  one  in  the  inertial  and  viscous  subranges.  For  an 
isotropic  normal  velocity  distribution,  the  structure  functions  of 
fluctuating  velocity  are  computed  and  found  to  be  nearly  space-time 
isotropic  in  the  former  subrange  and  exactly  so  in  the  latter.  The 
wavenumber -frequency  spectrum  of  energy  in  the  mean  rest  frame, 
E^(k,^) ,  however ,  in  the  regime  of  large  u)/vQk,  where  v  denotes  rms 

velocity,  is  not  given  correctly  by  the  local-convection  approxi¬ 
mation,  but  essentially  involves  dispersion.  Taylor's  hypothesis 
relating  wavenumber  (k,)  spectra  in  the  mean  rest  frame  to  frequency 
spectra  in  a  measurement  frame  with  velocity  -uQ  is  found,  in  the 
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p.  15,  line  preceding  Eq.  (2-14);  For  "energy  of"  read  "energy  by." 
p.  18,  line  (-8):  For  "P(v)v2dv"  read  "P(v)v2dv." 

p.  21:  In  Eq.  (2-38),  in  the  definitions  of  and  £.,  and  in  the  line 
following,  for  "v"  read  "v." 

p.  25,  second  line  before  Eq.  (2-52):  For  "v"  read  "v." 

p.  30,  (2-66):  For  "v^2  read  "v^2. 

p.  37,  line  (-9):  For  "relation  (2-64)"  read  "relation  (2-59)." 

p.  45,  line  6:  Insert  comma  after  "u  ." 

p.  45,  Eq.  (2-100):  For  "0"  read  "0." 

p.  46,  Eqs .  (2-104),  (2-105):  For  "P(v)M  read  "P(v)." 

p.  48,  line  2:  For  "for"  read  "For." 

p.  48,  lines  9  and  11:  For  "r"  read  "r"  and  for  "r^"  read  "r^"  in 
arguments . 

p.  49,  line  (-3):  For  "0^"  read  "0^." 
p.  49,  line  8:  For  "u^"  read  "u^." 
p.  50,  Eq.  (2-119):  Insert  (  after  [  in  second  line, 
p.  56,  line  12:  Close  parentheses  after  "respectively." 
p.  60,  line  (-5):  Close  parentheses  after  "(3-6)." 
p.  61,  line  (-5):  Close  parentheses  in  T(l/3)." 
p.  62,  line  (-10):  For  "7"  read  "7  ." 

p.  63,  Eq.  (3-18),  and  p.  64,  Eq.  (3-20):  Move  factor  "k2^m  ^ 

00 

ft  r  on  o  »» 

J  dk" ;  for  "y*  read  "7 
o  m 

p.  63,  Eq.  (3-18):  For  "0"  read  "0." 

p.  64,  Eq.  (3-21):  For  "b  (B/A  )2"  read  "b^  < 

m  o  o  m 

p.  64,  line  (-4):  For  "b^  and  b^  were"  read  "bQ  wa3." 
p.  65,  line  6:  For  "AT"  read  "AT  ." 
p.  68,  line  (-7):  For  "eo"  read  "S0*" 
p.  72,  Eq.  (4-5):  For  "u^K"  read  "uo*K." 
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p.  73,  line  6:  Read  "P(oj)  °c  <o 
p.  80,  line  1:  For  "k  «  1"  read  "k  «  L~*. 
p.  80,  line  3:  For  "F(K,<o)"  read  "P(K,u>)." 
p.  84,  line  12:  For  "K1"  read  "k^" 
p.  86,  line  7:  Delete  +  before  k,  . 

~  lm  ®  , 

p.  87,  Eq.  (4-52):  Should  read  "W  (£,t)  -  /  doa  e~lun  8  (C,ca)." 

p  -oo  p 

p.  90,  line  2:  Read  "k^  ^  <Vvq(»  K)." 

p.  94:  Insert  comma  after  'V.  line  10,  and  after  "frame," 

line  11.  1J 

p.  106,  line  (-2):  For  "through"  rend  "though." 
pp.  114-116:  Delete  these  pages. 

p.  117,  seventh  reference:  Following  "Bolt"  add  "Beranek  &  Newman 
Rept.  No.  1310." 


